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Abstract 

This  study  examined  three  aspects  of  animation  (movement,  color,  and  detail  of 
icons)  to  determine  which  one  (or  ones)  best  communicated  the  operation  of  a  simulation 
model.  The  procedure  was  done  in  the  context  of  using  animation  to  establish  a  model’s 
face  validity.  Movement,  color,  and  detail  of  icons  were  looked  at  individually  and  in 
combination.  The  ability  to  communicate  was  measured  both  subjectively  and  objectively. 
The  subjective  measures  were  a  selection  of  “best’*  and  “worst”  animation  types  where 
“best”  and  “worst”  referred  to  how  well  an  animation  communicated,  and  a  pairwise 
comparison  of  the  animation  types  which  resulted  in  preference  ratings  for  each  animation. 
There  were  seven  different  scenarios  containing  various  problems  with  the  system.  The 
objective  measures  were  subject  problem  identification  accuracy  and  time  delay  of  problem 
identification.  The  results  showed  that  movement  in  animations  was  always  preferred  to 
a  lack  of  movement  in  animations.  However,  movement,  color,  and  detail  of  icons  in 
combination  was  preferred  the  most.  Objectively,  movement  was  the  most  important 
aspect.  The  subjects  performed  equally  well  for  all  the  animations  with  movement  and, 
when  there  was  no  movement,  the  subjects’  performance  dropped  equally. 


USING  ANIMATION  IN  THE  VALIDATION  OF 
SIMULATION  MODELS 


I.  Introduction 

Background 

Over  the  past  two  decades  simulation  modeling  in  the  military  has  increased  sig¬ 
nificantly.  Today  the  military  uses  simulation  models  in  many  areas:  defense  planning, 
weapon  system  evaluation,  logistics,  and  the  like.  Models  have  become  essential  in  mil¬ 
itary  analysis  because  a  valid  simulation  model  can  point  out  potential  trouble  spots  or 
potential  solutionis  in  a  system  without  the  expense  in  time  and  money  of  operating  the 
actual  system.  Also,  simulation  models  are  valuable  tools  for  planning  and  evaluating 
future  systems. 


In  general,  models  are  abstractions  of  reality  that  attempt  to  capture  the  essence  of 
the  problems,  events,  or  actions  bong  considered.  Thus,  no  effort  is  made  to  include  every 
detail  in  a  model.  In  the  first  place  that  would  be  impossible,  and  secondly,  including  every 
detail  is  unnecessary.  Developing  a  model  is  like  an  artist  painting  a  portrait.  The  artist 
initially  draws  a  sketch  of  what  is  to  be  painted.  The  sketch,  like  the  model,  contains  the 
essentials  but  omits  the  details.  Just  as  the  sketch  cannot  be  a  substitute  for  the  portrait, 
a  model  cannot  be  a  substitute  for  reality.  However,  a  sketch  well  done  leads  to  a  better 
portrait,  and  an  accurate  model  leads  to  a  better  understanding  of  the  literal  or  proposed 
situation.  Models  are  developed  in  various  forms  (e.g.  mathematical,  verbal,  or  pictorial). 
Simulation  models  are  generally  considered  to  be  models  that  are  developed  as  computer 
programs. 

Although  simulation  models  are  used  extensively  in  decision  making  and  problem 
solving,  many  decision  makers  lack  confidence  in  simulation  model  results  (Sargent,  1991:37). 
If  models  are  to  be  legitimate  decision  aids,  it  is  critical  that  they  be  validated.  Simulation 
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model  validation  is  “substantiation  that  a  computerized  model  within  its  domain  of  appli¬ 
cability  possesses  a  satisfactory  range  of  accuracy  consistent  with  the  intended  application 
of  the  model”  (Schlesinger  and  others,  1979:103).  The  validation  process  consists  of  per¬ 
forming  tests  and  evaluations  during  the  development  of  a  simulation  model.  Thece  tests 
range  from  a  review  of  model  assumptions  to  detailed  statistical  procedures.  Although 
model  verification  is  closely  linked  with  model  validation,  the  two  should  not  be  confused. 
Simulation  model  verification  Is  the  process  of  ensuring  that  the  actual  computer  program 
is  error  free.  As  Whitner  and  Balci  state,  “Validation  deals  with  building  the  right  model, 
verification  deals  with  building  the  model  right”  (Whitner  and  Bald,  1989:559). 

Even  though  much  has  been  written  concerning  simulation  model  validation,  and 
many  validation  techniques  have  been  developed,  no  simulation  model  can  ever  be  fully 
validated.  A  successful  result  from  a  validation  procedure  can  only  increase  the  validity 
of  a  model.  Thus,  the  modeler  desires  to  tip  the  scales  of  validation  towards  complete 
validity  as  far  as  possible,  all  the  while  knowing  that  there  will  always  be  some  uncertainty 
weighing  down  the  other  side. 

Computer  animation  is  one  of  many  techniques  used  in  the  process  of  simulation 
model  validation  (and  verification).  Through  computer  animation  a  “model’s  operational 
behavior  is  displayed  graphically  as  the  model  moves  through  time”  (Sargent,  1991:39). 
Animation  is  becoming  more  popular  because  animation  software  has  dropped  In  price 
and  increased  in  quality.  Increased  computer  graphics  capability  allows  the  modeler  to  see 
the  simulated  operation  of  the  system  instead  of  just  a  printout  of  statistics.  In  addition, 
animation  can  be  used  to  enhance  a  model’s  credibility.  A  credible  model  is  a  model  the 
decision  makers  are  willing  to  use  because  they  have  confidence  in  the  model’s  results  (Sar¬ 
gent,  1991:37).  According  to  Law  and  Kelton,  the  ability  to  increase  a  model’s  credibility 
is  the  main  reason  for  animation’s  expanding  use  (Law  and  Kelton,  1991:241).  Thus,  the 
graphical  display  of  a  model  through  animation  can  add  credence  to  a  model  and  increase 
the  confidence  in  a  model. 


Problem 

The  use  of  computer  animation  in  simulation  model  validation  is  increasing;  how¬ 
ever,  there  are  few  guidelines  for  using  animation  in  conjunction  with  accepted  simulation 
validation  techniques.  One  primary  method  used  in  model  validation  is  face  validity.  “Face 
validity  is  asking  people  knowledgeable  about  the  system  whether  the  model  and/or  its  be¬ 
havior  is  reasonable”  (Sargent,  1991:39).  In  order  to  establish  the  face  validity  of  a  model, 
the  assumptions  and  operation  of  the  model  must  be  communicated  to  the  system  expert. 
Animation  is  one  tool  for  accomplishing  that.  Thus,  this  research  focused  on  determining 
which  aspects  of  animation  are  the  most  useful  for  communicating  the  operation  of  the 
model. 

Assumptions  and  Scope 

Simulation  Model.  The  simulation  model  used  for  this  research  was  taken  from  a 
SLAM  (Simulation  Language  for  Alternative  Modeling)  textbook.  A  simple  model  was 
chosen  so  that  some  of  the  basics  of  animation  could  be  examined  without  the  problems 
associated  with  a  complex  model.  The  simulation  models  one  bulldozer,  four  trucks,  and 
two  man-machine  loaders.  The  bulldozer  stockpiles  material  for  the  loaders.  The  loaders 
place  the  material  into  the  trucks.  Once  a  truck  is  loaded  it  hauls  the  material  away, 
dumps  the  load,  and  then  returns  for  another  load.  The  Truck  Hauling  Situation  model 
will  hereafter  be  referred  to  as  the  Loader  model  (Pritsker,  1986:237-242). 

Animation  Software.  The  model  was  animated  with  Proof  Animation.  Proof  Anima¬ 
tion  (hereafter  referred  to  as  Proof)  is  a  PC-based,  “post-processing”  animation  software. 
Post-processing  (or  playback)  means  the  animation  is  seen  after  the  simulation  is  run.  The 
events  or  state  changes  are  recorded  in  a  file  during  the  simulation  ran  and  then  “played 
back”  by  the  animator.  This  is  in  contrast  to  “concurrent”  animation.  Concurrent  means 
the  animation  is  seen  while  the  simulation  is  running  (Law  and  Kelton,  1991:241).  There 
are  advantages  and  disadvantages  to  each  type  of  animation.  Concurrent  animation  usu¬ 
ally  allows  the  analyst  or  user  to  immediately  see  the  result  of  changes  they  make  while 
the  simulation  is  running;  however,  certain  changes  cannot  be  made  while  the  simulation 


is  miming.  Also,  the  end  (or  any  other  part)  of  a  simulation  ran  can  only  be  seen  after 
viewing  all  that  comes  before.  That  is,  “fast-forwarding”  is  not  possible.  Post-processing, 
on  tb»  other  hand,  does  allow  one  to  view  to  any  part  of  a  simulation  run  without  waiting 
through  the  previous  portions  of  the  simulation,  but  changes  cannot  be  made  and  seen 
interactively  (Brunner  and  others,  1991:91). 

Scope  of  Animation.  Three  aspects  of  animation  were  considered: 

•  Movement 

•  Detail  of  Icons 

•  Color 

Each  of  the  above  was  looked  at  individually  and  in  combination.  Other  factors  that  could 
be  considered  include  graphs,  perspective  (as  in  two-dimensional  or  three-dimensional), 
concurrent  animation  versus  playback,  and  speed  of  animation. 

Approach 

The  main  reason  for  using  animation  is  that  it  improves  communication.  Thus,  the 
procedure  was  to  animate  the  Loader  model,  and  show  the  animations  to  volunteers.  They 
determined  how  useful  the  animations  were  in  enhancing  face  validity  and  which  animation 
type  was  most  beneficial.  This  was  measured  in  two  ways.  First,  each  animation  showed 
a  different  scenario.  The  subject  attempted  to  determine  if  there  were  any  efficiency 
problems  with  the  system  modeled  (queue  buildups,  long  idle  times,  and  the  like).  The 
time  delay  in  identifying  a  potential  problem  and  the  identification  error  rate  were  recorded. 
Second,  after  viewing  the  animations,  subjects  performed  a  pairwise  comparison  of  each 
animation’s  ability  to  communicate.  A  rating  for  each  animation  was  calculated  from  the 
pairwise  comparisons.  The  pairwise  comparison  procedure  is  called  the  Analytic  Hierarchy 
Process  and  is  discussed  in  the  next  chapter. 


Overview 


The  next  chapter  reviews  the  current  literature  concerning  simulation  model  valida¬ 
tion  and  animation.  Also,  a  section  discussing  the  Analytic  Hierarchy  Process  is  included. 
Chapter  3  gives  a  detailed  description  of  the  approach  to  the  research,  the  Loader  model, 
the  animations,  the  experimental  design  and  data  collection,  and  the  analysis  methods 
used.  The  fourth  chapter  is  a  discussion  of  the  results,  and  the  final  chapter  contains 
conclusions  and  recommendations.  The  data  collection  forms  are  shown  in  Appendix  A. 
Appendix  B  is  the  SAS  output  used  in  the  analysis,  and  Appendix  C  contains  the  SLAM 
and  Fortran  code. 


11.  Literature  Review 
Scope  end  Method  of  Presentation 

This  review  will  cover  information  about  validation  of  simulation  models,  the  use  of 
animation  in  validation,  and  the  Analytic  Hierarchy  Process  (AHP).  The  discussion  will 
begin  with  model  validation  in  general  and  include  animation’s  applicability  to  validation. 
A  review  of  information  dealing  with  the  use  of  animation  in  validating  models  will  follow, 
and  a  look  at  the  AHP  will  conclude  the  discussion. 

Discussion 

Model  Validation.  Since  the  1960’s  much  has  been  written  concerning  model  vali¬ 
dation.  Bald  and  Sargent  list  over  300  artides  in  a  comprehensive  bibliography  (through 
1983)  of  artides  dealing  with  validation  and  credibility  assessment  of  models  (Bald  and 
Sargent,  1984:15-27).  Carson  defines  validation  and  a  credible  model  as  follows: 

•  Validation:  “The  process  whereby  the  modeler  and  end  user  ask  the  questions:  How 
accurately  does  the  model  represent  reality?  Can  the  model  be  used  in  place  of  the 
real  system  for  the  purpose  of  making  decisions  concerning  the  real  system*  (Carson, 
1989:552)? 

•  Credible  Model:  A  model  “that  is  accepted  by  the  client  as  being  sufficiently  accurate 
to  be  used  as  an  aid  in  making  decisions”  (Carson,  1989:552). 

One  of  the  most  published  authors  in  the  area  of  validation  is  Sargent.  He  emphasizes  the 
relationship  between  the  purpose  of  a  model  and  the  validity  of  a  model.  A  model’s  validity 
can  only  be  determined  in  relation  to  the  purpose  of  a  model.  If  a  model  is  designed  to 
evaluate  several  different  problems,  it  must  be  validated  for  each  of  the  problems  consid¬ 
ered.  Also,  a  model  can  never  be  considered  completely  valid.  He  states,  “It  is  often  too 
costly  and  time  consuming  to  determine  that  a  model  is  absolutely  valid  over  the  complete 
domain  of  its  intended  applicability*  (Sargent,  1991:37).  There  are  only  degrees  or  levels 
of  validity  (Sargent,  1991:37). 
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In  his  discussion  of  validation  techniques  Sargent  distinguishes  between  subjective 
and  objective  (statistical)  validation  techniques.  Sargent  lists  15  subjective  techniques,  one 
being  animation.  He  states  that  the  model  validation  process  is  included  in  the  process 
of  model  development  and  is  concerned  with  four  different  evaluations:  conceptual  model 
validity,  computerized  model  verification,  operational  validity,  and  data  validity.  All  of  the 
validation  techniques  he  lists  are  useful  in  determining  operational  validity.  Operational 
validity  ensures  “the  model’s  output  behavior  has  the  accuracy  required  for  the  model’s 
intended  purpose  over  the  domain  of  its  intended  applicability”  (Sargent,  1991:39).  Most 
validation  occurs  during  operational  validity  testing  because  any  problems  found  in  the 
computerized  model  could  be  due  to  a  faulty  conceptual  model,  improper  programming, 
or  incorrect  data.  Thus,  operational  validity,  in  a  sense,  encompasses  conceptual  model 
validity,  computerized  model  verification,  and  data  validity  (Sargent,  1991:39-42). 

Bald  also  divides  validation  techniques  into  two  categories:  subjective  and  statisti¬ 
cal.  He  lists  13  subjective  validation  techniques  and  18  statistical  validation  techniques 
(Bald,  1989:67-68).  Animation  is  not  listed  as  a  validation  technique.  In  fact,  Bald  does 
not  mention  animation  at  all.  Banks  divides  validation  techniques  into  subjective  and 
statistical,  too  (Banks,  1989:550-551).  However,  he  does  not  list  animation  as  a  validation 
technique  either. 

A  three-step  approach  for  validation  is  offered  by  Carson.  He  places  face  validity 
at  the  beginning  of  his  validation  process.  Face  validity  is  a  process  which  involves  ques¬ 
tioning  those  knowledgeable  about  the  system  modeled  to  determine  whether  a  model  or 
its  operation  is  consistent  with  the  system  modeled.  Carson  then  recommends  validation 
of  the  model  assumptions  and  model  output.  In  contrast  to  Bald  and  Banks,  though, 
Carson  emphasizes  the  need  for  animation.  He  condudes  with  a  listing  of  verification  and 
validation  techniques  (Carson,  1989:555,557). 

Animation.  There  are  many  proponents  for  using  animation,  not  only  for  model 
validation,  but  also  for  the  whole  model  building  process.  Those  who  tout  animation  the 
most  are  those  who  sell  animation  software  (Brunner  and  others,  1991:90-94;  Kalasky  and 
Davis,  1991:123;  HoQocks,  1984:322-328;  Standridge,  1986, 121-143).  Users  of  animation 
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are  the  second  most  enthusiastic  proponents  of  animation  (Aiken  and  others,  1390:775-783; 
Johnson  and  Poorte,  1988:30-36;  Carson  and  Atala,  1990:798-801).  Academia,  although 
acknowledging  value  in  animation,  generally  stresses  other  methods  of  validation  or  does 
not  explicitly  mention  animation  at  all.  (Law  and  Kelton,  1991:242;  Bald,  1989:62-71; 
.  Banks,  1989:549-551) 

As  an  example  of  the  latter,  Law  and  Kelton  recognize  that  animation  is  increasing 
the  popularity  of  simulation  modeling.  However,  with  regard  to  model  validation,  they 
point  out  that  animation  can  show  that  a  model  is  not  valid,  but  it  cannot  show  that  a 
model  is  valid.  Also,  animation  cannot  replace  statistical  analysis  of  the  model  output. 
Another  drawback  is  that  animation  usually  cannot  be  watched  for  the  complete  run  of 
the  simulation.  Therefore,  an  error  could  occur  at  a  point  later  in  the  simulation  run  when 
there  is  no  one  watching.  A  final  drawback  is  “only  part  of  a  simulation  model’s  logic  can 
actually  be  seen  in  an  animation;  thus,  a  ‘correct’  animation  is  no  guarantee  of  a  valid  or 
debugged  model”  (Law  and  Kelton,  1991:242). 

I 

Animation  can  improve  face  validity  according  to  Schuppe.  However,  developing 

quality  animation  can  take  as  long  as  developing  the  model,  and  the  animation  should  be 

I 

verified.  That  is,  the  modeler  should  substantiate  that  the  animation  is  not  giving  a  false 
picture  of  what  the  simulation  model  is  doing  (Schuppe,  1991:523, 525). 

Canon  paints  a  much  brighter  picture  for  animation.  He  states  that  animation 
provides  a  degree  of  validity  that  previously  could  not  be  reached  for  simulations  involving 
items  travelling  through  space.  With  an  animator  spatial  relationships  an  easier  to  verify. 
He  recognizes  there  is  a  possibility  judgments  will  be  made  based  on  animation  before  the 
validation  process  is  complete.  However,  Canon  asserts  the  probability  of  this  happening  is 
not  high  because  those  who  use  simulations  “will  want  to  see  the  output  data  of  numerous 
long  runs  before  making  any  final  judgments”  (Carson,  1989:555). 

According  to  Johnson  and  Poorte,  animation  can  be  useful  for  debugging  and  verifi¬ 
cation,  validation,  analysis,  and  communication  and  presentation.  They  offer  a  hierarchical 
approach  for  using  animation  in  all  of  these  areas.  With  regard  to  validation,  they  believe 
“animation  can  provide  a  vital  link  between  the  modeler  and  expert”  (Johnson  and  Poorte, 
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1988:31).  Usually  the  system  expert  involved  in  conceptual  model  validation  is  unfamiliar 
with  simulation  code;  thus,  animation  improves  the  communication  between  the  modeler 
and  the  expert,  and  conceptual  model  validity  is  enhanced.  Nevertheless,  they  acknowl¬ 
edge  “animation  is  a  supplementary  tool  and  should  not  be  viewed  as  a  replacement  for 
standard  analysis  techniques*  (Johnson  and  Poorte,  1988:31,36). 

Cyr  states  that  animations  of  simulated  activities  offers  several  benefits: 

•  Animation  can  help  explain  simulations  (in  particular,  Monte  Carlo  simulations)  to 
upper-level  managers,  who  are  not  experts  in  mathematics  and  statistics. 

e  Animation  can  help  the  manager  and  others  understand  the  system  being  modeled. 

•  Animation  shows  actual  or  potential  system  problems  that  might  be  hard  to  under¬ 
stand  or  might  be  missed  by  looking  at  lists  of  numbers. 

•  Animation  opens  up  the  opportunity  to  interact  with  the  simulation. 

e  Animation  lowers  the  development  costs  of  models  by  uncovering  problems  that 
would  be  hard  to  see  without  animation. 

He  adds  that  the  value  of  a  simulation  models  increases  when  animation  is  added,  and  the 
additional  cost  is  outweighed  by  the  advantages  of  animation  (Cyr,  1992:1000, 1002). 

As  those  involved  with  the  making  and  marketing  of  animation  software,  Kalasky 
and  Davis  boldly  declare, 

In  recent  years  animation  has  become  a  requirement  of  the  simulation  process. 

One  of  the  reasons  for  this  requirement  is  that  numeric  summary  statistics  do 
not  necessarily  convey  information  about  the  dynamic  interactions  of  compo¬ 
nents  of  a  system.  Although  summary  statistics  are  a  crucial  part  of  evaluating 
the  performance  of  a  simulated  system,  it  is  only  through  animation  that  the 
analyst  can  easSy  identify  the  system  status  under  which,  for  example,  bottle¬ 
necks  occur.  (Kalasky  and  Davis,  1991:122) 

Concerning  validation,  they  agree  with  Johnson  and  Poorte  that  animation  improves  com¬ 
munication  between  the  system  experts  and  the  modelers.  In  addition,  they  state  that  sim¬ 
plifying  assumptions  can  be  seen  easier  with  animation.  (See  also  Standridge,  1986:121.) 
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lb  conclude  their  listing  of  the  uses  of  animation  they  give  the  standard  qualifier:  “anima¬ 
tion  cannot  replace  standard  statistical  analysis  techniques”  (Kalasky  and  Davis,  1991:122- 
123). 

Analytic  Hierarchy  Process.  The  AHP  is  a  decision  making  aid.  It  allows  a  subject 
to  subjectively  compare  alternatives.  A  subject  compares  every  alternative  with  every 
other  possible  alternative.  Thus,  the  AHP  is  a  pairwise  comparison.  Ratings  for  each 
alternative  can  be  calculated  from  the  pairwise  comparison  data.  This  is  in  contrast  to 
absolute  estimation  methods  (such  as  an  ordinal  ranking  of  alternatives).  For  example, 
the  pairwise  comparison  of  the  animations  allowed  a  subject  to  rate  a  specific  animation’s 
ability  to  communicate  the  model’s  operation  against  the  ability  of  every  other  animation. 
This  permitted  the  relative  merits  of  the  animations  to  be  quantified.  Since  there  were 
seven  animations,  the  subjects  performed  21  pairwise  comparisons  each.  (There  were  not 
28  comparisons  because  an  alternative  is  not  compared  with  itself.) 

The  AHP  was  developed  by  Thomas  Saaty  (Saaty,  1980).  The  advantages  of  the 
AHP,  as  summarized  by  Vidulich  and  Tsang  (1987),  follow. 

In  comparing  the  AHP  approach  to  the  absolute  estimation  methods,  Saaty 
(1980)  suggested  that  the  AHP  has  the  following  advantages:  (1)  While  the 
number  of  decisions  are  more  numerous  ([n(n  -  l)/2]  as  opposed  to  [n]  for 
n  conditions  to  evaluate),  they  are  simpler  because  the  subject  can  focus  on 
the  relationship  between  only  two  conditions.  (2)  The  comparison  of  each 
condition  to  every  other  prorides  a  great  deal  of  redundant  information  to 
improve  reliability.  (3)  It  is  possible  to  calculate  the  “consistency”  of  each 
subject’s  pairwise  mi  .ix  and  thereby  test  the  subject’s  ability  to  make  the 
discriminations  necessary.  (Vidulich  and  Tsang,  1987:1058) 

Given  n  alternatives,  the  pairwise  matrix  (or  judgment  matrix)  is  an  n  x  n  matrix  with 
each  row  and  column  headed  by  the  choices.  The  numerical  rating  from  each  pairwise 
comparison  is  placed  in  the  matrix,  and  the  principal  eigenvector  for  the  matrix  is  calcu¬ 
lated,  which  gives  the  rating  for  each  alternative  (Vidulich  and  Tsang,  1987:1058).  There 
are  three  steps  to  using  judgment  matrices:  collect  the  judgment  data,  construct  the  judg¬ 
ment  matrices,  and  calculate  the  ratings.  Vidulich  provides  a  good  overview  of  this  process 
(Vidulich,  1989:1407).  ,  7‘ 


Williair  s  and  Crawford  offer  another  method  of  calculating  ratings.  They  suggest 
using  a  geometric  means  approach  instead  of  an  eigenvector  approach  (Williams  and  Craw¬ 
ford,  1985:1).  They  state, 


The  geometric  mean  vector  v  =  v1,va,...,vn,  given  by 


which  satisfies  the  continuity  and  consistency  criteria  Saaty  uses  to  defend  the 
dominant  eigenvector,  has  several  other  desirable  traits:  In  certain  circum¬ 
stances,  it  is  statistically  optimal  and  gives  rise  to  an  estimate  of  scales  and  a 
measure  of  consistency  that  have  known  statistical  distributions.  In  empirical 
studies  reported  here  it  seems  to  do  as  well  as,  or  better  than,  the  eigenvector 
in  preserving  rank  order.  In  addition,  it  is  supported  by  a  literature  describing 
methods  of  handling  a  wealth  of  variations  of  the  problem,  including  missing 
data  and  multiple  judges.  (Williams  and  Crawford,  1985:vi) 


Another  issue  is  whether  judgment  matrices  are  consistent.  Vidulich  defines  a  consis¬ 
tent  matrix  as  one  where  there  are  “transitive  trends  among  related  judgments”  (Vidulich, 
1989:1407).  For  instance  if  alternative  A  is  preferred  twice  as  much  as  alternative  B, 
and  alternative  B  is  preferred  three  times  as  much  as  alternative  C,  then  alternative  A 
should  be  preferred  six  times  as  much  as  alternative  C.  If  these  types  of  relationships  hold 
within  the  matrix,  the  matrix  is  considered  consistent.  Inconsistency  increases  as  these 
relationships  are  violated  (Vidulich,  1989:1407).  Budescu  and  others  give  5*,  the  measure 
of  consistency  developed  by  Williams  and  Crawford. 
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where  is  the  pairwise  comparison  ratio  for  the  ith  and  y'th  alternatives  and  GMi  and 
GMj  are  the  geometric  means  of  rows  t  and  j  respectively  (Budescu  and  others,  1986:71). 


Bndescu  and  others  developed  a  consistency  criterion  for  the  geometric  means  ap¬ 
proach  during  their  comparison  of  the  two  methods  (Budescu  and  others,  1986:69-78).  In 
his  own  comparison  the  eigenvector  approach  and  geometric  mean  approach,  Vidulich  die- 


cusses  the  rale  developed  by  Budescu  and  others  and  reviews  the  comparison  of  the  two 
approaches  performed  by  Williams  and  Crawford  and  Budescu  and  others.  He  states 


Williams  and  Crawford  (1980)  did  a  Monte  Carlo  comparison  of  the  eigenvec¬ 
tor  and  geometric  mean  approaches  to  judgment  matrix  analysis.  They  used 
matrices  that  were  selected  from  known  distributions  with  varying  degrees  of 
error  perturbations  introduced.  Williams  and  Crawford  found  that  the  two  ap¬ 
proaches  gave  similar  results  when  perturbations  away  from  consistency  were 
minimal.  But  for  large  perturbations,  the  geometric  means  approach  produced 
resuits  that  deviated  less  from  the  known  distribution.  The  relative  advantage 
of  the  geometric  means  procedure  increased  with  both  matrix  size  and  error 
variability.  (Vidulich, 1989:1408) 

Matrices  of  randomly  generated  data  were  used  by  Budescu  and  others  to  compare 
the  two  means  of  calculating  ratings.  Budescu  and  others  showed  that  5 3  remained  stable 
as  matrix  size  increased.  Vidulich  suggested  at  the  end  of  his  study  that  the  geometric  mean 
approach  can  be  used  instead  of  the  eigenvector  approach.  Also,  Vidulich  recommended 
that  when  using  S2  as  &  consistency  criterion,  the  matrix  size  should  range  from  6  x  6  to 
10x10  (Vidulich,  1989:1410). 

Conclusion 

Many  techniques  are  offered  for  model  validation.  These  techniques  are  generally 
divided  into  subjective  and  objective  (statistical)  techniques.  Animation,  when  listed  for 
use  in  model  validation,  is  considered  to  be  a  subjective  technique.  Even  though  animation 
is  bring  used  in  simulation  modeling  more  and  more,  the  accepted  range  of  animation’s 
applicability  to  simulation  model  validation  varies  from  helpful  to  almost  indispensable. 
However,  animation  can  never  replace  statistical  validation  techniques. 

The  Analytic  Hierarchy  Process  is  a  useful  means  of  comparing  alternatives.  Al¬ 
though  there  is  debate  concerning  the  method  of  calculating  ratings  (eigenvector  versus 
geometric  means),  the  geometric  means  approach  is  easier  to  use  and  has  a  more  stable 
measure  of  inconsistency. 


III.  Methodology 


\ 


Approach 

The  experiment  consisted  of  subjects  viewing  seven  different  types  of  animations.  The 
subjects  determined  how  useful  the  animations  were  in  communicating  model  operation 
and  which  animation  type  was  most  beneficial.  This  was  measured  in  two  ways.  First, 
each  animation  showed  a  different  scenario.  The  subject  attempted  to  determine  if  there 
were  any  efficiency  problems  with  the  system  modeled  (queue  buildups,  long  idle  times, 
and  the  Qke).  The  time  delay  in  identifying  a  potential  problem  and  the  correctness  of 
problem  identification  were  recorded.  Second,  after  viewing  the  animations,  the  subject 
performed  a  pairwise  comparison  of  each  animation’s  contribution  to  face  validity.  A  rating 
for  each  animation  was  calculated  from  the  pairwise  comparisons  (AHP)  which  measured 
the  subject's  preference  for  each  animation  type.  Preparation  for  the  experiment  included 
selecting  a  simulation  model,  animating  the  model,  creating  system  problems  (different 
scenarios),  constructing  the  experimental  design,  and  preparing  data  collection  forms.  The 
next  step  was  to  perform  the  experiment  and  collect  data  by  showing  the  animations  to 
AFIT  faculty  and  student  volunteers.  The  final  step  was  analyzing  the  data. 

Initialization 

Simulation  Model  The  Loader  model  is  a  simple  SLAM  network  simulation  that 
models  a  loading  and  hauling  operation  for  480  minutes  (8  hours).  The  system  modeled 
consists  of  one  bulldozer,  four  trucks,  and  two  man- machine  loaders.  The  bulldozer  stock¬ 
piles  material  for  the  loaders.  Two  piles  of  material  must  be  stocked  prior  to  the  initiation 
of  any  load  operation.  In  addition  to  the  two  loads  of  material,  a  loader  and  an  unloaded 
truck  must  be  available  before  the  loading  operation  can  begin.  The  time  to  bulldoze  a  load 
is  Erlang  distributed  and  is  the  sum  of  two  exponentials  each  with  a  mean  of  4  minutes. 
The  loaders  are  modeled  as  servers  with  loading  time  for  server  1  exponentially  distributed 
with  a  mean  of  14  minutes  and  loading  time  for  server  2  exponentially  distributed  with  a 
mean  of  12  minutes.  When  a  truck  has  been  loaded,  it  hauls  the  material  to  the  dumping 
area,  dumps  its  loat,  then  returns  for  more  material.  Hauling  time  is  normally  distributed 
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with  a  mean  of  22  minutes  and  standard  deviation  of  3  minutes.  The  time  to  dump  is 
uniformly  distributed  between  2  and  8  minutes,  and  return  time  is  normally  distributed 
with  a  mean  of  18  minutes  and  standard  deviation  of  3  minutes.  The  loader  must  rest  5 
minutes  after  loading  a  truck  (Pritsker,  1986:237).  The  SLAM  network  code  and  a  SLAM 
network  diagram  are  shown  in  Appendix  C.  Figure  1  is  a  diagram  of  the  simulation. 


Truck  Return 


Loader  Return 

Figure  1.  Loader  Model  Diagram. 


The  Loader  model  was  changed  to  ease  the  process  of  animating.  In  the  original 
model  there  was  no  way  to  keep  track  of  the  trucks  and  loaders;  however,  that  information 
was  needed  for  the  animations.  Also,  the  original  model  allowed  two  trucks  to  be  loaded 
at  the  same  time  and  more  than  one  truck  to  dump  at  the  same  time.  Finally,  the  loader 
with  the  longest  idle  time  was  selected  when  a  loading  operation  could  begin. 

The  first  problem  was  solved  by  assigning  a  truck  number  (1, 2, 3,  or  4)  to  attribute 
10  of  the  truck  entities  and  assigning  a  loader  number  (5  or  6)  to  attribute  11  of  the  loader 
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entities.  The  loaders  and  trucks  were  placed  into  the  appropriate  queues  at  the  beginning 
of  the  simulation  using  ENTRY  statements.  Each  load  entity  had  the  following  attributes 
assigned: 

•  Attribute  1  -  Not  used  (originally  was  to  be  used  as  a  load  counter) 
e  Attribute  2  •  Next  load  arrival  time 

e  Attribute  3  -  Not  used 
e  Attribute  4  •  Loader  1  load  time 
e  Attribute  5  •  Loader  2  load  time 
e  Attribute  6  -  Loader  rest  time  (always  5) 

•  Attribute  7  •  Truck  hauling  time 

•  Attribute  8  -  Truck  dumping  time 

i 

•  Attribute  9  *  Thick  return  time 

Since  two  loads  were  accumulated  at  the  Accumulate  node,  only  the  attributes  of  the  second 

[ 

load  entity  arriving  were  saved.  When  two  loads,  a  truck,  and  a  loader  were  available 

! 

the  Select  node  summed  the  attributes  of  all  three  creating  a  new  entity  with  the  above 
attributes  plus  the  truck  number  and  loader  number  in  attributes  10  and  11  respectively. 
There  was  one  resource  representing  the  loading  area  so  the  new  entity  had  to  wait  until 
that  resource  (LOADSPOT)  was  available.  This  solved  the  problem  of  two  trucks  being 
loaded  at  once.  The  appropriate  loader  server  was  selected  based  on  Attribute  11.  The 
resource  LOADSPOT  was  then  freed,  and  the  entity  was  cloned  with  the  loader  being 
delayed  for  5  minutes.  All  of  the  attributes  of  the  loader,  except  for  Attribute  11  (loader 
number),  were  zeroed  out  before  the  loader  was  routed  back  to  the  loader  queue.  The 
second  entity  encountered  the  hauling  delay,  waited  for  the  resource  DUMPSPOT  (of  which 
there  was  one),  delayed  for  dumping,  freed  DUMPSPOT,  and  then  had  its  return  delay. 
Before  this  entity  was  routed  back  to  the  truck  queue,  all  of  its  attributes  except  Attribute 
10  (truck  number)  were  zeroed  out.  Therefore,  assigning  attributes  allowed  tracking  of 
•he  trucks  and  loaders,  and  the  resources  limited  loading  and  dumping  to  one  truck  at  a 
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time.  The  selection  of  loader  based  on  longest  idle  time  changed  to  alternating  loaders. 
Reasons  why  these  changes  were  needed  will  be  given  when  the  animations  are  discussed 
because  the  problems  were  not  encountered  until  the  animations  were  begun.  The  diagram 
in  Figure  1  was  still  correct  even  after  the  modifications.  These  modifications  significantly 
increased  the  length  of  the  network  code,  plus  Fortran  subroutines  were  needed  to  create 
the  trace  file  for  Proof.  The  modified  network  code  is  shown  in  Appendix  C. 

Animations.  Seven  different  animations  were  created  to  examine  movement,  color, 
and  detail  of  icons: 

•  M  -  Movement.  Simple  icons  that  move  but  do  not  change  level  of  detail  or  color. 

•  I  -  Icon.  Icons  exhibit  differing  levels  of  detail  but  do  not  move  or  change  color. 

•  C  -  Color.  Simple  icons  that  change  color  but  do  not  move  or  change  level  of  detail. 

•  MI  *  Movement  and  Icon.  Icons  move  and  change  level  of  detail  but  do  not  change 
color. 

•  MC  •  Movement  and  Color.  Simple  icons  that  move  and  change  color  but  do  not 
change  level  of  detail. 

a  Cl  -  Color  and  Icon.  Icons  change  color  and  level  of  detail  but  do  not  move. 

•  MCI  •  Movement,  Color,  and  Icon.  Icons  move  and  change  color  and  level  of  detail. 

Two  files  are  required  to  run  a  Proof  animation:  a  layout  file  (.LAY)  and  a  trace 
file  (.ATF'.  With  regard  to  the  layouts,  the  only  background  object  created  was  an  object 
representing  a  pile  of  material  at  the  dumping  area.  This  object  only  appeared  in  MI  and 
MCI.  Any  object  on  the  screen  that  could  move  or  change  status  in  some  way  required  a 
unique  number.  That  is  why  it  was  necessary  to  keep  track  of  the  trucks  and  loaders  in 
the  simulation  model.  Object  classes  were  created  that  represented  loaded  and  unloaded 
trucks  and  loaders.  When  a  particular  object  was  needed,  a  class  was  assigned  by  the 
trace  file  tc  the  object's  number.  The  object  that  appeared  on  the  screen  was  based  on  the 
class  assigned  to  the  object’s  number.  An  object’s  color,  speed,  or  travel  time  could  also 
be  changed  by  assigning  a  different  color  or  speed  to  the  object’s  number.  For  instance, 


“SET  10  COLOR  RED"  where  10  is  the  object  number,  or  “SET  10  TRAVEL  9.2"  where 
9.2  is  9.2  minutes  are  examples  of  trace  file,  commands  that  change  object  status.  Figure  2 
shows  the  different  objects  classes  used. 


□  □ 


Simple  Detailed  Track  Detailed  Track  Detailed  Track 
Track  Unloaded  Partially  Loaded  Fully  Loaded 


▲ 

Load 


□  □ 


Simple  Detailed  Loader  Detailed  Loader 
Loader  Unloaded  Loaded 


Figure  2.  Animation  Object  Classes. 


For  the  animations  with  movement,  paths  were  created  for  the  objects  to  follow. 
Allowing  two  trucks  to  load  or  dump  at  the  same  time  greatly  complicated  the  path 
structure  and  the  logic  for  path  usage.  That  is  why  the  trucks  were  limited  to  loading  and 
damping  one  at  a  time.  Also,  the  purpose  of  this  research  was  to  examine  animations,  not 
modeling  assumptions.  A  snapshot  of  an  animation  with  movement  is  in  Figure  3.  The 
icons  shown  are  the  ones  for  MI  and  MCI.  No  bulldozer  icon  was  used.  The  bulldozer  was 
represented  by  the  loads  arriving.  The  loading  times  were  divided  by  six  and  assigned  to 
a  portion  of  the  loading  process.  For  example,  if  the  loading  time  was  12  minutes,  the 
following  would  occur: 

t  • 

•  The  loader  would  leave  its  idle  position  and  travel  for  two  minutes  to  pick  up  a  load. 


•  The  loader  class  would  change  to  loaded,  and  the  loader  would  travel  for  two  minutes 
to  the  truck. 

•  The  loader  would  wait  at  the  truck  for  two  minutes  to  represent  the  transfer  of  the 
load  to  the  truck. 


..vC.~V  . 
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•  The  loader  class  would  change  to  unloaded  and  then  travel  for  two  minutes  to  pick 
up  the  second  load,  and  the  truck  would  change  class  to  show  one  load. 

•  The  loader  class  would  change  to  loaded,  and  the  loader  would  travel  for  two  minutes 
to  the  truck. 

•  The  loader  would  wait  at  the  truck  for  two  minutes  to  represent  the  transfer  of  the 
second  load  to  the  truck. 

e  The  loader  class  would  change  to  unloaded  and  then  travel  for  five  minutes  to  its  idle 
position,  and  the  truck  would  change  class  to  show  two  loads  and  begin  hauling. 

The  two  loads  “fell”  out  of  the  truck  onto  the  pile  when  dumping.  M  and  MC  did  not 
show  the  pile  of  material,  and  the  icons  were  the  ones  labeled  “Simple  Truck”  and  “Simple 
loader”  in  Figure  2;  however,  they  did  have  the  same  movement.  For  MC,  the  icons 
were  white  when  idle,  green  when  travelling  empty,  red  when  travelling  loaded,  pink  when 
partially  loaded  (trucks  only),  and  yellow  when  dumping  (trucks  only).  MCI  had  the  same 
coIot  scheme  with  the  icons  as  described  above.  The  icons  remained  red  for  M  and  MI. 


load  to  waiting  trade 


Load 

Queue 


Tract  Hauling 


Me  Loader 


Color  Key 
here  for 
MC  sod  MCI 


Track  Returning 

0 


Track  Damping 


il 


Pile  of  Material 


Figure  3.  Sample  Animation  With  Movement. 
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The  animations  without  movement  used  larger  stationary  icons.  The  four  tracks  were 
displayed  at  the  top  of  the  screen  and  the  loaders  at  the  bottom  of  the  screen.  Figure  4 
gives  a  representation  of  the  animations  without  movement  (C,  I,  and  Cl).  The  figure 
shows  the  icons  for  I  and  Cl.  As  with  M  and  MC,  the  icons  for  C  were  the  ones  labeled 
“Simple  Track”  and  “Simple  Loader”  in  Figure  2.  For  C  and  Cl,  the  color  changes  were 
the  same  as  for  MC  and  MCI,  and  the  icon  changes  for  I  and  Cl  were  the  same  as  for  MI 
and  MCI.  The  icons  in  I  remained  red.  There  was  no  representation  of  the  load  queue  in 
the  stationary  animations. 


Figure  4.  Sample  Animation  Without  Movement. 


Scenarios.  Problems  were  incorporated  into  the  experiment  to  motivate  the  subjects 
to  concentrate  on  the  animations  and  to  measure  the  performance  of  the  different  anima¬ 
tions.  The  subjects  were  the  “system  experts”  (owner,  operator,  or  manager  of  the  loading 

operation);  thus,  they  were  looking  for  efficiency  problems  with  their  system.  TVible  1  gives 

\ 

the  scenario  abbreviation,  the  simulation  model  modification,  and  the  associated  system 

1 

problems  for  each  scenario.  i 

■ 

,  l 

* 

Experimental  Design.  The  order  in  which  the  animations  were  viewed  was  random¬ 
ized  as  well  as  the  scenario  associated  with  each,  animation.  Therefore,  there  were  49 
different  animation  and  scenario  combinations.  Table  2  shows  the  experimental  design 


Table  1.  Scenario  Descriptions. 


|  Scenario 

Model  Modifications 

System  Problems 

LDEft 

Load  interarrival  time,  load¬ 
ing  time,  hauling  time,  damp¬ 
ing  time,  and  return  time 
unchanged.  (Original  Model 
with  Modifications) 

All  resources  are  used  ade¬ 
quately  with  a  slight  buildup 
of  loads. 

SLTK 

Haaling  time,  damping  time, 
and  return  time  doubled  for 
the  third  track.  (Slow  Track) 

One  truck  is  much  slower  than 
the  others,  and  there  is  a 
buildup  of  loads. 

SLLD 

Loading  time  doubled  for  the 
first  loader.  (Slow  Loader) 

One  loader  is  too  slow,  which 
creates  a  buildup  of  loads  and 
Idle  time  for  the  tracks. 

FT 

Hauling  time,  dumping  time, 
and  return  time  cut  in  half  for 
all  tracks.  (Fast  Trucks) 

There  are  too  many  trucks 
for  the  number  of  loaders  and 
bulldozers  or  not  enough  load¬ 
ers  and  bulldozers  for  the  num¬ 
ber  of  trucks. 

FL 

Loading  times  halved  for  both 
loaders.  (Fast  Loaders) 

There  are  too  many  loaders  for 
the  number  of  trucks  and  bull¬ 
dozers  or  not  enough  tracks 
and  bulldozers  for  the  number 
of  loaders. 

SL 

Load  interarrival  time  multi¬ 
plied  by  2.  (Slow  Loads) 

There  are  not  enough  bulldoz¬ 
ers,  which  creates  idle  time  for 
the  tracks  and  loaders. 

si  j 

Hauling  time,  dumping  time, 
and  return  time  doubled  for  all 
tracks.  (Slow  Trucks) 

There  are  too  few  tracks, 
which  creaks  a  buildup  of 
loads  and  loader  idle  time. 

constructed  for  50  subjects.  All  49  animation  and  scenario  combinations  had  been  selected 
by  run  32.  A  run  consisted  of  the  seven  animation  and  scenario  combinations  to  be  viewed 
by  a  subject. 

Data  Collection  Forms  and  Pictures. 

Forms.  Several  forms  were  needed  for  this  experiment,  all  of  which  are  given 
in  Appendix  A.  The  first  form  created  was  the  Scenario/ Animation  Viewing  form.  At  the 
top  were  places  for  the  subject’s  name,  the  date,  and  the  run  number.  Also,  there  was  a 
selection  for  student,  faculty,  or  other.  Following  this  there  were  seven  groups  of  lines  used 
for  listing  the  scenario,  animation,  problem  identification  time,  and  problem  observed.  At 
the  bottom  of  the  form  was  a  place  for  recording  the  subject’s  selection  of  “best’’  and 
“worst”  animation  type.  This  form  was  completed  by  the  experimenter  while  a  subject 
viewed  the  animations. 

The  second  form  was  the  pairwise  comparison  form.  This  form  consisted  of  three 
pages,  the  first  of  which  included  instructions  for  completing  the  form  and  a  listing  of  the 
animation  types.  The  second  page  gave  some  pairwise  comparison  examples,  and  the  last 
page  was  the  form  itself.  The  same  information  was  at  the  top  of  this  last  page  as  was  at 
the  top  of  the  Scenario/Animation  Viewing  form.  The  pairwise  form  contained  21  lines. 
This  allowed  for  the  comparison  of  each  animation  type  with  every  other  animation  type. 
No  animation  type  was  compared  with  itself.  Each  line  had  17  blanks  with  the  ninth  (or 
middle)  blank  marked  off  in  a  column  by  itself.  An  animation  type  was  listed  on  each  side 
of  the  line.  If  the  subject  considered  the  animation  types  to  be  equal  in  contributing  to 
face  validity  (that  is,  both  communicated  equally  well),  he  or  she  would  mark  the  middle 
column.  If  one  animation  type  was  preferred  over  the  other,  then  the  subject  would  mark 
in  a  space  on  that  animation  type’s  side.  The  more  the  subject  preferred  the  animation 
type,  the  further  out  to  the  side  he  or  she  would  mark. 

The  pairwise  comparison  form  provided  the  data  for  the  pairwise  matrix  (or  judgment 
matrix),  which  was  the  third  form.  The  pairwise  matrix  is  the  one  described  in  Chapter 
2.  The  top  of  the  form  included  the  subject’s  name,  the  date,  and  the  run  number.  The 


mation 


matrix  bad  the  rows  *nd  columns  headed  by  the  animation  types.  Ones  were  placed  on 
the  diagonal  since  it  was  assumed  that  an  animation  type  was  t'jual  with  itself.  The  lower 
portion  was  shaded  because  the  ji  element  is  simply  the  reciprocal  of  the  ij  element.  Places 
were  provided  to  the  side  of  the  matrix  for  the  calculated  ratings.  The  pairwise  matrix 
only  is  shown  in  Figure  5.  Finally,  a  description  of  the  model  was  written  for  the  subject 
to  read.  This  description  included  the  diagram  in  Figure  1,  a  discussion  of  face  validity 
and  the  subject’s  role  as  a  "system  expert”,  and  a  description  of  the  animation  types. 


Pictures.  Pictures  of  each  animation  type  were  used  at  the  beginning  to  explain 
each  animation  type,  and  the  subject  referred  to  the  pictures  when  completing  the  pairwise 
comparisons  of  the  animations.  The  pictures  were  8  x  10s  and  were  taken  at  the  same  point 
in  simulated  time  using  the  LDER  scenario. 


\ 


The  Experiment 

The  experiment  was  conducted  in  an  isolated  room,  and  a  “Do  Not  Enter”  sign 
was  placed  on  the  outside  of  the  door  to  preclude  interruptions.  Figure  6  shows  the 
physical  setup.  The  computer  was  against  the  wall  with  the  pictures  of  the  animation 
types  attached  to  the  wall  above  the  computer.  Control  of  the  computer  was  maintained 
by  the  researcher,  who  was  seated  next  to  the  subject.  This  allowed  the  researcher  to  start 
the  animations  and  the  stopwatch  at  the  same  time.  The  subject  was  verbally  given  the 
purpose  of  the  experiment  and  then  asked  to  read  the  description  of  the  model.  While 
the  subject  was  reading,  the  information  requested  at  the  top  of  each  form  was  filled  in 
by  the  researcher,  and  a  ran  number  was  assigned.  Once  the  description  was  read  and 
questions  were  answered  (if  any),  the  animations  were  described  using  the  pictures.  Each 
animation  was  viewed  for  one  minute,  which  equated  to  350  simulated  minutes.  The 
stopwatch  was  started  when  each  animation' began.  This  allowed  the  time  to  be  recorded 
if  the  subject  identified  a  problem.  If  no  problem  was  observed,  60  seconds  was  entered  for 
the  problem  identification  time  and  “no  problem  observed”  was  entered  for  the  problem 
observed  portion  of  the  form.  After  viewing  all  the  animations,  the  subject  performed  the 
pairwise  comparisons.  The  subjective  selection  of  the  animation  types  that  communicated 
the  “best”  and  the  “worst”  concluded  the  experiment.  The  total  time  required  was  25  to 
30  minutes.  The  pairwise  matrix  was  filled  in,  and  the  ratings  were  calculated  after  the 
subject  left. 

Analysis  Methods 

j 

The  analysis  was  divided  into  three  sections  based  on  the  type  of  data.  Summary 
statistics  such  as  numbers  and  percentages  were  calculated  for  the  subjective  evaluation 
data  (selection  of  “best”  and  “worst”  animations).  The  AHP  ratings  were  calculated  using 
the  SWORD  computer  program  provided  by  Vidulich  (SWORD,  1989).  This  program  also 
calculated  the  52  value  and  the  critical  5*  value.  Simple  statistics  were  calculated  for 
the  AHP  ratings;  plus,  Analysis  of  Variance,  Factor  Analysis,  and  Cluster  Analysis  was 
performed.  The  third  set  of  data  was  the  problem  data.  Simple  statistics  and  summary 
statistics  were  calculated  for  this  data  also.  In  addition,  the  data  was  analyzed  using  Anal- 
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Figure  6.  Room  Setup. 

you  of  Variance  and  Cluster  Analysis  from  the  animation  perspective.  Analysis  of  Variance 
wii  used  for  analyzing  the  problem  data  from  the  scenario  perspective.  The  animation 
perspective  looked  at  the  performance  of  the  subjects  when  viewing  the  animations  with 
regard  to  time  of  problem  identification  and  accuracy  of  problem  identification  by  isolating 
the  animation  type.  The  scenario  perspective  examined  the  animation  types  by  isolating 
the  scenario. 


IV.  Results  and  Discussion 

A  total  of  47  individuals  volunteered  to  ’/lew  the  animations.  Of  the  47,  41  were 
AF1T  students,  five  were  AFIT  faculty,  and  one  was  neither.  The  results  from  the  subjec¬ 
tive  evaluation,  including  the  “best”'  and  “worst”  selections  and  the  pairwise  comparison 
(AHP  ratings),  will  be  presented  first.  The  results  of  the  objective  evaluation  (problem 
identification  accuracy  and  times)  will  follow,  and  a  comparison  of  the  subjective  and 
objective  results  will  conclude  the  chapter. 

Subjective  Results 

" Best 9  and  "Worst9.  Only  three  animation  types  were  chosen  as  communicating  the 

"i 

best:  MCI,  MI,  and  M.  Two  animation  types  were  identified  as  worst:  C  and  I.  Table  3 
gives  the  results  for  each  animation  type. 

Table  3.  Number  of  “best”  and  “worst”  animation  types. 


Rating 

C 

I 

Cl  I  MC 

M 

MI 

MCI 

Best 

0 

o 

0 

0 

1 

9 

37 

Worst 

36 

0 

0 

0 

0 

0 

The  most  preferred  animation  type  was  MCI,  which  was  selected  as  the  best  37  times 
(79%).  Chosen  36  times  (77%)  as  the  animation  type  preferred  the  least  was  C. 

AHP  Ratings.  All  of  the  judgment  matrices’  S3  values  except  one  passed  the  con¬ 
sistency  criterion.  However,  the  judgment  matrix  that  did  not  pass  was  used.  This  is 
because  the  S3  valuo  was  only  slightly  above  the  critical  value.  Discussions  with  Vidulich 
revealed  that  slightly  inconsistent  matrices  did  not  affect  the  analysis  enough  to  warrant 
not  using  the  data  (Vidulich,  1992).  Table  4  gives  the  mean  and  standard  deviation  of  the 
ratings  for  each  animation  type  while  Figure  7  is  a  scatter  plot  of  the  ratings.  There  is 
a  dour  distinction  between  animation  types  without  movement  and  animation  types  with 
movement.  The  subjects  demonstrated  the  highest  preference  for  MCI.  However,  most 
Stated  that  movement  was  the  key  to  understanding  the  model 
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Table  4.  Animation  Rating  Means  and  Standard  Deviations. 

Statistic  |  C  1  I  Cl  M  MC  MI  MCI 

Mean  0.0287  0.0353  0.0561  0.1221  0.1693  0.2340  0.3544~ 

St  Dev  0.0195  0.0130  0.0243  0.0576  0.0478  0.0693  0.0658 


Avg  Rating 
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Animation  Type 


Figure  7.  Scatter  Plot  of  Animation  Ratings. 


AUF  AWOT'A.  This  Analysis  of  Variance  assumed  a  model  of  the  form 

y<>  =  /i+r<  +  flj+£,>  (2) 

where: 

•  ft  is  the  mean  of  the  responses 

•  Ti  are  constants  for  the  treatment  effects 

•  Bj  are  constants  for  the  block  effects 

•  By  are  independent  N( 0,  a3) 

•  i  *  1,...,7;  j  =  1,...,47 

Thus,  the  ratings  are  assumed  to  be  independent  and  normally  distributed,  with  mean 
and  constant  variance 

W;}  *  H  +  Ti  +  B, 

<r2{Yij}  =  <ra 

Each  animation  type  was  a  treatment  (s’  =  1,...,7),  and  each  subject  was  a  block  (j  = 
1,...,47).  So  this  model  determined  if  there  were  differences  among  the  animation  types 
and  differences  between  subjects.  Table  5  is  the  ANOVA  table. 


Thble  5.  Analysis  of  Variance  Table  for  Animation  Ratings. 


Source 

DF 

SS 

MS 

F 

P 

Treatment 

Block 

Error 

6 

46 

276 

4.05716 

3.429E*06 

0.73014 

0.67619 

7.453E-08 

0.00264 

255.61 

0.00 

0.0000 

1.0000 

Tbtal 

Grand  Average 

328 

1 

4.78731 

6.71428 

The  ANOVA  table  shows  that  there  is  a  difference  between  animation  types.  The 
interesting  result,  though,  i*  that  there  is  no  block  effect.  That  is,  there  was  no  significant 
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difference  between  subjects.  A  Tuiey  test  for  additivity  was  performed  to  test  for  interac¬ 
tion  effects  between  animation  type  and  subject.  Since  the  null  hypothesis  of  the  test  is  no 
interaction  effects  (treatment  and  block  are  additive),  the  F  value  of  0.04  and  the  P  value 
of  0.8370  indicated  there  was  no  animation  type  and  subject  interaction  effects.  Table  6 
gives  the  complete  results  of  the  test. 

Table  6.  Tukey’s  1  Degree  of  Freedom  Test  For  Additivity. 


Source 

DF 

SS 

F 

P 

Nonadditivity 

1 

1.126E-04 

0.04 

0.8370 

Remainder 

275 

0.73003 

Figure  8  shows  a  normal  probability  plot.  The  residuals  were  plotted  against  the 
expected  residuals  under  normality  to  check  the  normality  assumption.  The  almost  straight 
line  indicates  that  the  normality  assumption  is  valid.  Another  requirement  for  ANOVA 
is  equal  variances  for  the  animation  types  (treatments).  The  residuals  are  plotted  against 
the  animation  types  in  Figure  9.  The  plot  seems  to  show  different  variances.  However, 
there  is  a  good  scattering  of  the  residuals  in  Figure  10,  which  shows  the  residuals  versus 
the  subjects  (blocks). 

Since  the  plot  of  the  residuals  versus  treatment  for  (2)  seemed  to  indicate  different 
variances,  a  Bartlett’s  test  of  equal  variances  was  performed  after  dropping  the  block 
variable.  The  test  indicated  that  the  constant  variance  assumption  for  the  ANOVAs  was 
not  valid.  The  test  yielded  a  high  x2  value  (180.31)  and  a  P  value  of  0.0,  so  the  null 
hypothesis  of  equal  variances  for  each  animation  type  was  rejected.  The  results  of  the  test 
are  shown  in  I^tble  7.  However,  because  the  sample  sizes  are  the  same  for  each  animation 
type,  the  model  still  can  be  used.  Neter  and  others  assert  that  the  effects  of  unequal 
variances  are  minimized  when  the  sample  sizes  are  the  same.  They  state,  “When  the  error 
variances  are  unequal,  the  F  test  for  the  equality  of  means  with  the  fixed  ANOVA  model 
is  only  slightly  affected  if  all  factor  level  sample  sizes  are  equal  or  do  not  differ  greatly” 
(Neter  and  others,  1990:624). 

The  Tukey  method  of  multiple  comparisons  was  used  to  determine  if  the  differences 
between  the  mean  ratings  of  each  animation  type  were  statistically  significant.  The  means 
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Figure  8.  Normal  Probability  Plot  of  Residuals  vs  Expected  Residuals. 


Animation  I^pe 
(Treatment) 

figure  9.  Scatter  Plot  of  Residuals  vs  Animation  Type  (Treatment). 


were  tested  at  a  family  a  level  of  0.1.  The  differences  between  the  mean  ratings  of  the 
animations  without  movement  (C,  I,  and  Cl)  were  determined  not  to  be  statistically  signif¬ 
icant.  That  is,  at  the  90%  confidence  level  the  means  were  considered  to  be  the  same.  Each 
of  the  mean  ratings  for  the  animation  types  with  movement  (M,  MC,  MI,  and  MCI)  were 
statistically  different  from  each  other  and  the  group  of  animation  types  without  movement. 
Thus,  the  test  determined  that  there  were  five  groups:  MCI,  MI,  MC,  M,  and  (C,  I,  Cl). 

Since  there  was  no  block  effect,  and  there  was  no  treatment  and  block  interaction 
effects,  block  was  dropped.  So  the  model  became 

Yij  =  n  +  Ti  +  a;  (3) 

where: 

•  n  is  the  mean  of  the  responses 

•  T,  are  constants  for  the  treatment  effects 

•  Cij  are  independent  N( 0,  a3) 

a  ™  l|*"j  47 

Again,  the  ratings  are  assumed  to  be  independent  and  normally  distributed,  with  mean 
and  constant  variance 

E{Yh)  »  /i  +  Ti 
•*{*«}  =  ff3 

Table  8  is  the  ANOVA  table  for  (3).  This  ANOVA  also  indicated  that  there  is  an 
effect  based  on  animation  type.  That  is,  there  is  a  difference  in  preferences  for  animation 
types.  So  the  ANOVA  and  associated  tests  show  that  animations  without  movement  are 
preferred  the  least.  Since  these  animations  did  not  communicate  as  well  as  the  animations 
with  movement,  they  were  not  considered  to  be  high  contributors  to  understanding  the 
model.  Understanding  the  model  operation  aids  in  assessing  face  validity.  Preference  for 
animations  with  movement  varied  based  on  what  other  aspects  of  animation  were  included 
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(color  and  icon).  Movement  was  judged  as  a  necessary  contributor  to  face  validity,  but 
adding  other  aspects  were  judged  to  improve  the  performance  of  movement  in  communi¬ 
cating  .ne  operation  of  the  model. 

Table  8.  One-Way  Analysis  of  Variance  Table  for  Animation  Ratings. 


Source 

DF 

SS 

MS 

F 

P 

Between 

Within 

6 

322 

4.05716 

0.73014 

0.67619 

0.00226 

298.21 

0.0000 

Total 

328 

4.78731 

AHP  Factor  Analysis.  The  sample  covariance  matrix  was  used  for  the  fac¬ 
tor  analysis  since  all  the  ratings  were  in  the  same  units  (preference).  Three  factors  were 
retained  because  93%  o<  the  variance  was  explained  by  the  first  three  factors.  The  eigen¬ 
values,  proportion  of  variance  explained,  and  cumulative  proportion  of  variance  explained 
are  given  in  Table  9.  There  was  a  significant  drop  from  the  third  to  the  fourth  eigenvalue, 
and  only  a  small  increase  in  proportion  of  variance  explained.  The  three  factors  that  were 
kept  were  varimax  rotated  and  are  shown  in  Table  10. 


Table  9.  Eigenvalues  of  the  Animation  Ratings  Covariance  Matrix. 


Eigenvalue 

0.0081 

0.0046 

0.0021 

0.0008 

0.0002 

0.0001 

0.0000 

Proportion 

0.5088 

0.2872 

0.1325 

0.0501 

0.0131 

0.0082 

0.0000 

Cumulative 

0.5088 

0.7960 

0.9285 

0.9787 

0.9918 

1.0000 

1.0000 

Table  10.  Rotated  Animation  Rating  Factors. 


Animation 

Factor  1 

Factor  2 

Factor  3 

M 

-0.07000 

0.96439 

-0.18944 

C 

0.61113 

-0.00739 

-0.04354 

I 

-0.03858 

0.02393 

-0.27963 

>0.82203 

-0.06075 

-0.56269 

0.87741 

-0.36355 

-0.13631 

0.48971 

0.08657 

0.27445 

I  MCI 

-0.06252 

-0.54858 

0.82630 

The  animation  types  which  load  the  highest  on  Factor  1  are  C,  MI,  MC,  and  CL  The 
sign  of  MI  is  opposite  to  that  of  the  others,  though.  The  animation  types  C,  Cl,  and  MC 
include  color  and  are  positive  in  sign.  Although  MCI  contained  color  and  was  negative  in 
sign  like  MI,  it  did  not  load  high  on  Factor  1.  Movement  loads  the  highest  on  Factor  2,  and 
MCI  and  I  load  the  highest  on  Factor  3;  however,  *0.28  for  I  is  too  small  to  be  considered 
significant.  In  addition,  MCI  loads  significantly  on  Factor  2,  but  in  opposite  sign  to  M, 
and  MI  loads  significantly  on  Factor  3,  but  in  opposite  sign  to  MCI.  So  Factor  1  appears  to 
be  measuring  preference  for  color,  in  particular,  preference  for  movement  with  color  versus 
preference  for  movement  with  icon.  Either  a  subject  liked  color  or  not,  and  this  preference 
was  most  evident  when  color  was  combined  with  movement.  Factor  2  measures  preference 
for  movement  alone.  Preference  for  MCI  dropped  the  most  as  preference  for  M  alone 
increased.  The  third  factor  shows  preference  for  MCI.  As  the  rating  for  MCI  increased, 
the  rating  for  MI  decreased.  The  plots  of  the  factor  scores  confirm  these  conclusions. 

I 

Figure  11  shows  the  plot  of  Factor  1  scores  versus  Factor  2  scores.  The  outliers  for  Factor 

1  contrast  MC  and  C  with  MI.  For  those  on  the  left,  MC  and  C  were  rated  low,  and  MI 

I 

was  rated  high.  For  the  outlier  on  the  left,  MC  and  C  were  rated  high,  while  MI  was 
rated  low.  A  look  at  the  subjective  evaluations  also  confirms  this.  The  two  on  the  left 
that  rated  MC  and  C  very  low  chose  MI  as  the  best  animation  instead  of  MCI.  M  was 
rated  extremely  high  for  the  three  outliers  at  the  top  for  Factor  2,  and  M  was  rated  low 
for  the  outlying  point  at  the  bottom-  The  plot  of  Fhctor  1  versus  Factor  3  scores  is  shown 
in  Figure  12.  The  same  outliers  are  seen  for  Factor  1.  Factor  3’s  outliers  at  the  top  rated 
MCI  very  high,  while  the  point  at  the  bottom  rated  MCI  quite  low.  figure  13  contains 
the  plot  of  Factor  2  scores  versus  Factor  3  scores.  The  same  outliers  are  seen  here,  and, 
again,  the  scores  in  which  M  was  rated  high  are  significantly  different  from  the  group  of 
scores.  Therefore,  the  pairwise  comparison  measured  preference  for  color  with  movement 
versus  icon  with  movement,  preference  for  movement,  and  preference  for  movement,  color, 
and  icon. 

AHP  Ouster  Analysis.  Three  different  cluster  analysis  techniques  were  used: 
Complete  linkage,  Average  Linkage,  and  Ward’s  Minimum  Variance.  AH  three  techniques 
clustered  the  animations  without  movement  (C,  I,  and  Cl)  right  away.  There  were  differ- 
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Figure  13.  Plot  of  Factor  2  Scores  vs  Factor  3  Scores. 

ences  in  how  M,  MC,  MI,  and  MCI  were  clustered.  The  complete  linkage  technique  and 
Ward’s  minimum  variance  technique  clustered  M  and  MC,  then  MI  with  M  and  MC.  Next, 
M,  MC,  and  MI  were  clustered  with  C,  I,  and  Cl.  Finally,  at  a  much  farther  distance,  MCI 
was  clustered  with  all  the  other  animation  types.  Average  linkage  was  similar  except  that 
MI  was  clustered  with  MCI  instead  of  M  and  MC.  MI  and  MCI  were  then  clustered  with  the 
remaining  animation  types  at  approximately  double  the  distance  that  the  other  animation 
types  were  clustered.  Table  11  shows  the  dusters  and  the  distance  for  the  complete  linkage 
technique,  and  Figure  14  is  the  dendogram  for  the  complete  linkage.  A  dendogram  is  a 
graph  that  shows  the  groupings  and  the  level  at  which  they  were  grouped.  The  distance  is 
the  normalized  maximum  distance  between  two  dusters.  Table  12  and  Figure  15  show  the 
same  for  the  average  linkage  technique,  and  Table  13  and  Figure  16  gives  the  results  for 
Ward’s  minimum  variance  technique.  The  distance  measure  for  the  average  linkage  tech¬ 
nique  is  the  normalized  root-mean-square  (RMS)  distance,  which  is  “the  root-mean-square 
distance  between  pairs  of  objects  in  the  two  dusters  joined  with  one  object  from  each 
duster”  (SAS  Institute  Inc.,  1989:567).  Far  Ward’s  technique,  the  distance  measure  is  the 
Semipart ial  R-Squared.  This  is  the  decrease  in  the  proportion  of  variance  accounted  for 


•iter  the  two  groups  have  been  joined  and  is  eqnal  to  “the  between- cluster  sum  of  squares 
divided  by  the  corrected  total  sum  of  squares”  (SAS  Institute  Inc.,  1989:567). 

The  cluster  analysis  yielded  results  similar  to  Tukey’s  method  of  multiple  compar¬ 
isons.  There  was  very  little  difference  in  preference  for  animations  without  movement. 
That  is,  C,  I,  and  Cl  were  not  judged  to  be  high  contributors  to  face  validity.  With  re¬ 
gards  to  movement,  the  differences  in  preferences  for  M,  MC,  MI,  and  MCI  were  large 
enough  for  clustering  to  occur  after  C,  I,  and  Cl  were  clustered.  Also,  for  complete  linkage 
and  Ward’s  minimum  variance,  M,  MC,  and  MI  were  grouped  together  before  they  were 
grouped  with  C,  I,  and  Cl.  Thus,  movement  was  judged  as  significant.  As  with  the  factor 
analysis,  average  linkage  shows  that  icon  with  movement  was  considered  to  communicate 
better  than  color  with  movement.  This  is  seen  in  that  MI  was  clustered  with  MCI  instead 
of  M  and  MC.  The  primary  benefit  of  the  cluster  analysis  is  that  it  confirmed  that  all  of 
the  animation  types  without  movement  were  judged  to  be  basically  the  same.  C,  I  and  Cl 
were  considered  to  be  poor  animations  to  use  in  communicating  the  operation  of  a  model. 

Ifcble  11.  Complete  Linkage  Cluster  Analysis  Results. 


Table  12.  Average  Linkage  Cluster  Analysis  Results. 


.  /  Animations  Joined 

Distance 

C,I 

0.135478 

C,  L  Cl 

0.184487 

M,  MC 

0.529770 

C,  I,  Cl,  M,  MC 

0.662877 

MI,  MCI 

0.888342 

C,  I,  Cl,  M,  MC,  MI,  MCI 

1.311767 

Semi-Putial  R-Squarcd 


Table  13.  Ward’s  Minimum  Variance  Cluster  Analysis  Results, 


C  I  a  M  MC  MI  Ma 
Animation  Type 

Figure  16.  Ward’s  Minimum  Variance  Dendogram. 
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Objective  Results 

Problem  Results  From  Animation  Perspective. 

Summary  Statistics.  Given  an  animation  type,  two  items  were  looked  at  when 
the  problem  identification  data  was  analysed:  problem  identification  time  and  problem 
identification  accuracy.  The  AHP  ratings  measured  the  subjective  preference  for  an  anima¬ 
tion  type.  The  problem  data  measured  the  subjects'  objective  performance  when  viewing 
the  animations.  Figure  17  shows  the  average  time  (in  seconds)  that  a  potential  problem 
was  identified,  and  Table  14  gives  the  corresponding  numbers.  The  animations  with  move¬ 
ment  were  consistently  lower  than  those  without  movement.  The  percentage  of  problems 
correctly  identified  is  shown  in  Figure  18.  Again,  a  clear  difference  can  be  seen  between 
animations  without  movement  and  animations  with  movement  Problems  were  identified 
correctly  just  better  than  50%  of  the  time  when  the  animation  type  was  C,  I,  or  Cl,  whereas 
the  percent  correct  was  80%  or  better  for  M,  MC,  MI,  and  MCI. 

Table  14.  Problem  Identification  Time  Means  and  Standard  Deviations. 


Statistic 

C 

I 

Cl 

M 

MC 

MI 

MCI 

Mean 

47 

48 

46 

36 

30 

~33~ 

34 

St  Dev 

15 

14 

18 

L 21 

13 

15 

16 

Subject  performance  based  on  animation  type  was  also  examined  by  analyzing  the 
data  when  the  problems  were  correctly  identified  and  incorrectly  identified.  Given  that 
the  problem  was  correctly  identified,  the  animations  with  movement  (M,  MC,  MI,  and 
MCI)  overall  had  a  higher  number  of  correct  identifications  and  a  lower  average  time  in 
which  those  problems  were  recognized  while  the  animations  without  movement  (C,  I,  and 
Cl)  showed  just  the  opposite.  Figure  19  shows  the  results  graphically.  The  graph  for  the 
incorrect  identifications  is  shown  in  Figure  20.  Again,  with  regards  to  number  incorrect, 
the  subjects  performed  better  when  the  animations  had  movement.  However,  the  average 
time  of  an  incorrect  identification  was  similar  for  all  the  animation  types  with  M  the 
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Figure  17.  Average  Time  to  Identify  Potential  Problem. 
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Figure  18.  Percentage  of  Problems  Correctly  Identified. 
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Figure  19.  Number  and  A”erage  Time  of  Problems  Correctly  Identified. 
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Figure  20.  Number  and  Average  Time  of  Problems  Incorrectly  Identified. 
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Another  question  of  interest  was  how  many  subjects  responded  before  and  after  the 
animations  ended.  Since  the  animations  ran  for  60  seconds,  the  issue  was  how  many 
subjects  responded  within  60  seconds.  Figure  21  shows  the  results.  C,  I,  and  Cl  were 
about  even.  About  half  responded  within  60  seconds  and  about  half  after  60  seconds.  M, 
MC,  MI,  and  MCI  were  quite  different.  There  were  significantly  more  responses  within  60 
seconds  than  after  60  seconds. 

B  Ans  Before  60  Sec 
§3  Ans  After  60  Sec 


Animation 

Figure  21.  Number  of  Problem  Identifications  Before  and  After  60  Seconds. 

Therefore,  when  the  animations  contained  movement,  the  subjects  performed  better 
than  when  the  animations  did  not  contain  movement.  Unlike  the  AHP  ratings,  though, 
there  did  not  seem  to  be  any  significant  difference  between  M,  MC,  MI,  and  MCI. 

Animation  Perspective  ANOVA.  This  Analysis  of  Variance  assumed  the  same 
model  as  the  AHP  ANOVA,  which  is  given  in  Equation  (2).  In  this  case  the  times  Yij  are 
assumed  to  be  independent  and  normally  distributed,  with  mean  and  constant  variance 

—  P  +  Ti  +  Bj 
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As  before,  each  animation  type  was  a  treatment  (t  =  1. . . . ,  7),  and  each  subject  was  a  block 
(i  =  1, . .  .,47).  So  this  model  determined  if  there  were  differences  in  subject  performance 
given  an  animation  type  and  differences  in  performance  between  subjects.  Table  15  shows 
the  ANOVA  table. 

Tkble  15.  ANOVA  Table  for  Problem  Identification  Times  from  Animation  Perspective. 


Source 

DF 

SS 

MS 

F 

P 

Treatment 

6 

15437.2 

2572.87 

14.92 

0.0000 

Block 

46 

26524.0 

576.61 

3.34 

0.0000 

Error 

276 

47586.4 

172.45 

Total 

328 

89557.7 

Grand  Averagt 

1 

504500 

The  ANOVA  table  shows  that  there  is  a  difference  in  subject  performance  between 
animation  types,  and  there  is  a  difference  in  performance  between  subjects.  That  is,  there 
is  a  block  effect.  This  was  expected  since  there  was  a  variety  of  subjects,  each  with  a 
different  idea  of  what  constitutes  an  efficiency  problem.  Also,  some  were  familiar  with 
simulation  modeling  and  animation  and  some  were  not.  A  Tukey  test  for  additivity  was 
performed  to  test  for  interaction  effects  between  animation  type  and  subject.  Since  the 
null  hypothesis  of  the  test  is  no  interaction  effects  (treatment  and  block  are  additive),  the 
F  value  of  0.05  and  the  P  value  of  0.8308  indicated  there  was  no  animation  type  and 
subject  interaction  effects.  Table  16  gives  the  complete  results  of  the  test. 

Table  16.  Tukey’s  1  Degree  of  freedom  Test  For  Additivity. 


Source 

DF 

SS 

F 

P 

Nonadditivity 

1 

7.91337 

0.05 

0.8308 

Remainder 

275 

47588.5 

Figure  22  shows  a  normal  probability  plot.  The  residuals  were  plotted  against  the 
expected  residuals  under  normality  to  check  the  normality  assumption.  The  almost  straight 
line  indicates  that  the  normality  assumption  is  valid.  Another  requirement  for  ANOVA 
is  equal  variances  for  the  animation  types  (treatments).  The  residuals  are  plotted  against 
the  animation  types  in  Figure  23.  The  plot  seems  to  show  equal  variances,  and  there  is  a 


44 


good  scattering  of  the  residuals  in  Figure  24,  which  shows  the  residuals  versus  the  subjects 
(blocks).  Since  the  block  variable  could  not  be  dropped  because  of  the  block  effect,  a 
Bartlett’s  test  of  equal  variances  could  not  be  performed. 


Expected  Residuals 


Figure  22.  Normal  Probability  Plot  of  Residuals  vs  Expected  Residuals. 


The  Tukey  method  of  multiple  comparisons  was  used  to  determine  if  the  differences 
between  the  mean  problem  identification  times  of  each  animation  type  were  statistically 
significant.  The  means  were  tested  at  a  family  a  level  of  0.1.  The  differences  between 
the  mean  problem  identification  times  of  the  animations  without  movement  (C,  I,  and  Cl) 
were  determined  not  to  be  statistically  significant.  That  is,  at  the  90%  confidence  level  the 
means  were  considered  to  be  the  same.  Also,  the  differences  between  the  mean  problem 


identification  times  for  the  animation  types  with  movement  (M,  MC,  MI,  and  MCI)  were 
not  statistically  different  from  each  other;  however  they  were  statistically  different  from 

the  group  of  animation  types  without  movement.  Thus,  the  test  determined  that  there 

1 

were  two  groups:  (MCI,  MI,  MC,  M)  and  (C,  I,  Cl). 


Animation  Perspective  Cluster  Analysis.  The  duster  analysis  techniques  that 


were  used  with  the  AHP  ratings  were  also  used  here:  Complete  Linkage,  Average  Linkage, 


and  Ward’s  Minimum  Variance.  AH  three  techniques  clustered  the  animations  with  move¬ 
ment  first.  M  and  MCI  were  clustered,  as  were  MI  and  MC.  Then  all  four  were  clustered. 
C  and  Cl  were  clustered,  then  joined  with  I.  This  pattern  was  the  same  for  all  the  tech¬ 
niques.  Table  17  shows  the  clusters  and  the  distance  for  the  complete  linkage  technique, 
and  Figure  25  is  the  dendogram  for  the  complete  linkage.  Table  18  and  Figure  26  show 
the  same  for  the  average  linkage  technique,  and  Table  19  and  Figure  27  gives  the  results 
for  Ward’s  minimum  variance  technique. 

The  duster  analysis  yielded  results  similar  to  Tukey’s  method  of  multiple  compar¬ 
isons.  There  was  very  little  difference  in  subject  performance  for  animations  without 
movement,  and  very  little  difference  in  subject  performance  for  animations  with  movement. 
That  is,  movement  significantly  increased  the  subjects  ability  to  discern  the  operation  of  the 
system.  Also,  adding  color,  detailed  icons,  or  both  did  not  increase  the  subject’s  ability  to 
understand  the  operation  of  the  model,  as  measured  by  the  average  problem  identification 
time.  So  the  cluster  analysis  confirmed  that  all  of  the  animation  types  without  movement 
were  more  difficult  to  interpret  than  the  animations  with  movement,  and  movement  was 
the  most  important  communication  tool. 

'foble  17.  Complete  Linkage  Cluster  Analysis  Results  for  Problem  ID  Times  from  Ani¬ 
mation  Perspective. 


Animations  Joined 

Distance 

M,  MCI 

C.749877 

MI,  MC 

0.801696 

C,  Cl 

0.887343 

M,  MCI,  MI,  MC 

0.946333 

c,  a,  i 

0.974192 

M,  MCI,  MI,  MC,  C,  I,  Cl 

1.224518 

Problems  Results  Prom  Scenario  Perspective . 

Summary  Statistics.  As  with  the  problem  results  from  the  animation  perspec¬ 
tive,  two  data  items  were  examined:  problem  identification  time  and  problem  identification 
accuracy.  Whereas  the  problem  data  looked  at  from  the  animation  perspective  measured 
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Figure  25.  Complete  Linkage  Dendogram  for  Problem  ID  Times  from  Animation 
Perspective. 


Table  18.  Average  Linkage  Cluster  Analysis  Results  for  Problem  ID  Times  from  Anima¬ 
tion  Perspective. 


Animations  Joined 

Distance 

M,  MCI 

MI,  MC 

M,  MCI,  MI,  MC 

C,CI 

C,CI,I 

M,  MCI,  MI,  MC,  C,I,CI 

0.743121 

0.794472 

0.870167 

0.879348 

0.957954 

1.087015 

Tbble  19.  Ward’s  Minimum  Variance  Cluster  Analysis  Results  for  Problem  ID  Times 
from  Animation  Perspective. 


Animations  Joined 

Distance 

M,  MCI 

0.092038 

Ml,  MC 

0.105198 

C  ,CI 

0.128875 

M,  MCI,  MI,  MC 

0.153779 

c,cu 

0.160969 

M,  MCI,  MI,  MC,  C,  I,  a 

0.359141 

.48 
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the  performance  of  the  subjects  given  an  animation  type,  the  problem  data  analyzed  from 
the  scenario  perspective  shows  which  scenarios  were  most  recognizable  to  the  subjects. 
Figure  28  shows  the  average  time  (in  seconds)  that  a  potential  problem  was  identified,  and 
Table  20  gives  the  numbers.  Three  of  the  scenarios  (FL,  LDER,  and  SLTK)  appear  to 
have  higher  times  than  the  others.  However,  only  FL  stands  out  when  the  percentage  of 
problems  correctly  identified  is  examined.  These  are  shown  in  Figure  29. 

Table  20.  Problem  Identification  Time  Means  and  Standard  Deviations  from  Scenario 
Perspective. 


Statistic 

FL 

LDER 

SLTK 

ST 

SL 

FT 

SLLD 

Mean 

“46~ 

42 

41 

35 

IT 

36 

36 

St  Dev 

16 

17 

14 

20 

17 

16 

12 

a 

a 


AvgTime 


Figure  28.  Average  Time  to  Identify  Potential  Problem  from  Scenario  Perspective. 


arto  Perspective  ANOVA.  This  Analysis  of  Variance  assumed  the  same 
model  as  the  AHP  ANOVA  and  animation  perspective  ANOVA,  which  is  given  in  Equation 
(2).  In  this  lease  the  times  Yij  are  assumed  to  be  independent  and  normally  distributed, 
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Figure  29.  Percentage  of  Problems  Correctly  Identified  from  Scenario  Perspective, 
with  mean  and  constant  variance 


E{Yn)  =  p  +  Ti  +  Bj 

*  o'1 

In  this  case,  however,  each  scenario  was  a  treatment  (i  =  1, . . . ,  7),  and  each  subject  was 
a  block  (j  =  1,...,47).  So  this  model  determined  if  there  were  differences  in  subject 
performance  among  the  scenarios  and  differences  in  performance  between  subjects.  The 
purpose  of  this  ANOVA  was  to  determine  if  there  was  a  statistical  difference  in  subject 
performance  between  the  scenarios.  Table  21  is  the  ANOVA  table. 

Table  21.  ANOVA  Table  for  Problem  Identification  Times  from  Scenario  Perspective. 


Source 

DF 

SS 

MS 

F 

P 

Treatment 

6 

5096.44 

849.407 

4.00 

0.0008 

Block 

46 

26243.0 

570.500 

2.68 

0.0000 

Error 

276 

58650.9 

212.503 

Total 

328 

89990.4 

Grand  Average 

1 

502700 

51 


The  ANOVA  table  shows  that  there  is  a  difference  in  subject  performance  between 
scenarios,  and  there  is  a  difference  in  performance  between  subjects.  That  is,  there  is 
a  block  effect.  A  Tukey  test  for  additivity  was  performed  to  test  for  interaction  effects 
between  scenario  and  subject.  Since  the  null  hypothesis  of  the  test  is  no  interaction  effects 
(treatment  and  block  are  additive),  the  F  value  of  0.14  and  the  P  value  of  0.7114  indicated 
there  was  no  scenario  and  subject  interaction  effects.  Table  22  gives  the  complete  results 
of  the  test. 

Ifcble  22.  Tukey’s  1  Degree  of  Freedom  Test  For  Additivity. 


Source 

DF 

SS 

F 

P 

Nonadditivity 

1 

29.246 

0.14 

0.7114 

Remainder 

275 

58621.7 

Figure  30  shows  a  normal  probability  plot.  The  residuals  were  plotted  against  the 
expected  residuals  uuder  normality  to  check  the  normality  assumption.  The  almost  straight 
line  indicates  that  the  normality  assumption  is  valid.  Another  requirement  for  ANOVA 
is  equal  variances  for  the  scenarios  (treatments).  The  residuals  are  plotted  against  the 
scenarios  in  Figure  31.  The  plot  seems  to  show  equal  variances,  and  there  is  a  good 
scattering  of  the  residuals  in  Figure  32,  which  shows  the  residuals  versus  the  subjects 
(blocks).  As  with  the  ANOVA  for  *'  .*  problem  identification  times  from  the  animation 
perspective,  a  Bartlett's  test  of  equal  variances  could  not  be  performed  because  of  the 
block  effect. 

The  Tukey  method  of  multiple  comparisons  was  used  to  determine  if  the  differences 
between  the  mean  problem  identification  times  of  each  scenario  were  statistically  signifi¬ 
cant.  The  means  were  tested  at  a  family  a  level  of  0.1.  The  differences  between  the  mean 
problem  identification  times  for  FL,  LDER,  and  SLTK  were  determined  not  to  be  statis¬ 
tically  significant.  That  is,  at  the  90%  confidence  level  the  means  were  considered  to  be 
the  same.  Also,  the  differences  between  the  mean  problem  identification  times  for  SL,  FT, 
SLLD,  and  ST  were  not  statistically  different  from  each  other.  Thus,  the  test  determined 
that  there  were  two  groups:  (FL,  LDER,  SLTK)  and  (SL,  FT,  SLLD,  ST).  This  grouping 
had  been  suggested  by  Figure  28. 


Expected  Residuals 

Figure  30.  Normal  Probability  Plot  of  Residuals  vs  Expected  Residuals. 
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Figure  31.  Scatter  Plot  of  Residuals  vs  Scenario  (Treatment). 
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Figure  32.  Scatter  Plot  of  Residuals  vs  Subject  (Block). 

The  higher  average  problem  identification  times  for  FL,  LDER,  and  SLTK,  can  be 
explained.  The  times  for  FL  were  higher  on  the  average  because  this  scenario  was  an  under 
utilization  problem.  In  FL  the  loaders  were  faster,  which  caused  more  loader  idle  time. 
This  caused  a  higher  average  problem  identification  time  (and  lower  problem  identification 
accuracy)  because  fewer  subjects  recognized  idle  loaders  as  a  problem.  Many  thought 
there  was  no  problem  with  FL.  In  FT  the  trucks  were  faster,  which  caused  more  truck 
idle  time;  however,  most  recognized  this  under  utilization.  This  is  probably  because  the 
trucks  waited  in  line  for  loads.  So  two  or  more  trucks  waiting  to  load  (in  the  middle  of 
the  computer  screen)  were  more  obvious  than  two  loaders  (at  the  bottom  of  the  computer 
screen)  waiting  in  separate  locations.  Another  possibility  is  the  subjects  did  not  recognize 
the  idle  loaders  because  they  were  not  truly  “system  experts’*.  LDER  had  a  higher  average 
problem  identification  time  because  no  problem  could  be  seen  in  the  animations  without 
movement  (C,  I,  and  Cl).  The  only  problem  in  LDER  was  a  backup  of  loads,  which  was  not 
shown  in  C,  I,  and  Cl.  Thus,  usually  the  full  60  seconds  would  elapse  without  a  problem 
identification.  The  load  queue  was  not  included  in  the  icon-only  animation  (I)  because 
the  researcher  was  uncertain  how  to  show  the  load  queue  without  movement.  It  was  not 
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shown  in  the  color-only  animation  (C)  because  that  would  have  required  two  different  color 
keys  on  the  screen.  The  reasons  for  not  including  the  load  queue  in  C  and  I  also  apply  to 
Cl.  In  retrospect,  however,  there  should  have  been  some  representation  of  the  load  queue 
in  the  animations  without  movement.  The  average  problem  identification  time  for  SLTK 
was  higher  because  the  backup  of  loads  caused  by  the  one  slow  truck  showed  up  later  in 
the  animation.  Therefore,  the  problem  of  backup  of  loads  could  not  be  identified  as  early 
as  problems  in  the  other  scenarios.  The  problem  of  having  one  slow  truck  could  have 
been  identified  early,  though.  It  if.  interesting  that  very  few  subjects  recognized  one  truck 
being  slower  than  the  others  as  a  problem.  Most  subjects  focused  on  the  loading  area  only 
(for  animations  with  movement)  and  did  not  notice  that  one  truck  was  consistently  being 
passed.  That  is,  one  truck  was  consistently  taking  much  longer  to  haul,  dump,  and  return. 

Therefore,  looking  at  the  data  from  the  scenario  perspective  suggested  several  things. 
First,  animation  may  not  tell  the  whole  story.  Also,  several  animations  should  be  developed 
that  show  the  system  from  different  perspectives.  For  instance,  with  LDER,  the  animations 
without  movement  did  not  show  the  problem  with  the  system  (buildup  of  loads),  and 
with  FL  and  SLTK,  some  areas  of  the  computer  screen  were  not  given  as  much  attention 
as  others.  The  latter  implies  breaking  up  the  system  into  several  views.  That  is,  view 
animations  that  zoom  in  on  certain  areas  as  well  as  view  animations  that  give  an  overview 
of  the  complete  system. 

Comparison  of  Subjective  and  Objective  Results 

For  the  animations  without  movement  (C,  I,  and  Cl)  the  subjective  and  objective 
results  match  well.  These  animation  types  were  rated  as  poor  communicators  and  the 
subjects’  performance  when  viewing  the  animations  without  movement  agreed  with  those 
ratings.  There  was  a  clear  distinction  between  the  animations  with  movement  (M,  MC, 
MI,  and  MCI)  and  the  animations  without  movement  in  both  the  subjective  and  objective 
results.  M,  MC,  MI,  and  MCI  were  much  preferred  over  C,  I,  and  Cl,  and  the  subjects  per¬ 
formed  significantly  better  when  viewing  the  animations  with  movement.  However,  given 
that  there  was  movement  in  the  animation,  the  subjective  and  objective  results  differed. 
The  subjects  preferred  movement  as  the  primary  communication  tool,  but  the  subjects 
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rated  communication  higher  as  color  and  detailed  icons  were  added  to  the  animations. 
MCI  was  selected  79%  of  the  time  as  the  animation  that  communicated  the  best,  and  MCI 
was  rated  the  highest  in  the  pairwise  comparisons.  The  objective  results  showed,  though, 
that  there  was  no  difference  in  subject  performance  for  any  of  the  animations  with  move¬ 
ment.  Therefore,  even  though  the  subjects  preferred  a  more  detailed  animation,  movement 
was  the  primary  factor. 

Table  23  summarizes  the  comparison  of  the  subjective  and  objective  results  and 
summarizes  the  result  of  this  chapter  by  showing  the  animation  groupings  from  the  Tukey 
method  of  multiple  comparisons.  C,  I,  and  Cl  were  collectively  judged  to  communicate 
poorly,  and  the  subjects  did  not  perform  well  when  viewing  these  animation  types.  Al¬ 
though  movement  was  considered  by  the  subjects  to  be  important,  their  preferences  in¬ 
creased  as  color  and  detailed  icons  were  added.  With  regard  to  subject  performance, 
however,  there  was  no  difference  between  the  animations  with  movement  (M,  MI,  MC, 
MCI). 


V.  Conclusions  and  Rax  imendations 

Conclusions 

Summary.  As  stated  in  Chapter  lt  animation  is  a  communication  tool  that  can  help 
in  establishing  a  model’s  face  validity.  Since  communicating  the  operation  of  a  simulation 
model  to  the  system  experts  is  necessary  to  determine  face  validity,  this  research  looked 
at  three  aspects  of  animation  (color,  detail  of  icons,  and  movement)  to  determine  which 
ones  were  the  most  useful  for  communicating  the  operation  of  a  simulation  model.  This 
ability  to  communicate  was  measured  both  subjectively  and  objectively.  The  subjective 
measures  were  a  selection  of  “best”  and  “worst”  animation  types  where  “best”  and  “worst” 
referred  to  how  well  an  animation  communicated,  and  a  pairwise  comparison  of  the  ani¬ 
mation  types  which  resulted  in  preference  ratings  for  each  animation.  There  were  seven 
different  scenarios  containing  various  problems  with  the  system.  The  objective  measures 
were  subject  problem  identification  accuracy  and  time  delay  of  problem  identification. 

Results.  The  results  showed  that  movement  was  the  most  important  aspect  of  anima¬ 
tion.  Animations  with  movement  were  much  preferred  over  animations  withont  movement, 
and  the  subjects  identified  problems  more  accurately  in  less  time  when  viewing  anima¬ 
tions  with  movement  than  animations  without  movement.  Whereas  there  were  differences 
among  the  AHP  ratings  for  M,  MC,  MI,  and  MCI,  there  was  no  difference  among  these 
animations  with  regards  to  subject  performance.  Further  examination  of  the  AHP  ratings 
revealed  that,  given  movement,  detailed  icons  were  preferred  over  color.  Also,  those  who 
liked  detailed  icons  did  not  like  color,  and  those  who  liked  color  did  not  like  detailed  icons. 
However,  most  preferred  having  movement,  color,  and  detailed  icons. 

The  potential  for  bias  in  the  results  should  be  mentioned  here.  Looking  back,  it  is 
possible  that  movement  was  shown  to  be  most  important  because  of  the  scenarios  used. 
Since  the  problems  associated  with  each  scenario  were  primarily  utilization  problems,  they 
might  have  shown  up  easier  with  movement  than  with  color  or  detailed  icons.  Other  types 
of  problems  (such  as  a  truck  leaving  empty)  might  have  been  seen  more  quickly  with 
odor  or  detailed  icons  than  with  movement.  Nevertheless,  some  or  all  of  the  problems 
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associated  with  each  scenario  (except  for  the  LDER  scenario  as  discussed  in  Chapter  4) 
could  be  identified  in  each  animation  type. 

Additional  Observations.  The  simulation  model  used  in  this  experiment  was  chosen 
because  it  was  simple  and  concise.  However,  the  model  was  not  designed  with  animation  in 
mind.  Several  modifications  and  simplifications  were  made  to  the  model  so  that  is  could  be 
animated.  Looking  back,  a  couple  of  the  modifications  might  not  have  been  needed  if  the 
researcher  had  had  more  experience  with  animation  and  with  Proof.  Nevertheless,  some 
model  modifications  would  still  have  been  required.  The  primary  addition  to  the  model  was 
the  ability  to  keep  track  of  the  trucks  and  loaders.  Creating  the  animation  trace  files  would 
have  been  easier  if  each  separate  movement  of  an  entity  had  been  explicitly  modeled.  As 
described  in  Chapter  3,  the  loading  times  had  to  be  artificially  divided  to  account  for  the 
various  movements  required  by  the  loaders.  Therefore,  if  a  modeler  anticipates  animating 
a  model,  this  should  be  kept  in  mind  when  designing  and  coding  the  model.  The  modeler 
should  make  certain,  though,  that  the  system  being  modeled  determines  the  model  design 
and  not  animation  considerations. 

Finally,  even  though  the  initial  aim  of  this  study  was  to  examine  animation’s  role  in 
establishing  face  validity,  the  contribution  of  animation  to  face  validity  was  not  what  was 
actually  measured.  The  simulation  model  that  was  used  was  assumed  to  be  valid;  therefore, 
the  subjects  could  not  judge  the  model’s  face  validity.  In  the  context  of  face  validity,  the 
research  examined  which  aspects  of  animation  best  communicated  the  operation  of  the 
model.  So  the  result  that  using  movement  in  animations  is  important  applies,  not  only 
to  face  validity,  but  also  to  other  validation  and  verification  techniques.  In  addition,  this 
result  applies  to  any  other  areas  in  which  animation  could  be  used,  such  as  communicating 
the  model  to  a  decision  maker. 

Recommendations  for  Further  Study 

Several  aspects  of  animation  that  were  not  considered  in  this  study  (such  as  graphs 
and  speed  of  animation)  were  mentioned  in  Chapter  1.  The  aspects  of  animation  that  were 
not  examined,  plus  what  was  learned  during  this  research  suggest  the  following  studies: 
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•  Repeat  this  research  with  a  larger  simulation  model  of  an  actual  system.  That  is, 
investigate  movement,  color,  and  detail  of  icons  with  a  model  of  a  more  complex, 
real  world  system.  The  model  should  have  more  simultaneous  activities  and  a  larger 
number  of  entities.  Using  this  type  of  model  would  allow  real  “system  experts’*  to  rate 
the  animations.  Also,  the  results  from  this  type  of  study  would  assist  in  determining 
whether  the  results  of  this  study  hold  for  a  more  realistic  scenario. 

•  Investigate  the  use  of  graphs  alone  or  in  combination  with  other  aspects  of  animation. 
This  research  could  address  such  questions  as: 

1.  Is  a  graph  showing  queue  status  necessary  when  the  actual  entities  can  be  seen 
waiting? 

2.  What  information  can  be  displayed  with  graphs  that  can  not  be  displayed  with 
the  aspects  of  animation  investigated  here? 

3.  Is  the  unique  information  displayed  by  graphs  critical? 

4.  Do  graphs  improve  communication  or  become  distractions? 

•  Examine  the  impact  of  the  viewing  speed  of  the  animations.  Many  times  during 
the  viewing  of  the  animations  without  movement,  the  subjects  commented  that  the 
icons  or  colors  were  changing  too  fast.  How  would  have  the  subjects  performed  if  the 
animations  without  movement  had  been  slower?  Is  there  an  optimal  viewing  speed, 
and  can  it  be  determined? 

•  Study  the  usefulness  of  color  in  communicating  when  the  colors  have  well  known 
meanings  in  the  context  of  the  system  being  modeled.  For  instance,  in  this  study, 
the  colors  were  assigned  arbitrarily.  White  represented  an  idle  truck  or  loader,  and 
red  represented  a  loaded  truck  or  loader.  These  were  subjective  color  assignments. 
What  if,  in  the  system  modeled,  red  meant  stop  (as  on  a  traffic  light  or  stop  sign), 
or  red  meant  hot?  When  there  is  meaning  to  the  colors,  the  importance  of  color  as 
a  communication  tool  might  increase. 


Appendix  A.  Data  Collection  Forms 

This  appendix  shows  the  data  collection  forms  used.  Page  61  displays  the  form  osed 
to  document  which  scenario  and  animation  was  viewed,  the  time  required  to  identify  the 
problem,  and  the  problem  identified.  The  pairwise  comparison  data  collection  instructions, 
examples,  and  form  follow  on  page  62  through  page  64.  Page  65  shows  the  pairwise  matrix 
form.  The  model  description  ends  this  appendix  and  begin  on  page  66. 


i 

! 

! 


60 


Scenario/ Animation  Viewing 

Name: -  Date: -  Run: 

Student  _  Faculty  _  Other  — 

Scenario: _  Animation: _  Problem  ID  Time: _ 

Problem  Observed: _ 


Scenario: _  Animation: - 

Problem  ID  Time: _ 

Problem  Observed: 

Scenario: _ 

Animation: - 

Problem  ID  Time: _ 

Problem  Observed: 

Scenario: - 

Problem  Observed: 

Animation: - 

Problem  ED  Time: - 

Scenario: _ 

Problem  Observed: 

Animation: - 

Problem  ID  Time: _ 

Scenario: _ 

Problem  Observed: 

Animation: _ 

Problem  ID  Time: _ 

Scenario: _ 

Problem  Observed: 

Animation: - 

Problem  ID  Time: _ 

Best  Animation  for  Face  Validity:  . 
Worst  Animation  for  Face  Validity: 
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Animation  Pairwise  Comparison 


The  pairwise  comparison  of  the  animations  allows  yon  to  rate  an  animation’s  ability  to 
communicate  the  model’s  operation  against  the  ability  of  another  animation.  This  allows 
the  relative  merits  of  the  animations  to  be  quantified.  Since  there  are  seven  animations,  yon 
will  perform  21  pairwise  comparisons.  After  you  have  completed  the  first  few  comparisons, 
I  will  briefly  interrupt  you.  This  is  so  I  can  relate  back  to  you  one  or  two  of  your  responses 
to  verify  that  I  have  communicated  the  instructions  properly.  Feel  free  to  go  back  and 
erase  a  previous  response  if  you  think  it  needs  adjusting.  The  order  of  your  responses  is 
not  important.  What  I  need  is  a  completed  form  that  reflects  your  best  comparison  of  the 
animations. 

As  a  reminder,  the  animation  types  are: 

•  M  -  Movement.  Simple  icons  that  move  but  do  not  change  level  of  detail  or  color. 

•  I  -  Icon.  Icons  change  level  of  detail  but  do  not  move  or  change  odor. 

•  C  -  Color.  Simple  icons  that  change  color  but  do  not  move  or  change  level  of  detail. 

•  MI  -  Movement  and  Icon.  Icons  move  and  change  level  of  detail  but  do  not  change 
color. 

•  MC  •  Movement  and  Color.  Simple  icons  that  move  and  change  color  but  do  not 
change  level  of  detail. 

•  Cl  -  Color  and  Icon.  Icons  change  color  and  level  of  detail  but  do  not  move. 

•  MCI  •  Movement,  Color,  and  Icon.  Icons  move  and  change  color  and  level  of  detail. 

The  pictures  of  each  type  of  animation  are  available  for  you  to  look  at  while  you  are 
performing  the  comparisons.  The  next  page  gives  some  pairwise  comparison  examples, 
and  the  final  page  is  the  pairwise  comparison  form. 
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Pairwise  Comparison  Examples 


Animation  A  and  Animation  B  have  equal  contribution  to  face  validity. 


Animation  A  contributes  slightly  more  to  face  validity  than  Animation  B. 


»»  »>  »  >  as  <  «  «<  «« 


Animation  A  contributes  somewhat  more  to  face  validity  than  Animation  B. 


Animation  A  contributes  much  more  to  face  validity  than  Animation  B. 


Animation  B  contributes  slightly  more  to  face  validity  than  Animation  A. 


»»  »>  »  >  =  <  «  «<  «« 
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A  Thick  Hauling  Situation 

The  system  modeled  consists  of  one  bulldozer,  four  trucks,  and  two  man-machine  loaders. 
The  bulldozer  stockpiles  material  for  the  loaders.  Two  piles  of  material  must  be  stocked 
prior  to  the  initiation  of  any  load  operation.  In  addition  to  the  two  loads  of  material,  a 
loader  and  an  unloaded  truck  must  be  available  before  the  loading  operation  can  begin. 
There  is  room  in  the  loading  area  for  only  one  truck  and  one  loader.  After  a  truck  is 
loaded,  it  is  hauled  and  then  dumped.  It  must  be  returned  before  the  truck  is  available 
for  further  loading.  Also,  the  dumping  area  can  only  accommodate  one  truck.  Following  a 
loading  operation,  the  loader  must  return  to  its  idle  position  before  it  is  available  to  begin 
loading  again.  The  diagram  below  shows  the  flow  of  the  system. 


Truck  Return 


Loader  Return 


You,  as  the  owner  of  the  trucks,  loaders,  and  bulldozer,  requested  a  simulation  to  study 
the  efficiency  of  your  operation.  My  job  as  the  simulation  model  builder  is  to  communicate 
to  you  the  operation  of  the  model  so  that  you  can  determine  whether  the  model  accurately 
reflects  your  application.  This  procedure  is  establishing  the  face  validity  of  the  model. 
I  have  decided  to  use  animation  to  relate  how  your  system  is  modeled.  You  will  view 
seven  different  animations  of  the  model.  While  viewing  the  animations,  try  to  determine 
whether  there  is  a  problem  with  the  system.  Any  particular  animation  may  or  may  not 
exhibit  a  problem.  Problems  could  include  over  or  under  utilization  of  the  trucks,  loaders, 
or  bulldozer;  a  large  build-up  of  loads;  an  inadequate  number  of  trucks,  loaders,  and 


bulldozers;  or  a  combination  of  these.  After  all  the  animations  have  been  viewed,  you 
will  perform  a  pairwise  comparison  of  each  type  of  animation.  That  is,  you  will  compare 
how  each  type  of  animation  communicates  the  operation  of  the  system  (face  validity)  with 
every  other  type  of  animation. 

The  types  of  animation  are: 

•  M  -  Movement,  simple  icons  that  move  but  do  not  change  level  of  detail  or  color. 

•  I  -  Icon.  Icons  change  in  level  of  detail  but  do  not  move  or  change  color. 

•  C  •  Color.  Simple  icons  that  change  color  but  do  not  move  or  change  in  level  of 
detail. 

•  MI  -  Movement  and  Icon.  Icons  move  and  change  in  level  of  detail  but  do  not  change 
color. 

•  MC  •  Movement  and  Color.  Simple  icons  that  move  and  change  color  but  do  not 
change  in  level  of  detail. 

•  Cl  -  Color  and  Icon.  Icons  change  color  and  level  of  detail  but  do  not  move. 

•  MCI  -  Movement,  Color,  and  Icon.  Icons  move,  change  color,  and  change  level  of 
detail. 

The  above  abbreviations  will  be  used  to  identify  each  animation.  Using  pictures  I  will 
describe  the  layout  of  each  animation.  When  viewing  the  animations,  speak  out  as  soon 
as  you  think  you  have  discovered  a  problem  (such  as  saying  “I  got  it”).  Feel  free  to  ask 
questions  before  you  begin  viewing  the  animations.  Let  me  know  when  you  are  ready  to 
start. 


/ 
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Appendix  B.  SAS  Output 

This  appendix  contains  the  SAS  output  nsed  in  the  analysis.  The  SAS  command  hie 
begins  on  page  69,  and  the  output  listing  begins  on  page  71.  The  data  sets  used  are: 

•  GMEANS  -  AHP  ratings  for  each  animation  type. 

•  GMEANST  -  AHP  ratings  transposed.  That  is,  instead  of  the  ratings  for  each 
animation  type  being  in  columns,  the  ratings  are  in  rows.  This  transformation  was 
necessary  for  the  cluster  analysis. 

•  PROBS  •  Problem  data  from  animation  perspective.  Each  animation  type  had  an 
accuracy  variable  and  a  time  variable.  For  instance,  for  M,  MP  contained  the  prob¬ 
lem  accuracy  data  (0  for  incorrect,  1  for  correct),  and  MT  contained  the  problem 
identification  times. 

•  RPROBS  •  Problem  data  for  the  correct  identifications  only.  Dl  through  D7  were 
used  for  the  problem  identification  accuracy  data  because  that  data  was  not  needed 
(all  were  1). 

•  WPROBS  •  Problem  data  for  the  incorrect  identifications  only.  D8  through  D14  were 
used  for  the  problem  identification  accuracy  data  because  that  data  was  not  needed 
(all  were  0). 

■  INTIME  -  Problem  data  *n  which  the  subject  responded  before  60  seconds  had 
elapsed.  For  example,  MNA  contains  the  movement  animation  (M)  problem  ac¬ 
curacy  data,  and  MNT  contains  the  problem  time  data  for  M. 

•  OUTTIME  -  Problem  data  in  which  the  subject  responded  before  60  seconds  had 
elapsed.  D15  through  D21  were  used  for  the  problem  identification  time  data  because 
that  data  was  not  needed  (all  were  60). 

•  PROBST  -  Problem  identification  times  transposed  for  cluster  analysis. 

•  SCENARIO  -  Problem  data  from  scenario  perspective.  Each  scenario  had  an  accu¬ 
racy  variable  and  a  time  variable.  The  variables  follow  the  same  pattern  as  PR03S. 
For  example,  LDERP  contained  the  problem  accuracy  data  for  the  LDER  scenario, 
and  LDERT  contained  the  problem  identification  times  for  the  LDER  scenario. 
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OPTIOIS  LIBESIZE-78; 

DATA  GMEAHS ; 

IIPUT  M  C  I  MI  MC  Cl  MCI; 

XIICLUDE  MEAIS; 

DATA  GXEAIST; 

IIPUT  P1-P47; 

7IICLUDE  MEAIST; 

DATA  PROBS; 

IIPUT  MP  MT  CP  CT  IP  IT  MIP  MIT  MCP  MCT  CIP  CIT  MCIP  MCIT; 

XIICLUDE  PROBS; 

DATA  RPROBS; 

IIPUT  D1  MRP  D2  CRP  D3  IRP  D4  MIRP  DS  MCRP  D6  CIRP  D7  MCIRP; 

XIICLUDE  RPROBS; 

DATA  WPRQBS; 

IIPUT  D8  MVP  D9  CWP  DIO  IIP  Dll  MIVP  D12  MCVP  D13  CIVP  D14  MCIVP; 

XIICLUDE  HPROBS; 

DATA  IITIME; 

IIPUT  MIA  MIT  CIA  CRT  ISA  I  IT  MIIA  MIIT  MCMA  MOTT  CINA  CIBT  MCIIA  MCIIT; 
XIICLUDE  IITIME; 

DATA  OUTTIME; 

IIPUT  MOA  D15  COA  D16  IOA  D17  MIOA  D18  MCOA  D19  CIOA  D20  MCIOA  D21; 
XIICLUDE  OUTTIME; 

DATA  PROBST; 

IIPUT  T1-T47; 

XIICLUDE  PROBST; 

DATA  SCEIARIO; 

IIPUT  LDERP  LDERT  SLTKP  SLTKT  SLLDP  SLLDT  FTP  FTT  FLP  FLT  SLP  SLT  SI.'  STT; 
XIICLUDE  SCEIDAT; 

PROC  PRIIT  DATA-GMEAIS; 

PROC  MEAIS  DATA-GHEAIS; 

PMC  FACTOR  DATA-GMEAIS  IFACTORS-3  METHOD-PRIM  COV  ROTATE- VARI MAI 
OUT-SCORES; 

VAR  M  C  I  MI  MC  Cl  MCI; 

PMC  PRIIT  DATA-SCORES; 

VAR  FACTOR 1  FACT0R2  FACTOR3; 

PMC  PLOT  DATA-SCORES;  ; 

PLOT  FACT0R2*FACT0R1  FACTO R3-FACTOR1  FACT0R3*FACT0R2 ; 

PROC  CLUSTER  DATA-CHEAIST  METHOD-COMPLETE; 

PMC  TREE; 

PMC  CLUSTER  DATA-GMEAIST  METHOD-AVERAGE  PSEUDO; 

PROC  TREE; 

PROC  CLUSTER  DATA-GMEAHST  METHOD-HARD  PSEUDO; 

PMC  TREE; 

PROC  PRIIT  DATA-PMBS; 

VAR  MP  MT  CP  CT  IP  IT  MIP  MIT  MCp  MCT  CIP  CIT  MCIP  MCIT; 


PROC  HEARS  DATA-PROBS ; 

PROC  PRIIT  DATA-RPROBS; 

VAR  HRP  CRP  IRP  MIRP  MCRP  CIRP  HCIRP; 

PROC  HEARS  DATA-RPROBS; 

VAR  HRP  CRP  IRP  HIRP  MCRP  CIRP  HCIRP; 

PROC  PRIIT  DATA-VPROBS; 

VAR  MVP  CVP  IVP  HIWP  HCVP  CIWP  MCI VP; 

PROC  HEARS  DATA-HPROBS; 

VAR  MVP  CVP  IMP  MIVP  HCVP  CIWP  MCI VP; 

PROC  PRIIT  DATA-IHTIME; 

PROC  HEARS  DATA-IHTIME; 

PROC  PRIIT  DATA-OUTTIME; 

VAR  MOA  COA  IOA  MIOA  HCOA  CIOA  HCIOA; 

PROC  HEARS  DATA-OUTTIME; 

VAR  MOA  COA  IOA  HIOA  HCOA  CIOA  HCIOA; 

PROC  CLUSTER  DATA-PROBST  METHOD-COMPLETE; 

PROC  TREE; 

PROC  CLUSTER  DATA-PROBST  METHOD-AVEKASE  PSEUDO 
PROC  TREE; 

PROC  CLUSTER  DATA-PROBST  METHOD-HARD  PSEUDO; 
PROC  TREE; 

PROC  PRIIT  DATA-SCEHARIO; 

PR.C  HEARS  DATA-SCEHARIO; 


/  • 
/ 
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Th*  SIS  Systaa 

10:64  Thursday,  Decaabcr  17, 


1 

1992 


OfiS 

H 

C 

I 

MX 

CX 

HCI 

1 

0.072 

0.020 

0.033 

0.212 

0.200 

0.064 

0.399 

2 

0.072 

0.018 

0.076 

0.243 

0.197 

0.068 

0.326 

3 

0.206 

0.019 

0.028 

0.276 

0.168 

0.031 

0.276 

4 

0.104 

0.020 

0.033 

0.246 

0.162 

0.043 

0.402 

6 

0.091 

0.020 

0.061 

0.214 

0.180 

0.073 

0.371 

6 

0.117 

0.019 

0.038 

0.273 

0.111 

0.058 

0.384 

7 

0.109 

0.017 

0.028 

0.326 

0.084 

0.039 

0.397 

8 

0.173 

0.038 

0.024 

0.234 

0.137 

0.067 

0.338 

9 

0.037 

0.038 

0.039 

0.130 

0.231 

0.102 

0.426 

10 

0.111 

0.022 

0.030 

0.237 

0.141 

0.047 

0.413 

11 

0.047 

0.016 

0.033 

0.176 

0.159 

0.085 

0.484 

12 

0.119 

0.042 

0.019 

0.194 

0.213 

0.071 

0.342 

13 

0.074 

0.016 

0.024 

0.248 

0.181 

0.039 

0.418 

14 

0.166 

0.012 

0.036 

0.421 

0.083 

0.024 

0.267 

16 

0.136 

0.100 

0.040 

0.129 

0.166 

0.124 

0.306 

16 

0.266 

0.030 

0.020 

0.126 

0.139 

0.144 

0.287 

17 

0.117 

0.026 

0.063 

0.291 

0.148 

0.034 

0.321 

18 

0.128 

0.021 

0.048 

0.246 

0.112 

0.041 

0.405 

19 

0.116 

0.029 

0.038 

0.243 

0.174 

0.063 

0.347 

20 

0.088 

0.026 

0.026 

0.202 

0.168 

0.039 

0.463 

21 

0.145 

0.026 

0.06S 

0.244 

0.264 

0.066 

0.221 

22 

0.164 

0.021 

0.029 

0.223 

0.223 

0.028 

0.323 

23 

0.112 

0.022 

0.033 

0.336 

0.164 

0.048 

0.297 

24 

0.116 

0.022 

0.033 

0.287 

0.168 

0.039 

0.346 

25 

0.116 

0.016 

0.028 

0.278 

0.138 

0.044 

0.383 

26 

0.147 

0.020 

0.026 

0.341 

0.138 

0.033 

0.296 

27 

0.348 

0.017 

0.036 

0.241 

0.093 

0.034 

0.232 

28 

0.i26 

0.019 

0.030 

0.171 

0.240 

0.041 

0.374 

29 

0.087 

0.015 

0.038 

0.296 

0.137 

0.054 

0.373 

30 

0.107 

0.036 

0.016 

0.219 

0.267 

0.031 

0.326 

31 

0.131 

0.025 

0.037 

0.148 

0.16.. 

0.064 

0.434 

32 

0.096 

0.029 

0.066 

0.286 

0.164 

0.070 

0.301 

33 

0.033 

0.064 

0.021 

0.260 

0.238 

0.061 

0.332 

34 

0.072 

0.114 

0.024 

0.066 

0-272 

0.093 

0.360 

36 

0.078 

0.023 

0.029 

0.164 

0.165 

0.078 

0.473 

36 

0.106 

0.046 

0.039 

0.213 

0.162 

*5.068 

0.378 

37 

0.142 

0.047 

0.040 

0.242 

0.131 

0.058 

0.339 

38 

0.114 

0.018 

0.030 

b.260 

0.160 

0.043 

0.396 

39 

0.166 

0.023 

0.023 

0.360 

0.104 

0.046 

0.288 

40 

0.088 

0.026 

0.022 

0.131 

0.226 

0.043 

0.465 

41 

0.070 

0.042 

0.026 

0.162 

0.262 

0.063 

0.386 

42 

0.123 

0.016 

0.034 

01280 

0.148 

0.041 

0.369 

43 

0.153 

0.046 

0.062 

0.282 

0.139 

0.082 

0.236 

44 

0.097 

0.020 

0.044 

0.164 

0.193 

0.067 

0.416 

45 

0.099 

0.020 

0.038 

0.291 

0.144 

0.054 

0.364 

46 

0.280 

0.024 

0.049 

0,|l79 

0.216 

0.036 

0.216 

47 

0.123 

0.019 

0.030 

0.246 

0.156 

0.037 

0.389 
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The  SAS  System  3 

10:64  Thursday,  December  17,  1993 

Initial  Factor  Method:  Principal  Components 

Prior  Commonality  Estimates :  OXE 

Eigenvalues  o<  tha  Covariance  Matrix: 

Total  -  0.01687276  Average  »  0.00226764 


1 

2 

3 

4 

Eigenvalue 

0.0081 

0.0046 

0.0021 

0.0008 

Difference 

0.0035 

0.0026 

0.0013 

0.0006 

Proportion 

0.6088 

0.2872 

0.1325 

0.0501 

Cumulative 

0.5088 

0.7960 

0.9285 

0.9787 

6 

6 

7 

Eigenvalue 

0.0002 

0.0001 

0.0000 

Difference 

0.0001 

0.0001 

Proportion 

0.0131 

0.0082 

0.0000 

Cumulative 

0.9918 

1.0000 

1.0000 

3  factors  sill  be  retained  by  the  (FACTOR  criterion. 
Factor  Pattern 


FACT081 

FACT0R2 

FACTORS 

a 

-0.67140 

0.66869 

0.44361 

c 

0.30724 

0.47428 

-0.23686 

i 

-0.21608 

0.02492 

-0.18149 

HI 

-0.78181 

-0.69349 

-0.18454 

MC 

0.67950 

0.47349 

-0.60049 

Cl 

0.39917 

0.39295 

0.09424 

NCI 

0.79631 

-0.60630 

0.31687 

Variance  explained  by  each  factor 

FACTOR!  FACTORS  FACT0R3 
Weighted  0.008076  0.004659  0.002104 
Ohveighted  2  'TT774  1.636019  0.789035 

Final  Co— nnallty  Estimates  and  Variable  Weights 
Total  Commonality:  Weighted  -  0.014738  Obveighted  -  4.654823 

«  C  X  HI 
Commonality  0.970776  0.375436  0.080263  0.997622 
Weight  0.003314  0.000378  0.000168  0.004808 

ac  a  mci 
Commonality  0.920596  0.322632  0.987616 
Weight  0.002286  0.300589  0.004329 
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Rotation  Method:  Vaxisax 


11m  SIS  Sjstaa  4 

10:64  Thursday,  Dacaabar  17,  1992 


Orthogonal  Transformation  Matrix 


12  3 


1 

0.64283 

-0.62861 

0.65262 

2 

0.77336 

0.61765 

-0.14296 

3 

-0.32762 

0.68230 

0.74408 

Rotated  Factor  Pattern 

FACT0R1 

FACTORS 

FACT0R3 

N 

-0.07000 

0.96436 

-0.18944 

C 

0.61113 

-0.00739 

-0.04364 

1 

-0.03868 

0.02393 

-0.27963 

MI 

-0.82293 

-0.06076 

-0.56269 

MC 

0.87741 

-0.36356 

-0.13631 

Cl 

0.48971 

0.08667 

0.27446 

MCI 

-0.06262 

-0.64868 

0.82630 

Variance  explained  by  each  factor 

FACTOR!  FACT0R2  FACT0R3 
Weighted  0.006332  0.004709  0.004697 
Unweighted  2.070667  1.374910  1.209262 

Final  Commality  Eatiaataa  and  Variable  Weights 


Total  Commality: 

Weighted 

-  0.014738 

Unweighted  -  4.664828 

H 

C 

X 

MI 

Coomnality 

0.970776 

0.376436 

0.080253 

0.997622 

Weight 

0.003314 

0.000378 

0.000168 

0.004808 

MC 

a 

MCI 

Canswnality 

0.920696 

0.322632 

0.987615 

Weight 

0.002286 

0.000689 

0.004329 

Scaring  Coefficients  Estinsted  by  Regression 
Squared  Multiple  Correlations  of  the  Variablu  with  each  Factor 


FACT0R1  FACTORS  FACT0R3 

1.000000  1.000000  1.000000 
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Tbi  SIS  Sjritn 


10:64  Thursday,  December  17,  1992 


Rotation  Method:  Varinax 

Standardized  Scoring  Coeff  iciants 


FACT0R1 

FACT0R2 

FACT0R3 

R 

-0.06874 

0.80778 

0.28096 

C 

0.06222  -0.00810 

-0.02794 

I 

0.00302  -0.00551 

-0.01389 

MI 

-0.69868  -0.38613 

-0.62808 

HC 

0.48634  -0.32001 

-0.41238 

Cl 

0.04644 

0.03134 

0.03139 

MCI 

-0.36251  -0.14324 

0.83040 
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Ths  SAS  Sjstm  6 

10:64  Thursday,  Dscnbsr  17,  1992 


088 

FACT0R1 

FACT0R2 

FACT0R3 

1 

0.30699 

-0.86848 

0.24602 

2 

0.41129 

-0.87436 

-0.92362 

3 

-0.10170 

1.09098 

-0.90960 

4 

-0.66469 

-0.32118 

0.66797 

6 

0.23733 

-0.41422 

0.13600 

6 

-1.10447 

0.04215 

0.66772 

7 

-1.94219 

-0.23229 

0.51099 

8 

-0.26817 

0.96796 

0.31962 

9 

1.34803 

-1.12803 

0.77736 

10 

-0.66264 

-0.11766 

0.90990 

11 

-0.20344 

-0.90023 

1.86896 

12 

0.92024 

-0.06685 

-0.22662 

13 

-0.36267 

-0.98188 

0.36834 

14 

-2.16416 

0.16721 

-1.63606 

16 

1.44181 

0.94893 

0.26029 

16 

1.01814 

2.93499 

1.00966 

17 

-0.66993 

-0.21268 

-0.75004 

18 

-1.00116 

0.27298 

1.06662 

19 

0.01189 

-0.15666 

-0.23989 

20 

-0.42773 

-0.47789 

1.53791 

21 

1.46968 

-0.01932 

-2.39613 

22 

0.70139 

0.18718 

-0.63886 

23 

-0.74331 

-0.48370 

-1.40970 

24 

-0.67107 

-0.31938 

-0.46719 

26 

-0.90644 

-0.20446 

0.27272 

26 

-1.02391 

0.06770 

-1.16823 

27 

-0.48460 

3.88417 

0.16069 

28 

1.09930 

-0.13779 

0.13130 

29 

-0.96876 

-0.65489 

-0.10876 

30 

1.26348 

-0.74632 

-1.19473 

31 

0.26684 

0.50128 

1.79467 

32 

-0.24719 

-0.44104 

-1.07830 

33 

0.76966 

-1.80443 

-1.52314 

34 

2.82205 

-0.44477 

0.16445 

36 

0.08196 

-0.36963 

1.97096 

36 

0.04668 

-0.13166 

0.41209 

37 

-0.34220 

0.61716 

0.14301 

38 

-0.69833 

-0.17240 

0.62130 

39 

-1.46788 

0.34997 

-1.06101 

40 

0.87046 

-0.63616 

1.52719 

41 

1.68696 

-0.95398 

-0.04268 

42 

-0.70490 

-0.12446 

-0.10293 

43 

-0.01177 

0.64370 

-1.48063 

44 

0.64622 

-0.23917 

0.98802 

46 

-0.75006 

-0.47180 

-0.32710 

46 

1.48278 

2.48068 

-0.99428 

47 

-0.49877 

-0.06241 

0.46764 
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Tho  SIS  Syatui  7 

10:64  Thursday,  Dtcubn  17.  1992 
Plot  ot  FlCT0R2»FiCT0ai .  Lgend:  A  -  1  ot»,  B  -  2  ot»,  «tc. 

FACTORS  | 


4  ♦ 


A 


PACT0R1 


P1CT0R3 


Thm  SAS  Sjrttn  $ 

10:54  Thursday,  D*c«mb«r  17,  1992 
Plot  of  FACT0A3*PACT0R1 .  Lagand:  A  -  1  obo,  B  ■  2  obs,  otc. 
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0.0 


>0.6 
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Tha  3 AS  Syatas  10 

10:54  Thursday,  Dacasbar  17,  1992 
Coaplat*  Liakaga  Clustar  Analysis 


Eigantalnas  of  tha  Covarianca  Matrix 


Eigaatalna 

Oiffaranca 

Proportion 

Cumlativa 

1 

0.697871 

0.638318 

0.874651 

0.87465 

2 

0.059553 

0.037518 

0.074639 

0.94929 

3 

0.022035 

0.007156 

0.027617 

0.97691 

4 

0.014879 

0.012880 

0.018648 

0.99555 

6 

0.001999 

0.000462 

0.002606 

0.99806 

6 

0.001547 

0.001547 

0.001939 

1.00000 

7 

0.000000 

0.000000 

0.000000 

1.00000 

8 

0.000000 

0.000000 

0.000000 

1.00000 

9 

0.000000 

0.000000 

0.000000 

1.00000 

10 

0.000000 

0.000000 

0.000000 

1.00000 

11 

0.000000 

0.000000 

0.000000 

1.00000 

12 

0.000000 

0.000000 

0.000000 

1.00000 

13 

0.000000 

0.000000 

0.000000 

1.00000 

14 

0.000000 

0.000000 

0.000000 

1.00000 

16 

0.000000 

0.000000 

0.000000 

1.00000 

16 

0.000000 

0.000000 

0.000000 

1.00000 

17 

0.000000 

C.OOOOOO 

0.000000 

1.00000 

18 

0.000000 

0.000000 

0.000000 

1.00000 

19 

0.000000 

0.000000 

0.000000 

1.00000 

20 

0.000000 

0.000000 

0.000000 

1.00000 

21 

0.000000 

0.000000 

0.000000 

1.00000 

22 

0.000000 

0.000000 

0.000000 

1.00000 

23 

0.000000 

0.000000 

0.000000 

1.00000 

24 

0.000000 

0.000000 

0.000000 

1.00000 

25 

0.000000 

0.000000 

0.000000 

1.00000 

26 

0.000000 

0.000000 

0.000000 

1.00000 

27 

>.000000 

0.000000 

-.000000 

1.00000 

28 

-.000000 

0.000000 

-.000000 

1.00000 

29 

-.000000 

0.000000 

-.000000 

1.00000 

30 

-.000000 

0.000000 

-.000000 

1.00000 

31 

-.000000 

0.000000 

-.000000 

1.00000 

32 

-.000000 

0.000000 

-.000000 

1.00000 

33 

-.000000 

0.000000 

-.000000 

1.00000 

34 

-.000000 

0.000000 

-.000000 

1.00000 

35 

-.000000 

0.000000 

-.000000 

1.00000 

36 

-.000000 

0.000000 

-.000000 

1.00000 

37 

-.000000 

0.000000 

-.000000 

1.00000 

38 

-.000000 

0.000000 

-.000000 

1.00000 

39 

-.900000 

0.000000 

-.000000 

1.00000 

40 

-.000000 

0.000000 

-.000000 

1.00000 

41 

-.000000 

0.000000 

-.000000 

1.00000 

42 

-.000000 

0.000000 

-.000000 

1.90000 

43 

-.000000 

0.000000 

-.000000 

1.00000 

44 

-.000000 

0.000000 

-.000000 

1.00000 

45 

-.000000 

0.000000 

-.000000 

1.00000 

46 

-.000000 

0.000000 

-.000000 

1.00009 

47 

-.000000 

-.000000 

1.00000 

Kcnt-Haan-Squara  Total-Saspla  Standard  Satiation  -  0.130293 
Mean  Distanco  Batwaan  Obsanrations  ■  1.107356 
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The  SAS  Systea  11 

10:54  Thursday,  Decenber  17,  1992 


Coaplete  Linkage  Cluster  Analysis 


funber 

et 

Cl us tars 

Clusters  Joined 

Frequency 

oi  lee 

Cluster 

loraalized 

Marietta 

Distance 

C 

082 

083 

2 

0.154549 

6 

CL6 

□B6 

3 

0.214065 

4 

081 

OBS 

2 

0.604346 

3 

CL4 

084 

3 

0.866503 

2 

CL3 

CLfi 

6 

1.364139 

1 

CL2 

087 

7 

2.0628C8 
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The  SIS  System  12 

10:64  Thursday,  December  17,  1992 
Complete  Linkage  Cluster  Analysis 

Isas  of  Observation  or  Cluster 
OB  OB  OB  OB  OB  OB  OB 

1  6  4  2  3  6  7 
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Tha  SIS  Syatan  13 

10tS4  Thursday,  December  17,  1992 
Average  Linkage  Cluster  Analysis 


Eigenvalues  of  ths  Covariance  Matrix 


Eigenvalue 

Diffaranca 

Proportion 

Cumulative 

1 

0.697871 

0.638318 

0.874651 

0.87465 

2 

C. 069563 

0.037618 

0.074639 

0.94929 

3 

0.022036 

0.007156 

0.027617 

0.97*91 

4 

0.014879 

0.012880 

0.018648 

0.99555 

6 

0.001999 

0.000452 

0.002606 

0.99806 

6 

0.001547 

0.001547 

0.001939 

1.00000 

7 

0.000000 

0.000000 

0.000000 

1.00000 

8 

0.000000 

0.000000 

0.000000 

1.00000 

• 

0.000000 

0.000000 

0.000000 

1.00000 

10 

0.000000 

0.000000 

0.000000 

1.00000 

11 

0.000000 

0.000000 

0.000000 

1.00000 

12 

0.000000 

0.000000 

0.000000 

1.00000 

13 

0.000000 

o.ocoooo 

0.000000 

1.00000 

14 

0.000000 

0.000000 

0.000000 

1.00000 

16 

0.000000 

0.000000 

0.000000 

1.00000 

16 

n. 000000 

0.000000 

0.000000 

1.00000 

17 

0.000300 

0.000000 

0.000000 

1.00000 

18 

0.000000 

0.000000 

0.000000 

1.00000 

19 

0.000000 

0.000000 

0.000000 

1.00000 

20 

0.000000 

0.000000 

0.000000 

1.00000 

21 

0.000000 

0.000000 

0.000000 

1.00000 

22 

0.000000 

0.000000 

0.000000 

1.00000 

23 

0.000000 

0.000000 

0.000000 

1.00000 

24 

0.000000 

0.000000 

0.000000 

1.00000 

26 

0.000000 

0.000000 

0.000000 

1.00000 

26 

0.000000 

0.000000 

0.000000 

1.00000 

27 

-.000000 

0.000000 

-.000000 

1.00000 

28 

-.000000 

0.000000 

-.000000 

1.00000 

29 

-.000000 

0.000000 

-.000000 

1.00000 

30 

-.000000 

( . 000000 

-.000000 

1.00000 

31 

-.000000 

*\  000000 

-.000000 

1.00000 

32 

-.000000 

0.000000 

-.000000 

1.00000 

33 

-.000000 

0.000000 

-.000000 

1.00000 

34 

-.000000 

0.000000 

-.000000 

1.00000 

36 

-.000000 

0.000000 

-.000000 

1.00000 

36 

-.000000 

0.000000 

-.000000 

1.00000 

37 

-.000000 

0.000000 

-.000000 

1.00000 

38 

-.000000 

0.000000 

-.000000 

1.00000 

39 

-.000000 

0.000000 

-.000000 

1.00000 

40 

-.000000 

0.000000 

-.000000 

1.00000 

41 

-.000000 

0.000000 

-.000000 

1.00000 

42 

-.000000 

0.000000 

-.000000 

1.00000 

43 

-.000000 

0.000000 

-.000000 

1.00000 

44 

-.000000 

0.000000 

-.000000 

1.00000 

46 

-.000000 

0.000000 

-.000000 

1.00000 

46 

-.000000 

0.000000 

-.000000 

1.00000 

47 

-.000000 

# 

-.000000 

1.00000 

Root-Mean-Squara  Total-Saapla  Standard  Deviation  "  0.130293 
Root-Mean-Squara  Oistanca  Bataaan  Observation*  “  1.263238 
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The  SA 8  Sjatea  14 

10:64  Thuraday,  Dacambar  17,  1‘12 


Average  Linkage  Cluatar  Analyaia 
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of 
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Cluatars  Joined 
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r 
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2 

•6.18 
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6 

CL6 

0B6 

3 
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2.14 
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4 

OBI 
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2 

16.74 

. 
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3 

CU 

OS 

6 

7.99 

7.65 

0.662877 

2 

0B4 

0B7 

2 

10.07 

a 

0.888342 

1 

03 

02 

7 

. 

10.07 

1.311767 
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Tha  SiS  3y»t«« 

10:54  Thursday,  Dacanbar  17, 
Aaaraga  Liniaga  Clustar  Analysis 
■aaa  of  Obaarration  or  Clustar 


15 

1992 


OB 

1 


OB 

5 


OB 

2 


OB 

3 


6 


A 

a 

a 

r 

a 

S 


D 

i 

s 

t 

a 

n 

c 


1.4  ♦ 

I 

I 

II 

I  nmmnmnnramn 

|  TTTIXTTTniTtrrtTTTlITIT-r 

I  mnnmmnmmim 
1.2  +xmxxxxxmxxxxxxnxmx 
l  immmmnmmnii 
imnummxmuimn 
I  mnnnmmnamm 
I  miminmnnimnn 
I  mnnmmnmnirm 

|  TYTTTTTTTTITITITTTTIITTrf 

1  ♦TTTTITTTITTTmiTTniira 

|  nrmmimnmmnn 

I  innmnnmmmaa 
I  imininmnnininn 

[  TtTTTTTTTTIITlimllinri 

I  mnnmmnmnian 
o.8  ^mmjpmmninmnn 
i  mmnnmnnimnu 
I  mnmnmnnimnn 

I  iiiiiiiiniiiiiiiiniiiii 
t  Tmmnmnmmnn 


I 

I 

I 

0  ♦ 


t  lixum 

xxxxxmxxm 

v  o.6  +xxxxxxx 

a  I  xxxxm 

xmxxmxm 

a  mum 

xxxxxmxxm 

n  1. 

1. 

mxxmmn 

C  1. 

xxxxxmxxm 

1  1. 

mnmmn 

u  0.4  +. 

imimxxm 

s  1 . 

t  1. 

•  1* 

xmxxxxxxxn 

r  1  a 

xxmmxxxxx 

•  |a 

1. 

xxmmixm 

0.2  ♦. 

xxxmmxxxx 

1. 

xxxxm 

OB 

4 


QB 

7 


ixxxra 

nxxm 

mini 

mini 

innn 

ITTIIII 

xxxxm 

mini 

xxxxm 

xxxxm 

mnn 

imxxx 

xxxxm 

mini 

xxxxm 


85 


T&tJtLh  -  .&■ 1- 


The  8AS  System  16 

10:64  Thursday,  December  17,  1992 
Yard’s  Hiaisua  Variance  Clue  tar  Analysis 


Eigenvalues  of  the  Covariance  Matrix 


1 

Eigenvalue 

0.697871 

Difference 

0.638318 

Proportion 

0.874661 

Cumulative 

0.87465 

2 

0.069663 

0.037618 

0.074639 

0.94929 

3 

0.022036 

0.007166 

0.027617 

0.97691 

4 

0.014879 

0.012880 

0.018648 

0.99655 

6 

0.001999 

0.000452 

0.002506 

0.99806 

6 

0.001647 

0.001647 

C. 001939 

1.00000 

T 

0.000000 

0.000000 

0.000000 

1.00000 

8 

0.000000 

0.000000 

0.000000 

1.00000 

9 

0.000000 

0.000000 

0.000000 

1.00000 

10 

0.000000 

0.000000 

0.000000 

1.00000 

11 

0.000000 

0.000000 

0.000000 

1.00000 

12 

0.000000 

0.000000 

0.000000 

1.00000 

13 

0.002000 

0.000000 

0.000000 

1.00000 

14 

0.000000 

0.000000 

0.000000 

1.00000 

16 

0.000000 

0.000000 

0.000000 

1.00000 

16 

0.000000 

0.000000 

0.000000 

1.00000 

17 

0.000000 

0.000000 

o.oocooc 

1.00000 

18 

0.000000 

0.000000 

0.000000 

1.00000 

19 

0.000000 

0.000000 

0.000000 

1.00000 

20 

o.eocooo 

0.000000 

0.000000 

1.00000 

21 

0.000000 

0.000000 

0.000000 

1.00000 

22 

0.000000 

0.000000 

0.000000 

1.00000 

23 

0.000000 

0.000000 

o.oocooo 

1.00000 

24 

0.000000 

0.000000 

0.000000 

1.00000 

26 

0.000000 

0.090000 

0.000000 

1.00000 

26 

0.000000 

0.000000 

0.000000 

1.00000 

27 

-.000000 

0.000000 

-.000000 

1.00000 

28 

-.000000 

0.000000 

-.000000 

1.00000 

29 

-.000000 

0.000000 

-.000000 

1.00000 

30 

-.000000 

0.000000 

-.090009 

1.00000 

31 

-.000000 

0.000000 

-.000000 

1.00000 

32 

-.000000 

o.oocooo 

-.000000 

1.00900 

33 

-.000000 

0.000000 

-.000000 

1.00000 

34 

-.000000 

0.000000 

-.000000 

1.00000 

38 

-.000000 

0.000000 

-.000000 

1.00000 

36 

-.000000 

0.000000 

-.090000 

1.00000 

37 

-.000000 

0.000000 

-.000000 

1.00000 

38 

-.000000 

0.000000 

-.000900 

1.00000 

39 

-.000000 

0.000000 

-.000000 

1.00000 

40 

-.000000 

0.000000 

-.000000 

1.00000 

41 

-.000000 

0.000000 

-.000000 

1.000C0 

42 

-.0000^ 

0.000000 

-.000000 

1.00000 

43 

-.000000 

0.000090 

-.000000 

1.00000 

44 

-.000000 

0.000000 

-.000000 

1.00000 

46 

-.000000 

0.009000 

-.000000 

1.00000 

48 

-.000000 

0.000000 

-.000000 

1.00000 

47 

-.000000 

• 

-.000000 

1.09000 

ftoot-Rean-Square  Total-Saaple  Standard  Deviation  »  0.130293 
Root -Mean-Square  Distance  Betveea  Observations  "  1,263238 


86 


The  SIS  Systaa  17 

10:S4  Thursday,  Dececber  17,  19S2 


Hard’s  Hinisua  Variance  Cluster  Analysis 


T 

Pseudo  Pseudo  1 


X  Clusters 

Joined 

nzQ 

SPRSQ 

ASQ 

r 

t**3  a 

•  083 

033 

3 

0 . 003050 

0.994941 

46.3 

. 

S  CM 

cn« 

3 

0.004844 

0.990397 

61.4 

3.1 

4  081 

ass 

3 

0.044774 

0.943431 

14.7 

. 

3  CM 

024 

3' 

0.097934 

0.845487 

11.0 

2.1 

3  CM 

CLS 

« 

0.384481 

0.541304 

4.4 

7.4 

1  CL2 

0B7 

7 

0.541304 

0.000000 

. 

4.4 

87 


ixmxxx 

imxm 


zzxzxxxzxxxxx 

xxxnxraxra 

nmnrnm 

nmarnni 

xxmxi 

ixt  nix 


Th«  SAS  Syataa  19 

10:64  Thursday,  D«c«xb«r  17,  1992 


DBS  K?  HT  CP  CT  IP  IT  MP  HIT  KCP  KC?  CIP  CIT  KCIP  HCIT 
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17 
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33 

0 

60 
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36 
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23 
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49 
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14 
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24 
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39 
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31 

0 

13 
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26 
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30 

1 

26 

1 

27 

1 
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60 
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31 
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60 
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26 
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47 
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60 
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39 
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60 
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27 

1 

26 
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60 
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60 
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17 

1 

24 

0 

60 

1 

28 

1 

30 
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60 
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29 
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11m  SAS  Sy*t««  30 

10:64  Tfcnrftdajr,  D«c«mb«r  17,  1992 


1 

Hc«n 

Std  D*t 

Hl.nii.ta 

Kuiaa 

47 

0.8086106 

0.3977271 

0 

1.0000000 

47 

36.4893617 

14.6600292 

8.0000000 

60.0000000 

47 

0.6631916 

0.6026376 

0 

1.0000000 

47 

46.8086106 

16.2681244 

15.0000000 

60.0000000 

47 

0.6106383 

0.5052912 

0 

1.0000000 

47 

47.9787234 

13.9416320 

18.0000000 

60.0000000 

47 

0.8086106 

0.3977271 

0 

1.0000000 

47 

33.4256319 

14.9981806 

10.0000000 

60.0000000 

47 

0.9361702 

0.2470922 

0 

1.0000000 

47 

30.1063830 

12.8574696 

12.0000000 

60.0000000 

47 

0.6744681 

0.4997687 

0 

1.0000000 

47 

46.6631915 

18.0708862 

9.0000000 

60.0000000 

47 

0.8723404 

0.3373181 

0 

1.0000000 

47 

33.7446809 

15.9788829 

9.0000000 

60.0000000 

90 
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083 

RHP 

CRP 

IBP 

imtp 

RCSP 

CIKP 

KCIRP 

1 

33 

37 

40 

35 

31 

69 

42 

2 

40 

40 

40 

18 

38 

40 

47 

3 

27 

19 

37 

40 

40 

40 

57 

4 

40 

34 

40 

40 

30 

27 

40 

5 

40 

43 

58 

40 

30 

33 

34 

4 

43 

38 

33 

24 

16 

40 

32 

T 

35 

33 

40 

30 

.  14 

42 

40 

8 

30 

18 

33 

39 

20 

28 

60 

9 

36 

65 

50 

17 

40 

50 

67 

10 

30 

38 

28 

31 

33 

40 

28 

11 

38 

50 

47 

22 

27 

42 

27 

12 

24 

49 

31 

34 

28 

20 

23 

13 

25 

39 

60 

30 

57 

30 

28 

14 

30 

67 

40 

36 

31 

60 

24 

16 

30 

39 

60 

33  | 

18 

64 

30 

14 

28 

38 

36 

31  ! 

17 

23 

14 

17 

24 

36 

40 

46 

66 

39 

30 

18 

27 

20 

40 

62 

20 

26 

18 

19 

47 

40 

40 

IT  j 

40 

64 

24 

20 

25 

26 

40 

29 

16 

30 

18 

21 

30 

28 

40 

11 

29 

12 

39 

22 

33 

50 

37 

34 

16 

17 

60 

23 

37 

20 

38 

44  ! 

26 

60 

30 

24 

21 

41 

20 

17  ; 

31 

23 

24 

25 

23 

37 

a 

32 

39 

40 

20 

24 

26 

16 

a 

17 

24 

19 

47 

27 

22 

a 

a 

17 

12 

28 

19 

28 

26 

a 

a 

24 

30 

a 

27 

29 

8 

a 

a 

10  | 

24 

a 

22 

30 

27 

a 

a 

41 

16 

a 

18 

31 

19 

a 

a 

19 

29 

a 

20 

32 

42 

a 

35  ' 

28 

a 

16 

33 

42 

61 

25 

a 

19 

«4 

34 

34 

26 

a 

28 

36 

28 

28 

27 

a 

20 

34 

25 

a 

33 

47 

a 

29 

37 

27 

a 

a 

12 

26 

a 

21 

38 

21 

a 

a 

27 

19 

a 

44 

39 

a 

a 

a 

29 

a 

24 

40 

a 

a 

a 

14 

a 

63 

41 

a 

a 

a 

a 

30 

a 

9 

42 

a 

a 

a 

a 

44 

a 

43 

a 

a 

a 

38 

a 

44 

a 

a 

24 

a 

Tbs  >18  Sjatw  32 

10:64  Thursday,  Dscsabsr  17,  1992 


Tariabls  1 

■saa 

Std  Dot 

Rialma 

Karim 

HP 

» 

31.0799474 

19.2808707 

8.0000000 

60.0000000 

cu> 

26 

37.9923077 

13.1902042 

16.0000000 

60.0000000 

nr 

24 

44. >>416907 

13.3773966 

20.0000000 

60.0000000 

nv 

38 

30.8947348 

12.6401109 

10.0000000 

60.1009000 

RCtf 

44 

29.0000000 

11.4624326 

12.0000000 

60.0000000 

cup 

27 

40.2942943 

17.1672260 

12.0000000 

40.0000000 

■CUP 

41 

31.6686366 

14.2909233 

9.0000000 

60.0000000 

Th*  US  Systea  23 

10:54  Tfeorsday,  D»c«b*r  17,  1992 


DBS 

HVP 

CVP 

I¥P 

xm* 

new? 

CIV? 

HCIBP 

1 

60 

60 

67 

60 

60 

63 

60 

2 

60 

60 

60 

60 

19 

60 

60 

3 

60 

60 

62 

60 

60 

60 

60 

1 

60 

60 

60 

26 

22 

39 

6 

60 

60 

60 

60 

60 

60 

• 

60 

60 

31 

14 

• 

60 

9 

7 

60 

60 

18 

19 

• 

60 

• 

• 

60 

60 

60 

60 

• 

60 

9 

64 

60 

30 

38 

• 

60 

• 

10 

60 

60 

. 

• 

60 

• 

11 

60 

60 

• 

60 

12 

60 

19 

• 

60 

13 

60 

60 

• 

60 

14 

60 

39 

• 

• 

60 

• 

IS 

60 

46 

60 

• 

16 

60 

68 

• 

60 

• 

17 

• 

60 

60 

• 

60 

• 

18 

• 

20 

60 

• 

• 

9 

• 

19 

• 

60 

60 

• 

• 

60 

• 

20 

• 

60 

60 

• 

• 

9 

• 

21 

• 

60 

37 

• 

• 

• 

• 

22 

• 

• 

60 

• 

• 

• 

• 

23 

• 

• 

31 

• 

• 

• 

• 

93 


Th*  8AS  8yata»  24 
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Varlabla 

I 

Ha  an 

Std  Dav 

Minina 

Mariana 

IMP 

9 

69.3333333 

2.0000000 

64.0000000 

60.0000000 

CVP 

21 

68.0962381 

8.7287166 

20.0000000 

60.0000000 

iw 

23 

49.4782609 

14.6531708 

18.0000000 

60.0000000 

nv 

9 

44.1111111 

19.8899761 

14.0000000 

60.0000000 

RCRP 

3 

46.3333333 

23.6713610 

19.0000000 

60.0000000 

CHIP 

20 

62.6600000 

17.1901900 

9.0000000 

60.0000000 

HCT9P 

« 

48.0000000 

20.6710326 

9.000C300 

60.0000000 

Mjj— ,1  .  u  J  ■■  Kuaafi  Hgp 


Tha  SIS  ftystaa  26 
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OBS  HIA  JOT  CIA  CST  HA  1ST  HUA  HI  IT  KCIA  HC8T  CIIA  CIBT  HCIIA  HCIffT 


1 

1 

33 

1 

37 

0 

67 

1 

36 

1 

31 

1 

69 

1 

42 

2 

1 

27 

1 

19 

1 

40 

1 

18 

1 

38 

0 

63 

1 

47 

3 

1 

40 

1 

34 

1 

37 

1 

40 

1 

40 

0 

22 

1 

67 

4 

1 

43 

1 

43 

0 

52 

1 

26 

1 

30 

1 

27 

1 

34 

6 

1 

36 

1 

38 

1 

68 

1 

30 

1 

30 

1 

33 

1 

32 

6 

1 

30 

1 

33 

0 

31 

4 

39 

1 

16 

1 

42 

1 

60 

7 

1 

35 

1 

18 

1 

33 

1 

17 

1 

14 

1 

28 

1 

67 

8 

1 

30 

1 

66 

0 

18 

J 

31 

1 

20 

1 

42 

1 

28 

9 

1 

38 

1 

38 

1 

33 

1 

22 

1 

33 

1 

20 

1 

27 

10 

1 

24 

1 

60 

0 

30 

4 

34 

1 

27 

1 

30 

1 

23 

11 

1 

26 

1 

49 

1 

28 

0 

26 

1 

28 

1 

64 

0 

39 

12 

1 

30 

J 

39 

1 

47 

1 

30 

1 

67 

1 

23 

1 

28 

13 

1 

30 

1 

57 

0 

19 

1 

36 

1 

31 

1 

39 

1 

26 

14 

1 

28 

1 

39 

1 

31 

1 

33 

1 

18 

1 

26 

1 

30 

15 

1 

24 

1 

38 

0 

39 

1 

31 

1 

17 

1 

66 

1 

14 

16 

1 

27 

1 

36 

1 

36 

1 

46 

1 

56 

1 

SO 

1 

30 

17 

1 

47 

1 

20 

0 

46 

0 

14 

1 

20 

1 

12 

1 

18 

18 

1 

25 

1 

25 

1 

40 

0 

19 

1 

40 

1 

17 

1 

26 

19 

1 

30 

1 

28 

0 

68 

1 

62 

1 

16 

0 

9 

1 

18 

20 

1 

33 

1 

50 

1 

40 

1 

17 

1 

29 

0 

9 

1 

39 

21 

1 

37 

0 

20 

1 

37 

1 

29 

1 

15 

1 

23 

1 

30 

22 

1 

21 

1 

20 

0 

37 

1 

11 

1 

26 

1 

19 

1 

24 

23 

1 

23 

1 

41 

1 

38 

1 

34 

1 

31 

1 

28 

1 

20 

24 

1 

26 

1 

37 

0 

31 

1 

44 

1 

39 

a 

a 

1 

47 

26 

1 

22 

1 

15 

1 

20 

1 

17 

1 

24 

a 

a 

1 

19 

26 

1 

26 

a 

a 

a 

a 

1 

32 

1 

12 

a 

a 

1 

27 

27 

1 

8 

a 

a 

a 

a 

0 

38 

0 

19 

a 

a 

1 

22 

28 

1 

27 

a 

a 

a 

a 

1 

17 

1 

30 

a 

a 

1 

18 

29 

1 

19 

a 

a 

a 

a 

1 

17 

1 

24 

a 

1 

20 

30 

1 

42 

a 

a 

a 

a 

1 

26 

1 

16 

1 

16 

31 

1 

42 

a 

a 

1 

10 

1 

29 

1 

19 

32 

1 

34 

a 

a 

1 

41 

1 

28 

1 

28 

33 

1 

28 

a 

a 

a 

1 

19 

1 

25 

1 

20 

34 

0 

64 

a 

a 

1 

36 

1 

26 

1 

29 

36 

1 

26 

a 

1 

61 

1 

27 

1 

21 

36 

1 

27 

1 

34 

1 

47 

0 

9 

37 

1 

21 

1 

28 

1 

26 

1 

44 

38 

a 

1 

33 

1 

19 

1 

26 

39 

a 

a 

1 

12 

1 

29 

1 

63 

40 

a 

a 

1 

27 

1 

14 

1 

9 

41 

a 

a 

a 

a 

1 

30 

a 

a 

a 

42 

a 

a 

a 

a 

1 

44 

a 

a 

a 

43 

a 

a 

• 

, 

a 

1 

38 

a 

a 

a 

44 

• 

a 

• 

• 

• 

1 

26 

a 

95 
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Terlabl* 

i 

Kaaa 

Std  Dav 

Minima 

Maxima 

nu 

37 

0.972?730 

0.1643990 

0 

1.0000000 

HST 

37 

30.1351361 

8.7087509 

8.0000000 

54.0000000 

CIA 

25 

0.9600000 

0.2000000 

0 

1.0000000 

orr 

25 

35.2000000 

12.0968315 

15.0000000 

57.0000000 

m 

26 

0.6600000 

0.5066228 

0 

l.OOCOOOO 

ht 

25 

37.4000000 

11.1130654 

18.0000000 

68.COOOOOO 

HUA 

40 

0.9000000 

0.3038218 

0 

1.0000000 

HUT 

40 

28.7760000 

10.7881404 

10.0000000 

62.0000000 

HCIA 

44 

0.9772727 

0.1607557 

0 

1.0000000 

KCIT 

44 

28.0691813 

10.5000503 

12.0000000 

57.0000000 

cm 

23 

0.8260870 

0.3875534 

0 

1.0000000 

cut 

23 

30.4782609 

14.7490202 

9.0000000 

59.0000000 

HCIIA 

40 

0.9500000 

0.2207214 

0 

1.0000000 

hcht 

40 

29.1500000 

12.4808571 

9.0000C00 

57.0000000 

iSaassi 


aa  ~  ill'  ■ . '  J  .  -  ■•  -  ;  ■  -  -  -  ■--- 
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OBS 

?TOi 

COA 

IOA 

MICA 

RCOA 

CIOA 

NCIOA 

1 

1 

0 

c 

0 

0 

1 

1 

2 

0 

0 

1 

1 

1 

0 

0 

3 

0 

1 

1 

1 

0 

1 

1 

4 

1 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

« 

0 

0 

1 

0 

0 

1 

7 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

9 

0 

0 

1 

0 

io 

0 

0 

0 

1 

li 

0 

0 

0 

12 

0 

1 

0 

13 

0 

1 

1 

14 

a 

0 

1 

1 

15 

• 

0 

1 

0 

16 

• 

0 

0 

1 

17 

0 

1 

0 

18 

1 

1 

0 

19 

0 

0 

• 

0 

a 

20 

0 

0 

a 

0 

a 

21 

• 

0 

0 

• 

• 

0 

a 

22 

• 

0 

0 

• 

• 

1 

« 

23 

• 

• 

. 

• 

• 

0 

• 

24 

• 

. 

a 

• 

1 

a 

97 
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Variable  I  feu  Std  Dev  Rlalau  Harlan 


■u 

10 

0.2000000 

0.4219370 

1.0000000 

CM 

32 

0.0000001 

0.2042440 

1.0000000 

ztu 

22 

0.4545456 

0.5009472 

l.*WOOOO 

■nu 

r 

0.2057143 

0.4979600 

1.0000000 

mu 

0 

0.3333333 

0.5773603 

1.0000000 

cnu 

24 

0.3333333 

0.4915434 

1.0000000 

■cm 

7 

0.4295714 

0.6345225 

1.0000000 

The  SAS  System  29 

10; $4  Thursday,  Deceaber  17,  1992 
Coaplete  Linkage  Cluster  Analysis 


Eigenvalues  of  the  Covariance  Matrix 


Eigenvalue 

Difference 

Proper  don 

Cumulative 

1 

4024.21 

2288.74 

0.783063 

0.38305 

2 

1735.47 

47.18 

0.165196 

0.64826 

S 

1688.29 

297.21 

0.160704 

0.70896 

4 

1391.08 

360.47 

0.132413 

0.84137 

6 

1030.62 

394.68 

0.098101 

0.93947 

9 

636.94 

636.94 

0.060633 

1.00000 

7 

0.00 

0.00 

0.000000 

1.00000 

« 

0.00 

0.00 

0.000000 

1.00000 

9 

0.00 

0.00 

0.000000 

1.00000 

10 

0.00 

0.00 

0.000000 

1.C0000 

11 

0.00 

000 

0.000000 

1.00000 

12 

0.00 

0.00 

0.000000 

a  00000 

13 

0.00 

0.00 

0.000000 

1.90000 

14 

0.00 

0.00 

0.000000 

1.00000 

16 

0.00 

0.00 

0.000000 

1.00000 

19 

0.00 

0.00 

0.000000 

1.00000 

17 

0.00 

0.00 

0.000000 

1.00000 

18 

0.00 

0.00 

0.000000 

1.00000 

19 

0.00 

0.00 

0.000000 

1.00000 

20 

0.00 

0.00 

0.000000 

1.00000 

21 

0.00 

0.00 

0.000000 

1.00090 

22 

0.00 

0.00 

0.000000 

l.OOOCO 

23 

0.00 

0.00 

0.000000 

1.00000 

24 

0.00 

0.00 

0.000000 

1.00000 

26 

0.00 

0.00 

0.000900 

1.00000 

26 

0.00 

0.00 

0.000000 

1.00000 

27 

0.00 

0.00 

0.000000 

1.00000 

28 

0.00 

0.00 

0.000000 

1.00000 

29 

0.00 

0.00 

0.000000 

1.00000 

30 

-0.00 

0.00 

-.000000 

1.00000 

31 

-0.00 

0.00 

-.000000 

1.00000 

32 

-0.00 

0.00 

-.000000 

1.00000 

33 

-0.00 

0.00 

-.000000 

1.00000 

34 

-0.00 

0.00 

-.000000 

1.00000 

36 

-0.00 

0.00 

-.000000 

1.00000 

36 

-0.00 

0.00 

-.000000 

1.00000 

37 

-0.00 

0.00 

-.000000 

1.00000 

38 

-0.00 

0.00 

-.000000 

1.00000 

39 

-0.00 

0.00 

-.000000 

1.00000 

40 

-0.00 

0.00 

-.000000 

1.00000 

41 

-0.00 

0.00 

-.000000 

1.00000 

42 

-0.00 

0.00 

-.000000 

1.00000 

43 

-0.00 

0.00 

-.000000 

1.00000 

44 

-0.00 

o.oc 

-.000000 

1.00000 

46 

-0.00 

0.00 

-.000000 

1.00000 

49 

-0.00 

0.00 

-.000000 

1.00000 

47 

-0.00 

-.000000 

1.00000 

Boot-Maan-Square  Total-Saaple  Standard  Deviation  ■  14.96071 
■•an  Distance  Between  Obsecrations  •  143.9466 
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Complete  Linkage  Cluster  Analysis 


Busbar 

Frequency 

lonalised 

Si 

of  lev 

Has  Ira 

Clusters 

Clusters  Joined 

Cluster 

Distance 

6 

OBI 

DB7 

2 

0.749677 

CM 

006 

2 

0.801696 

on 

006 

2 

0.887343 

3 

CL6 

05 

4 

0.946333 

CL4 

003 

3 

0.974192 

1 

CL3 

02 

7 

1.224618 

0 


i 


f 


Tha  SIS  fjyataa  32 

10:64  Thursday,  Dtcnbtr  17,  1292 
Avaragn  Linkaga  Cl Qatar  inalyaia 


Eiganraluai  of  tha  CoTarlaaea  Matrix 


Eifanralna 

Diffaraoca 

Proportion 

Coxal at iva 

1 

4024.21 

2288.74 

0.383053 

0.38305 

2 

1736.47 

47.18 

0.185196 

0.64826 

3 

1888.29 

297.21 

0.180704 

0.70895 

4 

1391.08 

380.47 

0.132413 

0.84137 

1030.82 

384.88 

0.098101 

0.93947 

« 

836. 84 

836.94 

0.080533 

1.00000 

7 

0.00 

0.00 

0.000000 

1.00000 

8 

0.00 

0.00 

0.000000 

1.00000 

9 

0.00 

0.00 

0.009000 

1.00000 

10 

0.00 

0.00 

0.000000 

1.00000 

11 

0.00 

0.00 

0.000000 

1.00000 

12 

0.00 

0.00 

6.000000 

1.00000 

0.00 

0.00 

0.000000 

1.00000 

14 

0.00 

0.00 

0.000000 

1.00000 

16 

0.00 

0.00 

0.000000 

1.00000 

18 

0.00 

0.00 

0.000000 

1.00000 

17 

0.00 

0.00 

0.000000 

1.00000 

18 

0.00 

0.00 

0.000000 

1.00000 

19 

0.00 

0.00 

0.000000 

1.00000 

20 

0.00 

0.00 

0.000000 

1.00000 

21 

0.00 

9.00 

0.000000 

1.00000 

22 

0.00 

0.00 

'  0.000000 

1.00000 

23 

0.00 

0.00 

0.000000 

1.00000 

24 

0.00 

0.00 

0.000000 

1.00000 

28 

0.00 

0.00 

9.000000 

1.00000 

28 

0.00 

0.00 

0.000000 

1.00000 

27 

0.00 

0.00 

0.000000 

1.00000 

28 

0.00 

0.00 

0.000009 

1.00000 

29 

0.00 

0.00 

0.000000 

1.00000 

30 

-0.00 

0.00 

-.000000 

1.00000 

31 

-0.00 

0.00 

-.000000 

1.00000 

32 

-0.00 

0.00 

-.000000 

1.00000 

33 

-0.00 

0.00 

-.009000 

1.00000 

34 

-0.00 

0.00 

-.000000 

1.00000 

38 

-0.00 

0.00 

-.000000 

1.00000 

98 

-0.00 

0.00 

-.000000 

1.00000 

37 

-0.00 

0.00 

-.000000 

1.00000 

38 

-0.00 

0.00 

-.000000 

1.00000 

39 

•0.00 

0.00 

-.000000 

1.00000 

40 

-0.00 

0.90 

-.000000 

1.00009 

41 

-0.00 

0.00 

-.000000 

1.00000 

42 

-0.00 

0.00 

-.000000 

1.00000 

43 

•0.00 

0.00 

-.000000 

1.00000 

44 

•0.00 

0.00 

-.000000 

1.00000 

46 

-0.00 

0.00 

-.000000 

1.00000 

48 

-0.00 

0.00 

-.000000 

1.00000 

47 

-0.00 

a 

-.000000 

1.00000 
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ATerage  Linkage  Cloatar  Analysis 


lunber 

Frequency 

■oxaalixod 

of 

of  8a> 

Pseudo 

Pseudo 

8)0 

Clusters 

Clusters  Joined 

Cluster 

F 

t**2 

Distance 

6 

081 

QB7 

2 

1.87 

0.743121 

6 

084 

085 

2 

2.04 

. 

0.794472 

4 

06 

05 

4 

1.85 

1.66 

0.870167 

3 

QB2 

QB6 

2 

2.17 

• 

0.879348 

2 

03 

083 

3 

2.80 

1.26 

0.957954 
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7 

a 

2.80 

1.087015 
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Ward’s  Minima  Varianca  Clustar  Analysis 


Eiganvaluas  of  tha  Corariancs  Ratrix 


Eiganvalna 

Diffaranca 

Proportion 

Cuaulativa 

1 

4024.21 

2288.74 

0.383053 

0.38305 

2 

1736.47 

47.18 

0.165195 

0.54825 

3 

1688.29 

297.21 

0.160704 

0.70895 

4 

1391.08 

360.47 

0.132413 

0.84137 

5 

1030.62 

394.68 

0.098101 

0.93947 

5 

635.94 

635.94 

0.060533 

1.00000 

7 

0.00 

0.00 

0.000000 

1.00000 

• 

0.00 

0.00 

0.000000 

1.00000 

9 

0.00 

0.00 

0.000000 

1.00000 

10 

0.00 

0.00 

0.000000 

1.00000 

11 

0.00 

0.00 

0.000000 

1.00000 

12 

0.00 

0.00 

0.000000 

1.00000 

13 

0.00 

0.00 

0.000000 

1.00000 

14 

0.00 

0.00 

0.000000 

1.00000 

16 

0.00 

0.00 

0.000000 

1.00000 

16 

0.00 

0,00 

0.000000 

1.00000 

17 

0.00 

0.00 

0.000000 

1.00000 

16 

0.00 

0.00 

0.000000 

1.00000 

19 

0.00 

0.00 

0.000000 

1.00000 

20 

0.00 
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0.000000 

1.00000 
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1.00000 
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0.00 
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1.00000 

23 

0.00 

0.00 

0.000000 

1.00000 

24 

0.00 

0.00 

0.000000 

1.00000 

25 

0.00 

0.00 

0.000000 

1.00000 

26 

0.00 

0.00 

0.000000 

1.00000 

27 

0.00 

0.00 

0.000000 

1.00000 

28 

0.00 

0.00 

0.000000 

1.00000 

29 

0.00 

0.00 

0.000000 

1.00000 

30 

-0.00 

0.00 

-.000000 

1.00000 

31 

-0.00 

0.00 

-.000000 

1.00000 

32 

-0.00 

0.00 

-.000000 

1.00000 

33 

-0.00 

0.00 

-.000000 

1.00000 

34 

-0.00 

0.00 

-.000000 

1.00000 

35 

-0.00 

0.00 

-.000000 

1.00000 

36 

-0.00 

0.00 

-.000000 

1.00000 

37 

-0.00 

0.00 

-.000000 

1.00000 

38 

-0.00 

0.00 

-.000000 

1.00000 

39 

-0.00 

0.00 

-.000000 

1.00000 

40 

-0.00 

0.00 

-.000000 

1.00000 

41 

-0.00 

0.00 

-.000000 

1.00000 

42 

-0.00 

0.00 

-.000000 

1.00000 

43 

-0.00 

0.00 

-.000000 

1.00000 

44 

-0.00 

0.00 

-.000000 

1.00000 

45 

-0.00 

0.00 

-.000000 

1.00000 

46 

-0.00 

0.00 

-.000000 

1.00000 

47 

-0.00 

a 

-.000000 

1.00000 

Baot-flaan-Sqaara  Total-Saapls  Standard  Dariatica  «  14.95071 
Koot-Maan-Squara  Distanca  Bataaaa  Obaarrations  •  144.9525 
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Sard’s  Minina  Variants  Clnstar  Analysis 


Psando  Psando  i 


■CL  Clast  ara  Joinad 

fizq 

sptsq 

KSQ 

F 

ta*3 

<  081 

car 

3 

0.093038 

0.907963 

3.0 

6  084 
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3 
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3.0 
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Varlabla  1 

Haas 

Std  Dav 

Mlaiaua 

Karina 

UOP 

47 

0.6608611 

0.4711864 

0 

1.0000000 

Loonr 

47 

43.2137660 

17.4915336 

14.0000000 

60.0000000 

SLTXP 

47 

0.7659674 

0.4279763 

0 

1.0000000 

SLOT 

47 

41.2663191 

14.2806260 

19.0000000 

60.0000000 

5IXDP 

47 

0.8723404 

0.3373181 

0 

1.0000000 

SLUT 

47 

35.8936170 

12.1032407 

17.0000000 

60.0000000 

717 

47 

0.8086106 

0.3977271 

0 

1.0000000 

ITT 

47 

36.2766967 

15.7793235 

8.0000000 

60.0000000 

FU> 

47 

0.4680861 

0.5043749 

0 

1.0000000 

FLT 

47 

46.3617021 

15.9692987 

9.0000000 

60.0000000 

8U> 

47 

0.7872340 

0.4136881 

0 

1.0000000 

SLT 

47 

36.6695746 

16.7114600 

9.0000000 

60.0000000 

8TP 

47 

0.7021277 

0.4622673 

0 

1.0000000 

STT 

47 

34.9574468 

20.1461588 

12.0000000 

60.0000000 

Appendix  C.  SLAM  and  Fortran  Code 

SLAM  Code 

This  section  contains  the  SLAM  network  code  used.  The  first  model  shown  is  the 
original  model  (Pritsker,  1986:240).  Figure  33  shows  the  SLAM  network  diagram.  The 
next  seven  models  are  the  network  portion  of  the  modified  model  for  each  scenario.  The 
associated  Fortran  code  is  shown  in  the  next  section.  Figuren  page  113  shows  the  network 
diagram  for  the  LDER  scenario  model.  The  scenarios  were: 

•  LDER  -  Load  interarrival  time,  loading  time,  hauling  time,  dumping  time,  and  return 
time  unchanged.  (Original  Model  with  Modifications) 

•  SLTK  -  Hauling  time,  dumping  time,  and  return  time  doubled  for  the  third  truck. 
(Slow  Truck) 

•  SLLD  -  Loading  time  doubled  for  the  first  loader.  (Slow  Loader) 

•  FT  •  Hauling  time,  dumping  time,  and  return  time  cut  in  half  for  all  trucks.  (Fast 
Trucks) 

•  FL  -  Loading  times  halved  for  both  loaders.  (Fast  Loaders) 

•  SL  •  Load  interarrival  time  multiplied  by  2.  (Slow  Loads) 

•  ST  -  Hauling  time,  dumping  time,  and  return  time  doubled  for  all  trucks.  (Slow 
Trucks) 


Figure  33.  SLAM  Network  Diagram  for  Original  Model. 


Original  model. 


GEM, PRITSKER, TRUCK  HAULIMG.7/12/83,1, , ,Y/H; 
LIMITS, 3, 1,50; 

HETVORK; 

START  CREATE; 

ACT, , THOU. LE. 480; 

GOOM; 

ACT  .ERLHG  (4 . ,  2 . ) , ,  START ; 

ACT; 

ACCUM,2,2; 

LOAD  QOEUE(l) , , , ,BEGH; 

TRKS  QUEUE(2) ,4, , ,BEGN; 

LDER  QUEUE (3) ,2, , ,BEGN; 

BEGI  SELECT, ASM, LIT,, LOAD, TRKS. LDER; 

ACT/1 ,EXP0H(14.); 

ACT/2, EXPOS (1 2. ) ; 

GOOM; 

ACT, 5,, LDER; 


ACT,RI0RM(22.,3.); 


GOOM; 


CREATE  LOAD  TRAMS  ACTIONS 
STOP  IF  AFTER  8  HOURS 
ELSE 

BRAXCH  BACK  TO  START 

AND  CONTINUE 

ACCUMULATE  TVO  PILES 

QUEUE  OF  LOADS 

QUEUE  OF  TRUCKS 

QUEUE  0?  LOADERS 

ASM  OF  LOAD, TRKS.  AND  LDER 

L0ADER1  TIKE 

L0ADER2  TIME 

LOADER  RESTIN'*  TIKE 
TRUCK  HAULING  TIME 
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Figure  34.  SLAM 
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Modified  mode!  used  for  LDER  scenario.  All  times  were  unchanged  from  the  original 
model. 

QEI , CARPESTER . TRUCK  HAULXIG , 9/30/92 . 1 . . . Y/> . , . 72 ; 


LIMITS, 5, 14, 50; 

BITRY/2,0,0.0,0,0,0,0,0,0,1,0; 

EHttY/2, 0,0, 0,0, 0,0, 0,0, 0,2,0; 

EVTRT/2, 0,0, 0,0, 0,0, 0,0 ,0,3,0; 

DTRY/2,0,0, 0,0,0, 0,0, 0,0, 4,0; 

HTRY/3,0,0,0,0,0,0,0,0, 0,0,5; 

EITRY/3.0,0,0,0,0,0,0,0,0,0,6; 

VETV0RK; 

RES0URCE/L0 ADSP0T ( 1 ) ,4; 

RESOURCE/DUMPSPOT ( 1 ) ,6; 

START  CREATE; 

ASSIGI.ATRZB(l)  ■  IZ, 

ATRZBC2)  -  ERLVG(4. ,2. ) , 
ATRZB(4)  ■  EZP0K14.), 
ATRXB(S)  -  EXP0M(12.)f 
ATRZB(6)  >  5. 

ATRZBC7)  -  RV0RK(22. ,3.) , 
ATRXB(8)  -  UIFRM(2.,8.), 
ATBZB(9)  •  RH0RM(18. ,3.) ; 

ACT/1,. TIOW.LE. 480; 

0001; 

ACT/2.A?RZB(2),, START; 

ACT/3; 

ACCOM, 2, 2; 

LOAD  QUEUE(1),,,,BEGI; 

TRKS  QUEUE(2) , , , ,BEGI; 

LDER  QUEUE(3) , , , ,BEGI ; 

BEGI  SELECT, ASK/SUH,,, LOAD.TRKS, LDER; 

ACT/13; 

AVAIT(4) .LOADSPOT/l ,,1; 

ACT/4, ATRZB (4) .ATRZB(ll) .EQ.6.FRID; 

ACT/6, ATRZB (5) .ATRZB(ll) .EQ.8; 

FRID  FRB8.L0ADSP0TC1) ; 

0001; 

ACT/12,.. HAUL; 

ACT/8, ATRZB (6); 

ASSZGI,ATRZB(1)>0, 

ATRIB(2)-0, 

ATRXB(3)-0, 

ATRZB(4)*0, 

ATRIB(5)-0, 

ATRZB (6) «0, 


CREATE  LOAD  TRAISACTIOIS 


STOP  ZF  AFTER  8  HOURS 

WT-gg 

BRAICB  BACK  TO  START 
AID  COITXIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD.TRKS,  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTZIG  TIME 


I 

■A 

j 


i 

t 


/ 


ATRIB<7)«0, 
ATRIB (8) "0, 
ATRIB (9) »0, 
ATRIB C 10) -0; 

ACT/11,,,LDER; 

HAUL  GOOH; 

ACT/7, ATRIB (7); 

AHAIT(S) .DUHPSPOT/l; 

ACT/8, ATRIB (8) ; 

FREE.DUMPSPOTd) ; 

ACT/9, ATRIB (9); 

ASSIGI, ATRIB (l)-0, 
ATRIB (2) »0, 
ATRIB(3)*0, 
ATRIB(4)"0, 
ATRIB (5) »0, 
ATRIB(6)«0, 
ATRIB(7)-0, 
ATRIB (8) "0, 
ATRIB (9) -0. 
ATRIB(ll)-0; 

ACT/10, , ,TRKS; 

BIO; 

MOITR, TRACE, 0,480; 

FI*; 


TRUCK  HAULIIG  TIKE 
TRUCK  DUHPIIG  TIKE 
TRUCK  RETURI  TIKE 
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Modified  model  used  for  SLTK  scenario.  Attributes  7,  8,  and  9  (hauling,  dumping,  and 
return  time)  were  multiplied  by  2  for  the  third  truck.  This  was  accomplished  by  branching 
on  Attribute  10  (truck  number)  after  the  select  node. 

GER.CARPERTER, TRUCK  HAULIRG, 9/30/92,1, , ,Y/R, , .72; 


LIMITS, 8,14,60; 

ERTRY/2, 0,0, 0,0, 0,0,0, 0,0, 1,0; 

BITRT/2, 0,0, 0,0, 0,0, 0,0, 0,2,0; 
nrntY/2,0,0.0,0,0,0.0,0,0,3,0; 

ERTRY/2, 0,0, 0,0, 0,0, 0,0, 0,4,0; 
ERTRY/3,0,0,0,0,0,0,0,0,0,0,6; 
EITRY/3,0,0,0,0,0,0,0,0,0,0,6; 

IETU0RK; 

RESOURCE/LOADSPOT ( 1 ) ,4; 
RESOTJRCE/DUMPSPOT  ( 1)  .6; 

START  CREATE; 

ASSIGV.ATRIB(l)  -  II, 

ATRIBC2)  •  ERLHG(4. ,2.) , 
ATRIB(4)  >  EXP0NU4.), 
ATRIB(S)  •  EXP0M(12.) , 
ATRIBC6)  •  6, 

ATRIB(7)  -  RV0RM(22. ,3.) , 
ATRIB(B)  •  URFRM(2. ,6.) , 
ATRIB(9)  -  RR0RM(18. ,3.) ; 
ACT/1, .TIOW.LE. 480; 

G001; 

ACT/2,ATRIB(2),, START; 

ACT/3; 

ACCOM, 2, 2; 

LOAD  QUEUE(l) , , , ,BEGH; 

TRKS  QUEUE (2) , , , ,BEGR ; 

LDER  QUEUE (3) , , , ,BEGR ; 

BEGR  SELECT, ASK/SUM,,, LOAD, TRKS.LDER; 
ACT/13; 

GOOR; 

ACT, ,ATRIB(10) .RE.3.LSRD; 

ACT, ,A7RIB(10) .EQ.3; 

ASSIGR,ATRIB(7)  •  ATRIBC7)  *  2., 
ATRIBC8)  »  ATRIBC8)  *  2., 
ATRIBC9)  •  ATRIBC9)  •  2.; 
LSRD  AWAIT (4) .LOAOSPOT/l , , 1 ; 

ACT/4,  ATRIB  (4)  ,ATRIB(U)  .EQ.6.FRRD; 
ACT/6, ATRIB(S) .ATRIB(ll) .EQ.6; 

FRRD  FREE , LOADSPOT ( 1 ) ; 

GOOR; 

ACT/12,,, HAUL; 
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CREATE  LOAD  TRARSACTIORS 


STOP  IF  AFTER  8  HOURS 
ELSE 

BRARCB  BACK  TO  START 
ARD  CORTIRUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD. TRKS,  ARD  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


l 


ACT/6, ATRIB (6); 

ASSIGN, ATRIB(l)«0, 
ATRIB(2)-0, 
ATRIB(3)-0, 
ATRIB (4) -0, 
ATRIB(5)«0, 
ATRIB(6)*0, 
ATRIB (7) «0, 
ATRIB (8) »0, 
ATRIB(9)*»0, 
ATRIB ( 10) -0; 

ACT/11, , ,LDER; 

HAUL  COOI; 

ACT/7, ATRIB (7); 

AWAIT (5) .DUHPSPOT/l ; 

ACT/8, ATRIB (8); 

PREE,DUMPSPOT(l) ; 

ACT/9,  ATRIBO); 

ASSIGN, ATRIB (1)*0, 
ATRIB(2)«0, 
ATRIB (3) "0, 
ATRIB (4) “0, 
ATRIB (5) «0, 
ATRIB (6) *0, 
ATRIB (7) «0, 
ATRIB (8) *0, 
ATRIB(9)»0, 
ATRIB(11)*0; 

ACT/10,,, TR1CS; 

EID; 

MOMTR, TRACE, 0,480; 

PH; 
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LOADER  RESTING  TIKE 


TRUCK  HAULING  TIME 
TRUCK  DUMPING  TlhE 
TRUCK  RETURN  TIKE 


Modified  model  used  for  SLLD  scenario.  Attribute  4  (loading  time  for  first  loader)  was 
multiplied  by  2. 

GEI,CARPEITER, TRUCK  HAULIIG, 9/30/92, 1, , ,Y/I, , ,72; 

LIMITS, 6,14,60; 

EITRY/2, 0,0, 0,0, 0,0,0, 0,0, 1,0; 

EITRY/2, 0,0, 0,0, 0,0, 0,0, 0,2,0; 

EITRY/2, 0,0, 0,0, 0,0,0, 0,0, 3,0; 

EITRY/2,0,0.0,0,0,0,0,0,0,4,0; 

EMTRY/3, 0,0, 0,0, 0,0, 0,0, 0,0,5; 

E1TRY/3, 0,0, 0,0, 0,0,0, 0,0,0,6; 

IETV0RK; 

RESOURCE/LOAD SPOT ( 1 ) ,4; 

RESOURCE/DUMPSPOTd)  ,5; 

START  CREATE;  CREATE  LOAD  TRAISACTIOIS 

ASSICI, ATRIB(l)  >  II, 

ATRIB(2)  ■  ERLXG(4.,2.),  I 

ATRIBC4)  ■  EXP0IC14.)  *  2.. 

ATRIBC5)  •  EXP0K12.), 

ATRIBC6)  -  5, 

ATRIBC7)  -  RICRM(22. ,3.) , 

ATRIBC8)  ■  UHFRM(2. ,8.) ,  j 

ATRIBO)  •  R10RH(18.  ,3.) ; 

ACT/1,. Ti0W.LE.480;  STOP  IF  AFTER  8  HOURS 

COOK;  ELSE 

ACT/2, ATRIBC2),, START;  BRAICH  BACK  TO  START 

ACT/3;  AID  COITIIUE 

ACCOM.2,2;  ACCUMULATE  TWO  PILES 

LOAD  QUEUE(l) , , , ,BEGI;  QUEUE  OF  LOADS 

TRKS  QUEUE (2) , , , ,BEGI;  QUEUE  OF  TRUCKS 

LDER  QUEUE (3) , , , ,BEGI ;  QUEUE  OF  LOADERS 

BEGS  SELECT, ASM/SUM,,, LOAD, TRKS, LDER;  ASM  OF  LOAD, TRKS,  AID  LDER 

ACT/13; 

AWAIT(4) , LOADSPOT/1 , , 1 ; 

ACT/4, ATRIB (4) .ATRIB(ll) .EQ.S.FRID;  L0ADER1  TIME 

ACT/S, ATRIB(5) ,ATRIB(11) .EQ.6;  L0ADER2  TIME 

FRID  FREE.LOADSPOT(l); 

GOOV; 

ACT/12,,, HAUL; 

ACT/6, ATRIB(6);  LOADER  RESTIVG  TIME 

ASSIGI, ATRIB (1)*0, 

ATRIB (2) *0, 

ATRIB(3)«0, 

ATRIB(4)*0, 

ATRIB (6) *0, 

ATRIB(6)«0, 


STOP  IF  AFTER  8  HOURS 

KTRg 

BRAICH  BACK  TO  START 
AID  COITIIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF1  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD, TRKS,  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTIIG  TIME 


I 


ATRIB(7)»0, 
ATRIB(8)-0, 
ATRIB (9) -0, 
ATRIB(10)«0; 

ACT/1 1,,,LDER; 

HAUL  GOOV; 

ACT/7, ATRIB<7); 

AWAIT (5) ,DUKPSPOT/l ; 

ACT/8, ATRIB(8) ; 

FREE,DUMPSP0T(1) ; 

ACT/9, ATRIB (9); 

ASSIGX , ATRIB ( 1 ) *0 , 
ATaiB(2)«0, 
ATRIB(3)*0, 
ATRIB (4) ■O, 
ATRIB(5)«0, 
ATRIB (6) »0, 
ATRIB (7) *0, 
ATRIB (8) -0. 
ATRIB(9)*0, 
ATRIB(11)«0; 

ACT/10, , ,TRXS; 

EKD; 

MOITR, TRACE, 0,480; 

TIM; 


TRUCK  HAULIIG  TIME 
TRUCK  DUKPIIG  TIME 
TRUCK  RETORI  TIME 


Modified  model  used  for  FT  scenario.  Attributes  7, 8,  and  9  (truck  hauling,  dumping,  and 
return  times)  were  cut  in  half  for  all  the  truckB. 

GEM, C1RPEITER, TRUCK  HAULIIG, 9/30/92,1, , ,T/«, , ,72; 

LIMITS, 5, 14,50; 

EVTRY/2, 0,0, 0,0, 0,0, 0,0, 0,1,0; 

ERRY/2, 0,0, 0,0,0, 0,0,0, 0,2,0; 

EITRY/2, 0,0, 0,0, 0,0, 0,0, 0,3,0; 

EHRY/2, 0,0, 0,0, 0,0, 0,0, 0,4,0; 

EITRY/3,0,0,0,0,0,0,0,0,0,0,5; 

EITRY/3, 0,0, 0,0,0, 0,0,0, 0,0,6; 

IETV0RK; 

RESOURCE/LOADSPOTd)  ,4 ; 

RESOURCE/DUHPSPOTd) ,5; 

START  CREATE;  CREATE  LOAD  TRAISACTIORS 

ASSIGI, ATRIB d)  -  II, 

ATRIB(2)  *  ERLKG(4. ,2.), 

ATRXBC4)  -  EXP0M(14.), 

ATRIB (5)  «•  EXP0K12.), 

ATRIB (6)  -  5, 

ATRIBC7)  •  RI0RM(22. ,3.)  *  .5, 

ATRIB(d)  ■  UIFRM(2.,8.)  *  .5, 

ATRIB(9)  >  RK0RMC18. ,3.)  •  .6; 

ACT/1,, TI0W.LE.480;  STOP  IP  AFTER  8  HOURS 

GOOI;  ELSE 

ACT/2, ATRIBC2),, START;  BRAICH  BACK  TO  START 

ACT/3;  AID  COITIIUE 

ACCUM.2,2;  ACCUMULATE  TWO  PILES 

LOAD  QUEUE(l) , , , ,BEGI;  QUEUE  OF  LOADS 

TRKS  QUEUE (2) , , , ,BEGM;  QUEUE  OF  TRUCKS 

LDER  QUEUE(3) , , , ,BEGI ;  QUEUE  OF  LOADERS 

BEGI  SELECT, ASM/SUM,,, LOAD, TRKS, LDER;  ASM  OF  LOAD, TRKS,  AID  LDER 

ACT/13; 

AH  AIT(4) .LOADSPOT/ 1 , , 1 ; 

ACT/4, ATRIB (4) ,ATRIB(11) .EQ.5,FRID;  L0ADER1  TIME 

ACT/S,ATRIB(5) ,ATRIB(11) .EQ.6;  L0ADER2  TIME 

FRID  FREE.LOADSFOT(l) ; 

GOOI; 

ACT/12,,, HAUL; 

ACT/6, ATRIB (6);  LOADER  RESTIIG  TIME 

ASSIGI,  ATRIB  d)>*0, 

ATRIB (2) “0, 

ATRIB (3) “0, 

ATRIB (4) "0, 

ATRIB(B)*0, 

ATRIB (8) "0, 


STOP  IP  AFTER  8  HOURS 
ELSE 

BRAICH  BACK  TO  START 
AID  COITIIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASH  OF  LOAD, TRKS,  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTIIG  TIME 


JSb.,- .v  yd 


ATRIB (7) -0, 
ATRIB(8)"0, 
ATRIB(9)«0. 
ATRIB ( 10) "0; 

ACT/11 , , ,LDER; 

HAUL  GO 01; 

ACT/7, ATBIB(7) ; 

AVAZT(S) ,DUMPSP0T/1; 

ACT/8, ATRIB (8) ; 

FREE , DUMPSPOT ( 1 ) ; 

ACT/9, ATRIB(9) ; 

ASSIGI, ATRIB (1)*0, 
ATRIB(2)*0, 
ATRIB(3)«0, 
ATRIB(4)-0, 
ATRIB (5) -0, 
ATRIB(6)*0, 
ATRIB(7)-0, 
ATRIB (8) “0, 
ATRIB ( 9) *0, 
ATRIB(11)*0; 

ACT/10,,, TRKS; 

EVD; 

HOVTR, TRACE, 0,480; 

FI*; 


TRUCK  HAULIIG  TIKE 
TRUCK  DUKPIIG  TINE 
TRUCK  RETURN  TINE 


/  •' 


Modified  model  used  for  FL  scenario.  Attributes  4  and  5  (loader  loading  times)  were 
multiplied  by  0.  5. 


GEI.CARPEITER, TRUCK  HAULIIG, 9/30/92, 1, . ,T/I, , ,72; 
LIMITS, 5,14,60; 

EITRT/2, 0,0,0, 0,0,0, 0,0, 0,1,0; 

EITRY/2, 0,0, 0,0,0, 0,0, 0,0, 2,0; 

EITRT/2, 0,0, 0,0, 0,0,0, 0,0, 3,0; 

EITRT/2, 0,0, 0,0, 0,0, 0,0, 0,4,0; 
EITRY/3,0,0,0,0,0,0,0.0,0,0,5; 
EITRT/3,0,0,0,0,0,0.0,0, 0,0,6; 


IETV0RK; 

RESOURCE/LOADSFuT ()  )  ,4; 

RESOURCE/ DUMPSPOT(l) ,5; 

START  CREATE; 

ASSIGI.ATRIB(l)  -  II. 

ATRIBC2)  •  ERLSGC4..2.), 
ATRIBC4)  -  EXP0K14.)  *  .6, 
ATRIB(S)  -  EXP0K12.)  *  .5, 
ATRIB(6)  •  5, 

ATRIBC7)  -  RI0RM(22. ,3.) , 
ATRIB(8)  •  UIFRM(2.,8.>, 
ATRIB(9)  -  RIORMdS. ,3.); 

ACT/1,, TI0W.LE.480; 

0001; 

ACT/2, ATRlB(2) , .START; 

ACT/3; 

ACCUN.2,2; 

LOAD  QUEUE(l) , , , ,BEGI; 

TRKS  QUEUEC2) , , , ,BEGI ; 

LDER  QUEUE(3) , , , ,BEGI; 

BEGI  SELECT, ASM/ SUK, , .LOAD, TRKS, LDER; 

ACT/13; 

AWAIT (4) .LOADSPOT/l , , 1 ; 

ACT/4, ATRIB (4) ,ATRIB(11) .EQ.S.FRID; 

ACT/5, ATRIB(S) ,ATRIB(11) .EQ.6; 

FRID  FREE.LOADSPOTd) ; 

GOOI; 

ACT/12,,, HAUL; 

ACT/6,  ATMB(6) ; 

ASSIGI,ATRIB(l)-0, 

ATRIB  (2) -0, 

ATRIB (3) *0, 

ATRIB (4)*0, 

ATRIB(6)«0, 

ATRIB(6)"0, 


CREATE  LOAD  TRAISACTIOMS 


STOP  IF  AFTER  8  HOURS 

FT -SB 

BRAICH  BACK  TO  START 
AID  COHTIIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD, TRKS,  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTIIG  TINE 


v  122 


,  C-  3?  :v-. 


/ 


i 


L 


ATRIB (7) -0, 

ATRIB(8)»0, 

ATRIB (9) »0, 

ATRIB ( 10) -0; 

ACT/11, , ,LDER; 

HAUL  GO  01; 

ACT/7, ATRIB(7) ;  TRUCK  HAULING  TIME 

AWAIT (S) ,DUH?5P0T/1; 

ACT/8, ATRIB (8);  TRUCK  DUMPING  TIKE 

PREE.DUKPSPOT(l) ; 

ACT/9, ATRIB (9) ;  TRUCK  RETURN  TIME 

.  ASSIGN , ATRIB ( 1 ) *0 , 

ATRIB (2) °0, 

ATRIB (3) "0, 

ATRIB(4)»0, 

ATRIB (5) -0, 

ATRIB(6)«0, 

ATRIB (7)*0, 

ATRIB (8) *0, 

ATRIB(9)-0, 

ATRIB(ll)*0; 

ACT/10,,, TRKS; 

END; 

MONTR, TRACE, 0,480; 

FIN; 
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Modified  model  used  for  SL  scenario.  Attribute  2  (load  interarrival  time)  was  doubled. 


GEI , CARPEMTER , TRUCK  BAULIKG , 9/30/92 , 1 , , ,  Y/I ... 72 ; 


LIMITS, 5, 14, 50; 

EITRY/2.0,0,0,0, 0,0, 0,0,0, 1.0; 

EITRY/2,0, 0,0,0, 0,0, 0,0, 0,2,0; 
EITRY/2,0,0,0,0, 0,0, 0,0,0,3,0; 

EITRY/2,0, 0,0, 0,0,0, 0,0,0, 4,0; 

EITRY/3, 0,0, 0,0, 0,0, 0,0, 0,0,5; 

EITRY/3, 0,0, 0,0, 0,0, 0,0,0, 0,6; 

IETW0RK; 

RESOURCE/LOADSPOT (1) ,4; 
RESOURCE/DUHPSPOT  Cl) ,5; 

START  CREATE; 

ASSIGI.ATRIB(l)  -  II, 

ATRIB(2)  >  ERLMG(4. ,2.)  *  2., 
ATRIBC4)  ■  EZP0K14.), 
ATRIB(5)  -  EZP0IC12.), 
ATRIB(6)  -  5, 

ATRIBC7)  -  RI0RMC22. ,3.) , 
ATRIB(B)  -  UIFRM(2. ,8.) , 
ATRIBO)  ■  RI0RM(18.  ,3.) ; 
ACT/1,. TIOI.LE. 480; 

GOOI; 

ACT/2,ATRIB(2), .START; 

ACT/3; 

ACCUM.2,2; 

LOAD  QUEUE(l) , , , ,BEGI; 

TEES  QUEUEC2) , , , ,BEGI ; 

LDER  QUE0E(3) , , , ,BEGI ; 

BEGI  SELECT, ASM/SUM, , , LOAD, TRKS, LDER; 
ACT/13; 

AHAIT(4) .LOAOSPOT/l , , 1 ; 

ACT/4, ATRIB (4) .ATRIB(ll) .EQ.5.FRID; 
ACT/5, ATRIB(S) .ATRIB(ll) .EQ.6; 

FRID  FRBE.LOAOSPOT(l) ; 

GOOI; 

ACT/12, , .HAUL; 

ACT/6, ATRIBC6); 

ASSIGI, ATRIB (1)*0, 

ATRIB (2) "0, 

ATRIB(3)«0, 

ATRIB (4) *0, 

ATRIB (5)«0, 

ATRIB(6)»0, 

ATRIB (7)*0, 
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CREATE  LOAD  TRAMS ACTIONS 


STOP  IP  AFTER  8  HOURS 
EJSE 

B RAICH  BACK  TO  START 
AID  COITIIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD, TRKS,  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTIIG  TIME 


ATRIB(8)*0, 

ATRIB (9) -0, 

ATRIB ( 10) -0; 

ACT/11, , ,LDER; 

HAUL  GOON; 

ACT/7, ATRIB (7) ;  TRUCK  HAULING  TIKE 

AHAIT(S) .DUKPSPOT/l; 

ACT/8, ATRIB(8);  TRUCK  DUMPING  TIME 

FREE.DUMPSPOT(l); 

ACT/9, ATRIB (9);  TRUCK  RETURN  TIME 

ASSIGN, ATRIB (1)«0, 

ATRIB (2) »0, 

ATRIB (3) «0, 

ATRIB(4)«0, 

ATRIB (5) «0, 

ATRIB (6) -0, 

ATRIB(7)-0, 

ATRIB (8) «0, 

ATRIB (9) -0, 

ATRIB(ll)-0; 

ACT/10, , ,TRKS; 

END; 

MONTR, TRACE, 0,480; 

FIN; 
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Modified  model  used  for  ST  scenario.  Attributes  7,  8,  and  9  (truck  hauling,  dumping,  and 
return  times)  were  multiplied  by  2  for  all  trucks. 

GEI , CARPEMTER , TRUCX  HAULIIG , 9/30/92 , 1 ,, ,Y/I ,,, 72 ; 

LIMITS, 6,14, 60; 

E*TRY/2,0, 0,0,0, 0,0,0,0, 0,1,0; 

BMTRY/2, 0.0, 0,0,0, 0,0,0, 0,2,0; 

EITRY/2,0,0,0,0,0, 0,0,0, 0,3,0; 

B*TRY/2,0, 0,0,0, 0, 0,0,0, 0,4,0; 

HTRY/3,0, 0,0,0, 0,0,0, 0,0,0, 6; 

EITRY/3,0,0, 0,0, 0,0,0, 0,0, 0,6; 

IETV0RK; 

RESOORCE/LO  ADSPOT  ( 1 ) ,  4 ; 

RESOORCE/DUMPSPOTd)  ,6; 

START  CREATE;  CREATE  LOAD  TRAMS ACTI OHS 

ASSIGI.ATRIB(l)  -  II. 

ATRIB (2)  -  ERLIG(4. ,2.) , 

ATRIB  (4)  -  EXP0K14.), 

ATRIB(5)  -  EXP0K12.), 

ATRIBC6)  •  5, 

ATRIBC7)  •  RV0RM(22. ,3.)  *  2., 

ATRIB(8)  -  0IFRM(2. ,8.)  *  2., 

ATRIB(9)  -  RIQRM(18. ,3.)  *  2.; 

ACT/1,. TIOW.LE. 480;  STOP  IF  AFTER  8  HOURS 

GOOI •  ELSE 

ACT/2, ATRIB (2),, START;  BRAICH  BACK  TO  START 

ACT/3;  AID  COMTIIUE 

ACCOM, 2, 2;  ACCUMULATE  TWO  PILES 

LOAD  QUEUE(l) , , , ,BEGK;  QUEUE  OF  LOADS 

TRKS  QUEUE(2) . , , ,BEGI ;  QUEUE  OF  TRUCKS 

LDER  QUEUEC3) , , , ,BEGV;  QUEUE  OF  LOADERS 

BEGI  SELECT, ASM/SUM,,, LOAD, TRKS, LDER;  ASM  OF  LOAD, TRKS,  AID  LDER 

ACT/13; 

AWAIT(4) .LOADSPOT/l , , 1 ; 

ACT/4, ATRIB(4) .ATRIB(ll) .EQ.6,FRXD;  L0ADER1  TIME 

ACT/5, ATRIB (6) ,ATRIB(11) .EQ.6;  L0ADER2  TIME 

FRKD  FREE , LO ADSPOT ( 1 ) ; 

GOOI; 

ACT/12,,, HAUL; 

ACT/6, ATRIBC6) ;  LOADER  RESTIIO  TIME 

ASSIGI, ATRIB (i)"0, 

ATRIB(2)*0, 

ATRIB (3) "0, 

ATRIB (4)*0, 

ATRIB ( 6) *0, 

ATRIB (6) «0, 


STOP  IF  AFTER  8  HOURS 
ELSE 

BRAICH  BACK  TO  START 
AID  COMTIIUE 
ACCUMULATE  TWO  PILES 
QUEUE  OF  LOADS 
QUEUE  OF  TRUCKS 
QUEUE  OF  LOADERS 
ASM  OF  LOAD, TRKS.  AID  LDER 


L0ADER1  TIME 
L0ADER2  TIME 


LOADER  RESTIIG  TIME 


1 


ATRIB (7)«0, 
ATRIB(8)«0, 
ATRIB (9) -0, 
ATRIB ( 10) *0; 

ACT/1 1...LDER; 

HAUL  COOK; 

ACT/7, ATRIB (7); 

AHAIT(S) ,DUMPSP0T/1 ; 

ACT/8, ATRIB (8); 

PREE,DUMPSPOT(l); 

ACT/9, ATRIB (9); 

ASSIGI, ATRIB (l)«0, 
ATRIB (2) «0, 
ATRIB(3)*0, 
ATRIB (4) -0, 
ATRIB (5) -0, 
ATRIB (6) "0, 
ATRIB(7)-0, 
ATRIB (8) -0, 
ATRIB (9) -0, 
ATRIB(ll)-0; 

ACT/10, , ,TRXS; 

EHD; 

NOITR, TRACE, 0,480; 

PH; 


TRUCE  BAULXXG  TIKE 
TRUCK  DUHPIIG  TIME 
TRUCK  RETURI  TIME 


Fortran  Code 

This  section  contains  the  Fortran  Code  used  in  this  experiment.  These  five  sets  of 
code  were  used  with  the  SLAM  network  code  to  create  the  unsorted  Proof  trace  files  for  the 
LDER  scenario.  The  Fortran  code  for  the  other  scenarios  is  the  same  except  for  different 
trace  file  names.  The  final  set  of  code  is  the  program  used  to  sort  the  trace  files. 


C  Capt  Mika  Carpenter,  G0R-93N,  6  Oct  92 

C  This  program  generates  the  unsorted  trace  file  for  the  color  only 
C  animation  of  the  Truck  Hauling  simulation.  The  program  PRFSRTC  must 
C  be  run  to  generate  the  sorted  trace  file  (LBERC.ATF) 

C 

PROGRAM  HAIR 
DIMEHSIOM  ISET(IOOOO) 

IICLUDE  ’ SLAMiDIR : PARAM . IKC  * 

COMMOH/SCOMl/ATRIBdOO),  DD(iOO),  ODL(iOO),  DTHQW,  II.  MFA, 
♦MSTOP.HCLHR,  RCRDR,  IPRRT,  RIRUH,  ISSET,  MTAPE,  SS(IOO), 

♦SSLCiOO)  .THEXT,  TROW,  XXdOO) 

COMMON  QSET( lOOOO) 

BQUIVALEHCE  (SSET(i) ,QSET(1)) 

RHSET- lOOOO 

ICRDR-5 

RPRET-6 

■TAPE-7 

■PLOT-2 

0PERC10 .PIIE- ’LDERC .TR », STATUS- 'UIKHOW ' ) 

CALL  SLAM 

STOP 

KID 

C  Subroutine  IRTLC  generates  the  loaders  and  the  trucks  and  places 
C  them  on  the  screen. 

SUBROUTIIE  IRTLC 

C0HM0I/SCQN1/ATRIB(100)  ,DD<100)  ,DDL(100)  ,DTHOW,II,MFA,MSTOP,RCLRR 
1,HCRDR,HPRHT,IKEUH, RHSET, HTAPE.SSdOO)  .SSLdOO)  ,THEXT,THOW, XXdOO) 
C0MM0R/UC0RM1/MYARRAT (5) 
n(D  -  iooo. 

HRITEd0.909)XId)  *  0.0 
DO  10  I  ■  1,4 

IF  (I  .EQ.  1)  THEM 
■RITEdO, 901)1*4 

901  FORMATdX,’ CREATE  Loader  »,I2) 

■RITE (10, 902) I *4 

902  F0RMAT(1X , 'PUCE  ',12 , •  OR  Loaderl * ) 

EMDIF 

IF  (I  .EQ.  2)  THER 
MRITE(10,90S)I*4 

905  FORMATdX, 'CREATE  Loader  ',12) 

WRITE (10, 900) 1+4 

FORMATdX , 'PUCE  »,  12 ,'  OH  Loader2 » ) 

BRDIF 


900 


WUTE(  10,907)  I 

907  F0RMAT(1X,  'CREATE  Truck  ',  12) 

VRXTE(10,908)I,X 

908  FORMAT (IX, 'PLACE  ’,12,'  01  Truck', II) 

XF(X  .IE.  4)  THES 

VRITE(10,909)I  *1.5 

909  FORMAT (IX, 'TIME  *,F9.4) 

EIDXF 

10  COITXVUE 
RETORI 
EVD 

SUBROUTIXE  OTPtJT 
WRITE  (10, 10) 

10  FORMAT(IX.'ESO') 

CLOSE (10) 

RETORI 

nro 

C  Subroutine  OMOIT  is  tbs  user  written  SLAM  truce  that  writes  the 
C  unsorted  trace  file  after  X1TLC  writes  the  statements  to  initialize 
C  the  animation. 

SOBROOTXIE  OMOIT(XT) 

COMMOI/SCOMi/ATRlB(iOO) ,00(100) ,DDL(100) .DTIOW.II.MFA.MSTOP.ICLIR 
l,ICRDR,IPRJfT,IIRUIfIISET,ITAPE,SS(100) ,SSL(100) ,TIEXT,TIOW,XX(100) 
COMMOI/OCOIH1/MYARRAT(6) 

CHARACTERS  LABL.WLBL 
IITEGER  IMUMB, ACTMUM.LOADIOM 
REAL  LOADTIME 

XFOlOtf  .EC.  0.0)  GO  TO  999 
LABL-IILBL(IDUM) 

C  Returns  current  node  label 
ACnOM«IROMB  (  XDOM) 

C  Returns  current  activity  nunber 
XF(IT  .EQ.  0)  GO  TO  100 
C  Go  to  100  if  a  node  label  is  hit 
IF (IT  .EQ.  -6)  GO  TO  999 
C  OMOIT  code  that  is  not  needed 

C  The  following  statsaents  skip  over  activities  that  are  not 
C  important  for  the  animation. 

XF(ACTIUM  .EQ.  *1)  GO  TO  999 
IF(ACTIUM  .EQ.  1)  GO  TO  999 


IFUCTIUM  .EQ.  2)  GO  TO  999 
IFUCTIUM  .EQ.  3)  GO  TO  999 
IFUCTIUM  .EQ.  7)  GO  TO  999 
IFUCTIUM  .EQ.  10)  GO  TO  999 
IFUCTBUH  .EQ.  11)  GO  TO  999 
IFUCTIUM  .EQ.  12)  GO  TO  999 
IFUCTIUM  .EQ.  13)  GO  TO  999 
WRITE(10,900)TIOW 
C  Writ*  oat  correct  time 
900  F0RMAT(1X, 'TIME  ’.F9.4) 

IF((ACTSUM  .EQ.  4)  .OR.  (ACTIUM  .EQ.  6))  THEM 
C  Loading  activites 

IF(ACTIUM  .EQ.  4)  THEM 
C  Loader  1 

LOADTIME  ■  ATUIB(4)/6. 

C  Divide  loading  tine  by  6  to  allow  for  state  changes  within 
C  the  animation 

LOADEUM  •  5 
EIDIF 

IFUCTIUM  .Eq.  6)  THEM 
C  Loader  2 

LOADTIME  *  ATRIB(S)/6. 

LOADEUM  -  6 
EIDIF 

ATRIB(li)  «  LOADEUM 
WRITE (10, 800) LOAD HUM 
C  Traveling  eapty 

800  FORMAT (IX, 'SET  ’,12,’  COLOR  GREEB’) 

WRITE (10, 805) T30W  ♦  LOADTIME 
806  FORMAT (IX, ’TIKE  »,F9.4) 

WRITE (10, 815) LOADHUH 
C  Traveling  loaded 

816  F0RMAT(1X, 'SET  ’,12,'  COLOR  RED') 

WRITE(10,806)TIQW  ♦  (3.  ♦  LOADTIME) 

C  Increment  time 

WRITE (10, 800) LOADEUM 
C  Traveling  empty 

WRITE(10,830)ATRXB(10) 

C  Track  partially  loaded 

830  F0RMAT(1X, » SET  ’,F3.0,»  COLOR  PIIK») 

WRITE (10, 805 )THOW  ♦  (4.  •  LOADTIKS) 

C  Increaent  time 

WRITE (10, 816) LOADEUM 
C  Traveling  loaded 
EIDIF 


IF(ACttUN  /iQ.  6)  THE! 

C  Hauling  activity 

LOADIUM  -  ATRIB(ll) 

WRITE (10, 800) LOAD BUM 
C  Traveling  empty 

WRITE(  10 ,856)  ATRIB(IO) 

C  Track  trailing  loadad 
856  FORMAT (IX, ’SET  ’,F3.0,»  COLOR  RED’} 

WRXTE(10,605)TI0W  ♦  ATRIB(6) 

C  Increment  tin* 

WRITE (10, 880) LOADIUM 
C  Loader  idle 

880  FORMAT (IX, 'SET  ',12,'  COLOR  WRITE') 

EMDIF 

IF(ACTIUM  .EQ.  8)  THE! 

C  Duaping  Activity 

WRITE(10,865)ATRIB(10) 

C  Track  dumping 

865  F0RMAT(1X, 'SET  '.F3.0,'  COLOR  YELLOW') 

EIDIF 

IF (ACTION  .EQ.  9)  THEI 
C  Truck  return  activity 

WRITE(10,875)ATRIB(10) 

C  Track  traveling  empty 

876  FORMATOX.’SET  *,F3.0,*  COLOR  GREET) 

EIDIF 
IT  •  1 
RETTJRI 

C  The  following  statements  skip  over  nodes  that  are  not 
C  important  for  the  aniaation. 

100  IFOABL  .EQ.  *  »)  GO  TO  999 

XFCLABL  .EQ.  ’LOAD')  GO  TO  999 
IFOABL  .EQ.  'LDER')  GO  TO  999 
IF(LABL  .EQ.  »BEGI»)  GO  TO  999 
IFOABL  .EQ.  'HAUL')  GO  TO  999 
IFOABL  .EQ.  'FRID')  GO  TO  999 
IFOABL  .EQ.  'STAR')  GO  TO  999 
WRITE(10,900)TIOW 
C  Write  oat  current  time 

XF(LABL  .EQ.  *TRKS»)  THEI 
WRITE(i0,915)ATRIB(10) 

C  Truck  idle  t 

915  FORMAT (IX, 'SET  '.F3.0,'  COLOR  WHITE') 


C  Capt  Nik*  Carpenter,  G0R-93M,  6  Oct  92 

C  This  prograa  generate*  the  unsorted  trac*  file  for  the  color  and  icon 
C  animation  of  th*  Truck  Hauling  simulation.  Th*  prograa  PRFSRTCI  aust 
C  ba  run  to  ganarat*  tha  sorted  trac*  file  (LDERCI.ATF) 

C 

PROGRAM  MAIM 
DIMEMSIOI  MSET( 10000) 

IICLUDE  ’ SLAM t DIR : PARAM . XXC  * 

CONNOM/SCONl/ATRIBClOO) ,  DD(100),  DOL(IOO).  DTMOW,  II,  MPA, 
♦MSTOP.ICLIR,  ICROR,  IPRIT,  MIRUX,  MXSET,  IT APE,  SS(IOO), 

♦SSL  (100) ,TVEZT,  TMOW,  XI (100) 

COMMOM  QSET(IOOOO) 

EQUITALEICE  (ISET(l) ,QSET(1)) 

MISBT-10000 

ICROR-S 

IPRIT-6 

■TAPE-7 

■PLOT-2 

OPEV ( 10 , PILE- ' LOERCI . TR *. STATUS- ' OlIXIOHI ' ) 

CALL  SLAM 

STOP 

EHO 

C  Subroutine  IITLC  ganarat**  tha  loader*  and  th*  truck*  and  place* 

C  thea  on  th*  screen. 


SUBROUTIIE  IITLC 

C0MM0I/SC0M1/ATRIB ( 100) ,DD(100) ,DDL<100) .DTMOW, II, NFA.MSTOP.MCLMR 
1.ICRDR, IPRIT, RIRUI,RISET,ITAPE,SS(100) , SSL (100) ,TVEXT,TIOW,XX(100) 
COHHOI/UCOIMl/MTARRAYCe) 

XX(1)  •  1000. 

WRm(10,909)XX<l)  *  0.0 
DO  10  I  j-  1,4 
IP  (I  .EQ.  1)  THE! 

WRITB(10, 901)1*4 
F0RNAT(1X, ‘CREATE  Loader  ’,12) 

WRITE (10, 902) 1*4 

FORMAT (IX, 'PLACE  ',12,'  01  Loaderl’) 

■  udif  \' 

IF  (I  .EQ.  2)  THEM 


WRITEdO, 905)1+4 
FORMAT (IX, 'CREATE 


FORMAT (IX, ’CREATE  Loader  ’,12) 

WRITE (10, 906) 1+4 

FORMAT (IX, ’PLACE  ’,12,’  01  Load*r2’) 
EIDIF 


WRITE (10, 907) I 

907  F0RHAT(1X, ’CREATE  Truck  *,I2) 

WRITECIO, 908)1. I 

90S  FOmT(lX, 'PLACE  ',12, »  01  Truck*, II) 

IF(I  .SB.  4)  THEM 
WRITECIO, 909)1  *  1.5 
909  F0MAT(1X,'TIKE  »,F9.4) 

UDIF 
10  COITIIUE 
RETURI 
HD 

SOBROUTIWE  OTPUT 
WRITECIO. 10) 

10  F0RHAT(1X,»ESD') 

CLOSE(IO) 

RETURI 

EVD 

C  Subroutine  UNbST  is  the  uaar  written  SLAM  trace  that  writes  the 
C  unsorted  trace  file  after  IBTLC  writes  the  stateaants  to  initial!*'* 

C  the  animation. 

SUBROUTIKE  OHOIT(IT) 

COmOV/SCOHl/ATRIBdOO)  .DDdOO)  .DDLdOO)  .DTWOW.II.NFA.NSTOP.WCLKR 
1,ICRDR,WPRMT,IHRUW,IISET,WTAPE,SS(100)  ,SSL(100)  ,TVEXT,TBOW,XX(1CO) 
COHMOW/UCOIH 1 ,‘HYARRAY (5) 

CBARACTERM  LABL.RHLBL 
IITEGER  SmJMB.ACTSUM.LOADSUM 
REAL  LOADTIKE 

IF (TWO W  .EQ.  0.0)  CO  TO  999 
LABL-IBLBL(IDOM) 

C  Returns  current  node  label 
ACTWUM-WmmBClDUM) 

C  Returns  current  activity  nuaber 
IF(IT  .EQ.  0)  GO  TO  100 
C  Go  to  100  if  at  a  node 

IFUT  .EQ.  -5)  GO  TO  999 
C  OKCIT  code  not  needed  for  aniatation 

C  The  following  stateaents  skip  over  activities  that  are  not 
C  important  for  the  animation. 

nCACTWOX  .EQ.  -1)  GO  TO  999 
IFCACTWOM  .EQ.  1)  GO  TO  999 


XFUCTIUM  .EQ.  2)  GO  TO  999 
XFUCTIUM  .EQ.  3)  GO  TO  999 
XFUCTIUM  .EQ.  7)  GO  TO  999 
XFUCTIUM  .EQ.  10)  GO  TO  999 
XFUCTIUM  .EQ.  11)  GO  TO  999 
XFUCTIUM  .EQ.  12)  GO  TO  999 
XFUCTIUM  .EQ.  13)  GO  TO  999 
WRITE (10, 900 )TMQW 
C  Writ*  out  currant  tin* 

900  FORHATdX, ’TIME  \F9.4) 

XFCUCTIUM  .EQ.  4)  .OR.  UCTIUM  .EQ.  6))  THEI 
XFUCTIUM  .EQ.  4)  THEI 
C  Loader  1 

LOADTIME  •  ATRIB(4)/6. 

LOADIUM  -  5 
EIDIF 

XFUCTIUM  .EQ.  5)  THEI 
C  Loader  2 

LOADTIME  ■  ATRXB(5)/6. 

LOADIUM  -  6 
EIDIF 

ATRIB(ll)  -  LOADIUM 
WRirE(10,eOO)LOADIUM 
C  Traveling  eaptj 

800  FORMATClX, ’SET  »,I2,*  COLOR  GREEI’) 

WRITE ( 10, 805 )TIOW  ♦  LOADTIME 
C  Increment  tin* 

806  FORMATClX, 'TIME  *,F9.4) 

WRITE (10, 8 14) LOADIUM 
C  Chang*  icon  to  Loaded 

814  FORMATClX, *SET  M2,*  CLASS  Loaded*) 
WRXTE(10,81S)LOAD!UM 

C  Traveling  loaded 

815  FORMAT (II, 'SET  ’,12,'  COLOR  RED*) 

WRITE (10, 805) TIOV  ♦  (3.  *  LOADTXHE) 

C  Increment  tin* 

WRITE (10, 826) LOADIUM 
C  Chang*  icon  to  Unloaded 
WRITE (10, 800) LOADIUM 
C  Traveling  eapty 

825  FORMATClX, *SET  *,X2,*  CLASS  Loader*) 

WRIiE(10,829)ATRIB(10) 

C  Chang*  truck  icon  to  1  load 
829  FORMATClX, *SET  *,F3.0,»  CLASS  Tloadl*) 

WRXTE(10,830)ATRXB(10) 
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C  Tr aval lag  partially  loaded 
830  F0RMAT(1X. ‘SET  ‘,F3.0,‘  COLOR  PIRK‘) 
WRITE (10, 805 )TEOW  ♦  (4. .  *  LOADTIKE) 

C  Incr*a*nt  tin* 

WRITS (10, 8 14) LOAD HUM 
C  Chang*  icon  to  Loaded 

WRITE (10, 815)L0AD SUM 
C  Traveling  loaded 

EVDIF 

IF(ACTBUH  .EQ.  8)  THEH 
LOADBOH  -  AVRIB(li) 

WRITE (10, 325) LOAD SUM 
C  Chang*  icon  to  Unleaded 

WRITE (10, 800 )L0ADSUM 
C  Traveling  eepty 

WRITE ( 10 , 854) ATRIB ( 10) 

C  Chang*  truck  icon  loaded 
864  FORHATdX, 'SET  ‘.F3.0,‘  CLASS  Tload2‘) 
WRITE ( 10 , 855) ATRIB ( 10) 

C  Truck  traveling  loaded 
856  F0RMAT(1X,‘SET  ‘.P3.0,‘  COLOR  RED') 

WRITE (10, 805) TWOW  ♦  ATRIB (6) 

C  Increaent  tin* 

WRITE (10, 880) LOAD HUM 
C  Loader  idle 

880  FORHATdX, ‘SET  ',12, '  COLOR  WHITE*) 

EKDIP 

IF(ACTWUH  .EQ.  8)  THE* 

WRITE (10 , 829) ATRIB(IO) 

C  Truck  partially  loaded 

WRITE ( 10 , 866) ATRIB (10) 

C  Truck  dueping 

866  FORHATdX, ‘SET  ',F3.0,‘  COLOR  YELLOW) 

WRITE(10,80S)T*OW  ♦  (ATRIB(8)/2. ) 

C  Increment  time 

WRITE(10,864)ATRIB(10) 

C  Chang*  truck  icon  unloaded 
864  FORHATdX, ‘SET  ‘,F3.0,f  CLASS  Truck') 

WRITE(10,865)ATRIBd0) 

C  Keep  truck  yellow 

EIDIF 

IFUCTFUM  .EQ.  9)  THEM 
WRITE(10,875)ATRIB(10) 

C  Truck  traveling  eapty 

878  FORHATdX, 'SET  »,F3.0,‘  COLOR  GREER') 
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SVDXF 
XT  ■  1 
RETORI 

C  Tbs  following  atataxanta  skip  0T«r  nodes  that  are  not 
C  important  for  tbs  animation. 

100  IFCLABL  .EQ.  »  ’)  GO  TO  999 

IFCLABL  .EQ.  ’LOAD’}  GO  TO  999 
IFCLABL  .EQ.  »LDER»)  GO  TO  999 
IFCLABL  .EQ.  ’BEGS')  GO  TO  999 
IFCLABL  .EQ.  ‘HAUL’)  GO  TO  999 
XFCLABL  .EQ.  ’FRIO’)  GO  TO  999 
XFCLABL  .EQ.  ’STAR’)  GO  TO  999 
WUTEC10.900)THOM 
C  Writs  out  currant  time 

XFCLABL  .EQ.  >TRKS>)  THEI 
WRITE C 10, 915) ATRX8 CIO) 

C  Truck  idle 

916  F0RMATC1X, ’SET  ’.F3.0,’  COLOR  WHITE’) 

EMDIF 
999  IT  ■  1 
RETURI 
EID 
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C  Capt  Mika  C&rponter ,  G0R-93M,  6  Oct  92 

C  This  program  generates  the  unsorted  trace  file  for  the  icon  only 
C  initiation  of  the  Truck  Hauling  simulation.  The  prograa  PRFSRTI  must 
C  be  run  to  generate  the  sorted  trace  file  (LDERI.ATF) 

C 

PROGRAM  MAIM 
DIMEMSIOH  MSET( 10000) 

'  HCLUDE  ' SLAM9DIR : PARAM . IHC ' 

COMMOI/SCOMl/ATRIBClOO) ,  DD(100),  DDL(IOO),  DTMOW,  II,  MPA, 
♦MSTQP.MCLMR,  SCilDR,  MPRMT,  HIRUH,  INSET,  MTAPE,  SS(100), 

♦SSL (100)  ,TNEXT,  THOW,  XX(100) 

COMMON  QSET(IOOOO) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

nSET- 10000 

MCRDR-5 

MPRMT-6 

MTAPE-7 

MPLOT-2 

OPEN ( 10 , FILE" *  LDERI . TR ’ .STATUS-' UNKNOWN') 

CALL  SLAM 

STOP 

END 

C  Subroutine  INTLC  generates  the  loaders  and  the  trucks  and  places 
C  then  on  the  screen. 

SUBROUTINE  INTLC 

COMMOM/SCQMl/ATRIBC  100)  ,DD(100)  .DDL(IOO)  ,DTKOW,II,MFA,MSTOP.MCLIIR 
l.MCRDR.NPRJrr.NHRUN.NNSET, MTAPE, SS(IOO)  ,SSL(100)  ,TNEXT,TNOW,H(100) 
C0MM0N/UC0HM1/MY ARRAY (5) 

XX(1)  -  1000. 

WRITE(10,909)XX(1)  *  0.0 
DO  10  I  -  1,4 

IF  (Z  .EQ.  1)  THEM 
WRITEdO, 901)1+4 

901  FORMATdX, 'CREATE  Loader  ',  12) 

WRITEdO,  902)1+4 

902  FORMATdX, 'PLACE  ',12,*  OM  Loaderl*) 

EMDIF 

IF  (I  .EQ.  2)  THEM 
WRITEdO,  90S)  1+4 

905  F0RMAT(1X, 'CREATE  Loader-  *,I2) 

WRITEdO, 906)1+4 

906  FORMATdX, 'PUCE  ',12,*  OM  Loader2') 

EMDIF 


\ 
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WRITE(10, 907)1 

907  F0RMAT(1X,' CREATE  Truck  ',12) 

WRITE (10, 908) I, I 

908  FORMAT (IX , ' PLACE  »,I2,'  01  Truck *,I1) 
XF(X  .IE.  4)  THEM 

WRITE (10. 909) I  *  l.S 

909  FORMAT (IX, 'TIME  \F9.4) 

EIDXF 

10  COITXIUE 
RETORI 

hd 

SUBROUTXIE  OTPUT 
WRITE  (10, 10) 

10  FORMAT (IX, *EID*) 
a0SE(10) 

RETURI 

HD 


C  Sub rout in*  UMOIT  is  ths  ussr  srittsn  SLAM  trie*  that  vrltss  tbs 
C  unsort ad  traca  fils  aftsr  IITLC  vritss  ths  stataaents  to  initializa 
C  ths  animation. 

SOBROUTIIE  UMOIT (IT) 

C0NM0I/SC0M1/ATRXB ( 100) ,DD(100) .DDL(IOO) ,DTIOW,II,MFA,MSTOP,ICLIR 
l,ICRDR.IPRIT,IIRUI,nSET,ITAPE,SSdOO)  ,SSL(100)  ,TIEXT,TBOW,XX(100) 
C0MM0I/UC0IM1/MYARRAY (5) 

CHARACTER*4  LABL.IILBL 
XITEGER  IIUMB , ACTIUM , LOADIUM 
REAL  LOADTIME 

IF(TIOW  .EQ.  0.0)  CO  TO  999 
LABL-IILBL(IDUM) 

C  Rstums  currant  nods  labsl 
ACTIUM-IIUMB(IDUM) 

C  Rstums  currant  activity  nuabsr 
ZF(IT  .EQ.  0)  00  TO  999 
C  If  at  nods  ratum 

2F(IT  .EP.  -5)  00  TO  999 
C  UMOIT  cods  not  nssdsd 

C  Ths  following  statsasnts  skip  ovsr  activitiss  that  ara  not 
C  important  for  ths  aninatlon. 

IF (ACTIUM  .EQ.  -1)  00  TO  999 
IF (ACTIUM  .EQ.  1)  00  TO  999 


140 


/ 


■  -  ■  •»"  •  r 


IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 
IFCACTNUH  .EQ. 


2)  GO  TO  999 

3)  GO  TO  999 
7)  GO  TO  999 

9)  GO  TO  999 

10)  GO  TO  999 

11)  GO  TO  999 

12)  GO  TO  999 

13)  GO  TO  999 


WRITECl0,900)TN0W 


C  Writ*  out  currant  tima 
900  FORMATClX, ’TIKE  ’,F9.4) 

IF((ACTirUM  .EQ.  4)  .OR.  (ACTVUM  .EQ.  6))  THEN 
C  If  loadar  activity 

IFCACTNUH  .EQ.  4)  THEN 
C  Loadar  1 


LOADTIHE  -  ATRIB(4)/6. 

LOADNUM  •  6 
END  IF 

IFCACTNUH  .EQ.  6)  THEN 
C  Loadar  2 

LOADTIME  -  ATRIBCS)/6. 

LOADNUM  -  6 
ENDIF 

ATRIBCll)  -  LOADNUM 
WRITE CIO, 80S )TNOV  ♦  LOADTIME 
C  Incraaent  tina 
805  FORMAT C IX, ’TIME  *,F9.4) 

WRITE C 10, 8 IS) LOADNUM 
C  Changa  loadar  icon  to  Loadad 
816  FORMATClX, ’SET  ’,12,’  a ASS  Loadad’) 

WRITEC10,80S)TNOW  ♦  C3.  *  LOADTIME) 

C  Incraaant  tina 


WRITE  C 10 , 82S ) LOADNUM 
C  Changa  loadar  icon  to  Unloadad 
826  FORMATClX, ’SET  ’,12,’  CUSS  Loadar’) 

WRITEC10,830)ATRIBC10) 

C  Changa  truck  icon  partially  loadad 
830  FORMATClX, ’SET  ’,F3.0,’  Class  Tloadl’) 

WRITE C 10, 80S )TNOW  ♦  C4.  •  LOADTIME) 

C  Incrasant  tins 

WRITE (10, 8 IS) LOADNUM 
C  Changa  loadar  icon  to  Loadad 
1IDXF 

IFCACTNUH  .EQ.  8)  THEN 
C  Hauling  activity 
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LOADIUM  -  ATRIB(ll) 

WUTEdO ,  855)  ATRIB  ( 10) 

C  Chang*  truck  icon  Loaded 
855  F0RMAT(1X, 'SET  \F3.0,’  CLASS  Tload2’) 
WRITE ( 10 , 825) LOAD NUM 
C  Chang*  loader  icon  to  Unloaded 
EVDIF 

1F(ACTSUM  .EQ.  8)  THEM 
C  Dumping  activity 

WRITE(10,830)ATRIB(10) 

C  Chang*  truck  partially  loaded 

VR1TE( 10 , 805)TVOW  ♦  (ATRIB (8) /2.) 

C  Increaent  time 

1IRITE(  10 , 865)  ATRIB  (10) 

C  Change  truck  to  empty 
865  F0RHAT(1X, 'SET  ',F3.0,»  CLASS  Truck’) 
EMDIF 
999  IT  -  1 
RETURI 
EVD 
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C  Capt  Mike  Carpenter,  G0R-93M,  6  Oct  92 

C  This  program  generates  the  unsortrd  trace  file  for  the  movement  and 
C  color  and  the  movement,  icon,  and  color 

C  animations  of  the  Truck  Hauling  simulation.  The  program  PRFSRTKC 
C  must  be  run  to  generate  the  sorted  trace  files 
C  (LDERMC . ATF  and  LDERMCI.ATF) 

C 

PROGRAM  MAIM 
DIMEMSIOM  NSET( 10000) 

INCLUDE  ' SLAMlDIR : PARAM .ISC’ 

COMMOH/SCOM1/ATRIB(100) ,  DD(100),  DDL (100) ,  DTNOW,  II,  MFA, 
♦MSTOP.NCLNR,  NCRDR,  MPRNT,  NNRUH,  NNSET,  ITAPE,  SS(100), 

♦SSL(IOO)  ,TNEXT,  TNOW,  n(100) 

COMMON  QSET(IOOOO) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

INSET** 10000 
NCRDR-5 i 
NPRNT-6 
NTAPE-7 
MPLOT-2 

OPEN ( 10 .FILE- *  LDERMC . TR  * , STATUS- » UNKNOWN  * ) 

CALL  SLAM 

STOP 

END 

C  Subroutine  INTLC  generates  the  loaders  and  the  trucks  and  places 
C  them  on  the  screen. 

SUBROUTINE  INTLC 

COMMON/SCOMl/ATRIBClOO)  ,DD<100)  ,DDL(100)  ,DTNOW,II,MFA,MSTOP.NCLMR 
1 , NCRDR  .NPRNT  .NNRUN ,  NNSET ,  NTAPE  ,SS  (100)  ,SSL(100)  .TIEXT  ,TN0W ,  XX  (100) 
COMMON  AJCOHM 1 /MY ARRAY ( 5 ) 

II  •  6 

xxd)  -  iooo. 

WRITE(10,909)XX(1)  *  0.0 
DO  10  I  -  1,4 
IF  (I  .EQ.  1)  THEN 
II  -  II  ♦  1 
WRITE (10, 90 1)1+4 

901  F0RMAT(1X,» CREATE  Loader  M2) 

WRITE (10, 902) 1+4 

902  F0RMAT(1X, 'PLACE  »,I2,»  ON  Getloadl') 

WRITEdO,  903)11 

903  F0RMAT(1X, 'CREATE  Load  ',13) 

WRITE(10, 904)11 


-i'  w  1 WT******* 1 1 1'  !r^y^??y!R!WW  i.l*’Ft'^tP»Ll  fPl.UI  .1 1 ' 


904  F0RMAT(1X , 'PLACE  M3,*  01  Loadpath’) 

EIDIF 

IF  (I  .EQ.  2)  THEN 
WRITE (10, 905) 1+4 

905  F0RMAT(1X, ’CREATE  Loadar  M2) 

WRITE (10, 906) 1+4 

906  F0RMAT(1X, ’PLACE  M2,*  0V  Gatload2’) 

EIDIF 

WRITE (10, 907) I 

907  FORMAT (IX, ’CREATE  Truck  ’,12) 

WRITE(10,908)I 

908  FORMAT ( IX , *  PLACE  ’,12,*  01  Truckload’) 

IF(I  .IE.  4)  THE! 

WRITE (10, 909) I  *  1.5 

909  FORMAT (IX, ’TIME  *,F9.4) 

EIDIF 

10  COITIIUE 
RETURI 
EID 

•  I 

I 

SUBROUTINE  OTPUT 
WRITE (10, 10) 

10  FORMAT(IX.’EID’) 

CLOSE(IO) 

RETURI 

EID  { 

■  ■ 

C  Subrout in •  UMOIT  is  tho  uaar  writtan  SLAM  traco  that  writas  tha 
C  unaortad  tract  f ila  altar  IITLC  writas  tha  stataaants  to  initializa 
C  tha  animation. 

SUBROUTINE  UMOIT(XT) 

COMMOI/SCOM1/ ATRIB ( 100) ,DD(100) ,DDL(100) ,DTIOW,II,MFA.MSTOP.ICLIR 
1 ,ICRDR,IPRIT , MIRUI , IISET , ITAPE , SS ( 100) , SSL ( 100) , TIEXT ,TNOW , XX ( 100) 
COMMOI/UCOIM1/NYARRAY (5) 

CHARACTER+4  LABL.INLBL 
INTEGER  HUMB.ACTIUM.LOADIUM 
REAL  LOADTIME 

XF(TIOW  .EQ.  0.0)  GO  TO  999 
LABL-IILBL(IDUM) 

ACTIUM-IIUMB(IDUM) 

IF(IT  .EQ.  0)  GO  TO  100 
IF(IT  .EQ.  -5)  GO  TO  999 

C  Tha  following  stataaanta  skip  owar  actlvitiaa  that  ara  not 
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C  important  for  the  animation. 

IFUCTSUM  .EQ.  -1)  GO  TO  999 
IFUCTSUM  .EQ.  1)  GO  TO  999 
IFUCTSUM  .EQ.  2)  GO  TO  999 
IFUCTSUM  .EQ.  3)  GO  TO  999 
IFUCTSUM  .EQ.  7)  GO  TO  999 
XFUCTSUM  .EQ.  10)  GO  TO  999 
XFUCTSUM  .EQ.  11)  GO  TO  999 
XFUCTSUM  .EQ.  12)  GO  TO  999 
XFUCTSUM  .EQ.  13)  GO  TO  999 
WRITE(10,900)TSDW 
C  Writ*  out  currant  time 
900  FORMAT (IX,  'TIME  \F9.4) 

IF((ACTHUM  .EQ.  4)  .OR.  (ACTSUM  .EQ.  5))  THEN 
C  Loading  activity 

IFUCTSUM  .EQ.  4)  THEN 
C  Loader  1 

LOADTIKE  -  ATRIB(4)/6. 

LOADSUM  -  S 
WRITE (10, 790) LOADSUM 
C  Send  Loader  to  get  load 
790  FORMAT (IX, ’PLACE  ’,X2,»  OS  Getloadl’) 

ESDIF 

IFUCTSUM  .EQ.  6)  THEM 
C  Loader  2 

LOADTIME  -  ATRIB(5)/6. 

LOADSUM  >6 
WRITE(10,795)LOADSUM 
C  Sand  Loader  to  get  load 

79S  FORMAT( IX, ’PLACE  ’,12,’  OR  Getload2’) 

ESDIF 

ATRIB(ll)  -  LOADSUM 
WRITE(10,799)LOADSUM 
C  Traveling  empty 

799  F0RMAT(1X, ’SET  ’,12,’  COLOR  GREES’) 

WRITE ( 10 , 800) LOADSUM , LOADTIME 

C  Sat  travel  time 

800  F0RMAT(1X, ’SET  »,X2,»  TRAVEL  \F9. 4) 

WRITE (10, 80S )TSOV  4  LOADTIME 

C  Increment  time 

806  FORMAT (IX, ’TIME  ’.F9.4) 

WRITE(10,810)ATRIB(1)  -  1 
C  Loader  picks  up  load 
810  FORMAT( IX, ’DESTROY  ’,F3.0) 
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WRITEdO, 815)L0ADHUM 
C  Chang*  loadar  to  Load*d 
816  FORMAT (IX,’ SET  M2,'  CLASS  Load*d') 

WRITEdO, 816)L0ADHUM 
C  Traveling  loaded 

816  FORMAT ( IX SET  »,I2,»  COLOR  RED*) 

WRITE (10, 820) LO ADIUM 
C  Plac*  load  in  truck 

820  FORMAT (IX, 'PLACE  ',12,'  01  Dunpload’) 

WRITE ( 10 , 800) LO ADIUM , LOADTIME 
C  S*t  travel  tin* 

WRITE(10,805)TMOW  ♦  (3.  •  LOADTIME) 

C  Increment  tin* 

WRITEdO,  b25)L0ADIUM 
C  Chang*  loader  icon  to  unloaded 
826  FORMAT (IX, 'SET  ',12,'  CUSS  Loader') 

WRITE(10,799)LOADSUM 
C  Traveling  empty 

WRITE(10,830)ATRIB(10) 

C  Chang*  truck  icon  to  partially  loaded 

830  FORMAT (IX, 'SET  '.F3.0,'  CLASS  Tloadl') 

WRITE(10,831)ATRIB(10) 

C  Truck  partially  loaded 

831  F0RMAT(1X, 'SET  *,F3.0,'  COLOR  PIIK») 

WRITEdO,  US)  LOADIUM 

C  Send  loader  to  get  new  load 
836  F0RMAT(1X,'PUCE  »,I2,»  01  Getnewload') 

WRITE(10 . 800)LOADHUM .LOADTIME 
C  Set  loader  travel  tine 

WRITE(10,806)TIOW  ♦  (4.  *  LOADTIME) 

C  Increment  tin* 

WRITEdO  ,810)  ATRIB  ( 1 ) 

C  Loader  pick*  up  aecond  load 
WRITEdO, 816)L0ADIUM 
C  Chang*  loader  to  loaded 
WRITE(10,816)L0ADRUM 
C  Traveling  loaded 

WRITE (10, 820) LO ADIUM 
C  Plac*  load  in  truck 

WRITEdO , 800)  L01DNUM ,  LOADTIME 
C  Set  travel  tine 
EIDIF 

IF(ACTliUM  .EQ.  6)  THEM 
C  Hauling  activity 

LO  ADIUM  •  ATRIBdl) 
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WRITE (10, 825) LOADWUH 
j  Change  loader  icon  to  unloaded 
WRITE (10, 799) LOADHUM 
C  Traveling  empty 

IF (LOAD HUM  .EQ.  5)  THEI 
WRITE (10,840) LOAD HUM 
C  Return  loader  1 

840  FORMAT (IX,’ PLACE  ',12,'  OH  Returnl’) 

EHDIF 

IF(LOADHUM  .EQ.  6)  THEM 
WRITE(10,84S)L0ADHUM 
C  Return  loader  2 

845  F0RMAT(1X, ’PLACE  M2,'  OH  Retum2‘) 

EHDIF 

WRITE ( 10 , 850) LOADHUM , ATRIB ( S ) 

C  Set  travel  time 

850  FORMAT (IX,  ’SET  M2,'  TRAVEL  *,F3.0) 

WRITE ( 10 , 855 ) ATRIB ( 1 0) 

C  Change  truck  to  loaded 

855  FORMAT (IX. 'SET  ’.F3.0.’  CLASS  Tload2*) 
WRITE(10,856)ATRIB(10) 

C  Truck  traveling  loaded 

856  FORMAT (IX, ’SET  ’,F3.0,’  COLOR  RED’) 
WRITE(10,860)ATRTB(10) 

C  Start  truck  hauling 

860  FORHATdX. ’PLACE  ’.F3.0,’  OH  Truckduop’) 
WRITE(10,875)ATRIB(10) ,ATPIB(7) 

C  Set  truck  travel  time 

WRITE(10,805)TNOW  ♦  ATRIB (6) 

C  Increment  time 

WRITE (10,880) LOADHUM 
C  Loader  idle 

880  FORHATdX, ’SET  ’.12,’  COLOR  WHITE’)  ( 
EHDIF 

IF(ACTHUM  .EQ.  8)  THEH  | 

C  Dumping  activity 

xxd)  -  xxd)  ♦  i 

WRITE(10,861)ATRIB(10) 

C  Start  truck  dumping 

861  F0RMAT(1X, ’PLACE  \F3.0,’  OH  Dumping’) 
WRITE ( 10 , 862) ATRIB(IO) ,ATRIB(8)/2. 

C  Set  travel  time 

862  FORMAT (IX, ’SET  ’.F3.0,’  TRAVEL  ’,F9.4) 
WRITE (10, 805) THOW  +  (ATRIB(8)/2. ) 

C  Increment  time 
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WRITE( 10 , 865)  ATRIB(IO) 

C  Change  truck  to  empty 
865  F0RMAT(1X, ' SET  ’.F3.0,'  CLASS  Truck’) 
WRXTE(10,863)ATRIB(10) 

C  Truck  duaping 

863  FORMAT (IX, 'SET  ’,F3.0,'  COLOR  YELLOW) 

VRITE(10,866)XX(1) 

C  Place  two  loads  in  truck 

888  FORMAT ( IX ,  ’ CREATE  Twoloads  ’ , F6 . 0) 

WRITE(10, 867)11(1) 

C  Start  two  loads  duaping 

867  FORMAT ( IX , ' PLACE  \FS.O,*  OR  Loadduap') 

EVDIF 

IFUCTEUM  .EQ.  9)  THEM 
C  Truck  return  activity 

WRITE(10,869)ATRIB(10) 

C  Truck  traveling  empty 

869  FORMAT (IX.' SET  \F3.0,'  COLOR  GREER') 
URITE(10,870)ATRIB(10) 

C  Start  truck  returning 

870  FORMAT (IX, 'PLACE  '.F3.0,'  OR  Truckret') 
WRITE(10,875)ATRIB(10) ,ATR1B(9) 

C  Set  truck  travel  tine 

875  F0RMAT(1X,'SET  '.F3.0,'  TRAVEL  \F9.4> 

ERDIF 
IT  -  1 
RETURI 

C  The  following  statements  skip  over  nodes  that  are  not 
C  iaportant  for  the  animation. 

100  XF(LABL  .EQ.  ’  ’)  GO  TO  999 

XF(LABL  .BQ.  'LOAD')  GO  TO  999 
.  IF(LABL  .EQ.  »U> ER*>  GO  TO  999 
‘  IF(LABL  .EQ.  'BEGR')  GO  TO  999 
IF(LABL  .EQ.  'HAUL’)  GO  TO  999 
XF(LABL  .EQ.  'FRED')  GO  TO  999 
WRITE (10, 900) TROW 
C  Write  out  current  tiae 

XF(LABL  .EQ.  'STAR')  TBER 
IX  ■  II  ♦  1 
WRITE(10,905)XI 
C  Create  new  load 
905  F0RMAT(1X. 'CREATE  Load  »,I3) 

WRITE (10, 9 10) II 


.^i'<ij.fyypf,iui.j!j|ii(|jjli).i,j; 


Vi*'-’--  V'  ‘--^v.--v  -<  ■■  -  • 


. . .  lll'JIJIJIp 


C  Sand  load  for  pickup 

910  FORMAT ( IX ,  ’ PLACE  *,13,*  OH  Loadpath’) 
EHDIF 

IFCLABL  .EQ.  ’TRKS’)  THEM 
WRITE (10 , 91 1 ) ATRIB (10) 

C  Truck  idle 

911  FORMAT (IX,’ SET  ’,F3.0,’  COLGR  WHITE’) 
WRITE ( 10 , 91S) ATRIB ( 1 0) 

C  Placa  truck  in  loading  queue 
916  FORMAT (IX, ’PUCE  »,F3.0,’  OH  Truckload’) 

WRITE( 10 , 920) ATRIB ( 10) 

C  Sat  truck  travel  tins 

920  FORMAT (IX, ’SET  \F3.0,*  TRAVEL  2’) 

EHDIF 
999  IT  ■  1 
RETURN 
END 
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C  Capt  Nika  Carp ant nr,  G0R-93N,  6  Oct  92 

C  This  program  generates  the  unsort ad  tr&ca  file  for  the  movement  only 
C  and  the  movement  and  icon  animations  of  the  Truck  Hauling  simulation. 
C  The  program  PRFSRTMI  must  be  run  to  generate  the  sorted  trace  files 
C  (LDERM.ATF  and  LOERMI.ATF) 

C 

PROGRAM  MAIM 
DIMENSION  ISET(IOOGO) 

INCLUDE  *  SLAMtDIR :  PARAM .  INC  ’ 

COMMOR/SCOMl/ATRIBClOO) ,  DD(100),  DDL(IOO),  DTHOW,  II.  MPA. 
♦MSTOP.ICLKR,  ICROR.  NPRNT,  INRUN,  UNSET,  NT APE,  SS(IOO), 

♦SSL(IOO) ,TNEXT,  TNOW,  XX(100) 

COMMON  QSET(IOOOO) 

EQUIVALENCE  (NSET(l) ,QSET(i» 

INSET-10000 

NCRDR-5 

NPRNT-6 

■TAPE-7 

NPLOT-2 

OPEN (10, FILE- 1 LDERMI.TR* .STATUS- 'UNKNOWN') 

CALL  SLAM 

STOP 

END 

C  Subroutine  IMTLC  generates  the  loadors  and  the  trucks  and  places 
C  them  on  the  screen. 

SUBROUTINE  INTLC 

C0MM0N/SC0M1/ATRIB ( 100) ,DD(100) ,DDL(100) ,DTNOW,II,KFA,MSTOP,NCLHR 
1,ICRDR,NPRNT,NNRUN, INSET, XTAPE,SS (100) ,SSL(100),TNEXT,TNOW,ZX(100) 
COMMOI/UCONMi/NY ARRAY (S) 
ii  -  e 

n(i)  -  looo. 

WRITE(10,909)ZX(1)  *  0.0 
DO  10  I  •  1.4 
IF  (I  .EQ.  1)  THEN 
II  -  II  4  1 
WRITE (10, 901) 1+4 

901  FQRMAT( IX,* CREATE  Loader  ’,12) 

WRITE (10, 902) 1+4 

902  FORMATUX, 'PLACE  ',12,'  ON  Getloadl') 

WRITEdO, 903)11 

903  F0RMAT(1X, 'CREATE  Load  »,  13) 

WRITEdO, 904)  II 

904  F0RMATC1X, 'PLACE  ',13/  ON  Loadpath') 


;»all!.»y-1,  |i.ii,»u.ji^iiiiii!m>iiBi.p'gii".').i>.»->,w>ii-.11.1"-  ’ mm?  t^Bggpw^!Wrl,,iwsiiM 


EHDIF 

IF  (I  .EQ.  2)  THEM 
WRITE(10, 905)1*4 

905  FORMAT (IX,* CREATE  Loader  *,I2) 

«RITE(10,906)I+4 

906  FORMAT ( IX , *  PLACE  \I2,*  ON  Getload2’) 

EIDIF 

WRITE (10, 907) I 

907  FORMATC IX,' CREATE  Truck  M2) 

WRITE (10, 908) I 

908  FORMAT (IX, * PLACE  »,I2,»  ON  Truckload*) 

IF(I  .ME.  4)  THEN 

WRITE (10, 90S) I  *1.6 

909  F0RMAT(1X, ’TIME  *,F9.4) 

EMDIF 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  OTPUT 
WRITE (10, 10) 

10  F0RMAT(1X,’END’) 

CLOSE (10) 

RETURN 

END 

C  Subroutine  UHONT  ia  the  user  written  SLAM  trace  that  writes  the 
C  unsorted  trace  file  after  INTLC  writes  the  atatenents  to  initialize 
C  the  animation. 

SUBROUTINE  UMONT(IT) 

COMMON/SCOM 1/ATRIB ( 100)  ,DD(100)  ,DDL(100)  ,DTNOW,II,MFA,MSTOP,HCLHR 
1,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100)  .SSL (100)  ,TNEXT,TNOW,XX(100) 
COMMON/UCONMl/MYARRAY (5 ) 

CHARACTERS  LABL.NNLBL 
INTEGER  NNUKB , ACTNUM , LO ADNUM 
REAL  LOADTIKE 

IF(TNOW  .EQ.  0.0)  GO  TO  999 
LABL-NVLBL(IDUM) 

C  Returns  current  node  label 
ACTNUM-NNUMB(IDUM) 

C  Returns  current  activity  number 
IF (IT  .EQ.  0)  GO  TO  100 
C  Go  to  100  if  at  node 

IF(IT  .EQ.  -6)  GO  TO  999 
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C  UMONT  cods  that  is  not  needed 


C  The  following  statements  skip  over  activities  that  are  not 
C  important  for  the  animation. 

IPUCTMUM  .EQ.  -1)  GO  TO  999 
IF (ACTION  .EQ.  1)  GO  TO  999 
IFUCTMUM  .EQ.  2)  GO  TO  999 
IFUCTMUM  .EQ.  3)  GO  TO  999 
IFUCTMUM  .EQ.  7)  GO  TO  999 
IFUCTMUM  .EQ.  10)  GO  TO  999 
IFUCTMUM  .EQ.  11)  GO  TO  999 
IFUCTMUM  .EQ.  12)  GO  TO  999 
IFUCTMUM  .EQ.  13)  GO  TO  999 
WRITE (10, 900) TSOW 
C  Write  out  current  time 
900  FORMAT (IX, ’TIME  ’.F9.4) 

IF((ACTIUH  .EQ.  4)  .OR.  (ACTNUM  .EQ.  6))  THEM 
C  Loader  activity 

IFUCTMUM  .EQ.  4)  THEM 
C  Loader  1 

L0ADTIH2  -  ATRIB(4)/6. 

LOADMUM  »  5 
WRITE(10,790)LOAOMUM 
C  Send  loader  to  get  load 
790  F0RMAT(1X, 'PLACE  ',12,*  OM  Getloadl’) 

EMDIF 

IFUCTMUM  .EQ.  5)  THEM 
C  Loader  2 

LO AOTIME  •  ATRIB(5)/6. 

LOADMUM  •  6 
WRITE ( 10 , 7S5) LOADMUM 
C  Send  loader  to  get  load 
796  FORMAT( IX, 'PLACE  ',12,'  OH  Getload2>) 

EMDIF 

ATRIB(ll)  -  LOADMUM 
WRITEdO ,  800)  LOADMUM ,  LOADTIME 
C  Set  travel  time 

800  FORMAT (IX, 'SET  ’,12,'  TRAVEL  ’,F9.4) 

WRITE (10, 805) TIOW  ♦  LOADTIME 
C  Increment  time 
805  F0RMAT(1X,'TIME  »,F9.4) 

WRITE(10, 810) ATRIB(l)  -  1 
C  Place  load  on  loader 

810  FORMATdX, 'DESTROY  ’,F3.0)  ' 


\ 


fkWIPV-T- 


WRITE(10,81S)LOADNUM 
C  Change  loader  to  loaded 
815  FORMAT (IX, ‘SET  ‘ ,12, *  CLASS  Loaded’) 

WRITE (10, 820) LOADNUH 
C  Place  load  in  truck 

820  F0RMAT(1X, ‘PLACE  ’,12,’  ON  Dumpload’) 

WRITE ( 10 , 800) LOADNUM , LOADTIME 
C  Set  travel  time 

WRITE(10,805)TNOW  +  (3.  *  LOADTIME) 

C  Increment  time 

WRITE  ( 10 , 825)  LOADNUM 
C  Change  loader  to  unloaded 
825  FORMAT (IX, ’SET  ’,12,’  CLASS  Loader’) 
WRITE (10,830) ATRIB ( 10) 

C  Change  truck  to  partially  loaded 
830  F0RMAT(1X, ’SET  ’,F3.0,’  CLASS  Tloadl’) 

WRITE  (10,835)  LOADNUM 
C  Send  loader  to  get  new  load 
835  FORMAT (IX, ’PLACE  ’,12,’  ON  Getnewload’) 
WRITE ( 10 , 800)LOADNUM , LOADTIKE 
C  Set  travel  time 

WRITE (10, 805 )TN0W  ♦  (4.  *  LOADTIME) 

C  Increment  time 

WRITE( 10,810) ATRIB (1) 

C  Place  load  on  loader 

WRITE(10,815)L0ADNUM 
C  Change  loader  to  loaded 

WRITE ( 10 , 820) LOADNUM 
C  Place  load  in  truck 

WRITE (10,800) LOADNUM .LOADTIME 
C  Set  travel  time 
ENDIF 

IFUCTNUM  .EQ.  6)  THEN 
C  Hauling  activity 

LOADNUM  -  ATRIB (11) 

WRITE( 10 , 825) LOADNUM 
C  Change  loader  to  unloaded 

IF(LOADNUM  .EQ.  5)  THEN 
WRITE( 10 , 840)L0ADNUM 
C  Send  loader  1  to  rest  spot 
840  FORMAT  (IX, ’PUCE  ’,12,’  ON  Returnl’) 

ENDIF 

IF (LOADNUM  .EQ.  6)  THEN 
WRITE (10, 845) LOADNUM 
C  Send  loader  2  to  rest  spot 
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045  FORMATClX, ’PLACE  ’,12,’  01  !Uturn2*) 

HOXF 

VRXTE(1O,8SO)LOADKUN,ATRXB(0) 

C  Sat  travel  tiaa 

850  FORMATClX, ’SET  *,12,*  TRAVEL  *,F3.0) 

WRITEC10,855)ATRIBC10) 

C  Chang*  truck  to  loadad 
8S8  FORMATClX, * SET  *,F3.0,*  aASS  Tload2*) 
WRITE CIO, 860) ATRIBC 10) 

C  Start  truck  hauling 

860  FORMATClX, ’PLACE  \F3.0,*  OM  Truckdump’) 
WRITEC 10, 876) ATRIBC 10) .ATRIBC7) 

C  Sat  traval  tin# 

EXDXF 

XFCACTWUM  .EQ.  8)  THEM 
C  Dumping  activity 

XXCl)  •  MCI)  ♦  1 
WRITEC10,861)ATRXB(10) 

C  Start  truck  duaping 

861  FORMATClX, ’PUCE  \F3.0,’  01  Duaping’) 
WRITE CIO, 862) ATRIBC 10) ,ATRIB(8)/2. 

C  Sat  truck  traval  tiaa 

•62  FORMATClX, ’SET  ’,F3.0,*  TRAVEL  ’,F9.4> 

WRITEC10,805)TWOW  ♦  CATRXBC8)/2. ) 

C  Xncraant  tiaa 

WRITE C 10, 865) ATRIBC 10) 

C  Changa  truck  to  unloaded 
866  FORMATClX, ’SET  ’.F3.0,*  CUSS  Truck’) 

WRITEC 10,666) XXCl) 

C  FI  xa  tvo  loads  in  truck 

866  FORMATClX, ’CREATE  Tvoloada  *,F5.0) 

WRITEC 10,867) XXCl) 

C  Start  t*o  load#  dumping 
•67  FORMATClX. 'PUCE  ’.F5.0,’  01  Loadduap') 

EVDXF 

XFCACTWUM  .EQ.  9)  THEM 
C  Saturn  activity 

WRITEC 10, 670) ATRIBC 10) 

C  Start  truck  returning 

•70  FORMATClX, ’PUCE  ’.F3.0,’  01  Truckrat’) 

WRITEC10,876)ATRIBC10) .ATRIBC9) 

C  Sat  truck  traval  tiaa 
•76  FORMATClX, ’SET  ’,F3.0,»  TRAVEL  \F9.4) 

HOXF 

n  •  i 


154 


/ 


RETURN 


C  The  following  statements  skip  over  nodes  that  are  not 
C  important  for  the  animation. 


100 


IFCLABL  .EQ.  *  ’) 

IF(LAEL  .EQ.  'LOAD') 
IFCLABL  .EQ.  ’ LPER * ) 
IF(LABL  .EQ.  ’BEGN’) 
IFCLABL  .EQ.  'HAUL’) 
IFCLABL  .EQ.  'FUND') 
VRITECl0.900)TNOW 
IFCLABL  .EQ.  'STAR’) 
II  -  II  ♦  1 
WRITEC10,905)II 


GO  TO 
GO  TO 
GO  TO 
GO  TO 
GO  TO 
GO  TO 

THEN 


999 

999 

999 

999 

999 

999 


C  Create  nos  load 
905  FORMATC IX, ’CREATE  Load  ’,13) 
WRITEC10, 910)11 


C  Send  load  to  load  queue 

910  FORKATC IX, ’PLACE  ’.13,’  ON  Loadpath’) 

END  IF 

IFCLABL  .EQ.  ’TRKS’)  THEN 
WIITEC10,91S)ATRIBC10) 

C  Place  truck  in  loading  queue 
915  F0RMATC1X, ’PLACE  ’,F3.0,»  ON  Truckload’) 
HRITECl0,920)ATRIBCl0) 

C  Set  truck  travel  tine 

920  FORMATC IX, ’SET  ’.F3.0,’  TRAVEL  2') 

ENDIF 
999  IT  •  1 
RETURN 
END 
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This  is  a  sample  of  the  programs  used  to  created  sorted  Proof  trace  files.  Only  one  is 
listed  since  the  only  differences  in  the  programs  are  the  input  and  output  file  names.  This 
particular  program  took  one  unsorted  trace  file  and  created  two  sorted  Proof  trace  files, 
one  for  M  and  one  for  Ml.  Thus,  the  trace  files  for  M  and  MI  were  exactly  the  same.  Also, 
the  trace  files  for  MC  and  MCI  were  the  same.  Differences  in  the  animations  were  made  by 
changing  the  animation  layout  files.  All  of  the  object  classes  for  M  and  MC  were  changed 
to  simple  icons.  For  example,  when  the  trace  file  changed  the  object  class  to  loaded,  no 
change  was  seed  on  the  screen  for  M  and  MC  because  all  the  object  classes  were  the  simple 
icons.  (See  Figure  2  on  page  17).  Therefore,  each  scenario  had  to  be  run  only  five  times 
instead  of  seven  to  created  the  seven  Proof  trace  files. 

PROGRAM  SORT 

INTEGER  I , J , K , NUMLINES , LINE , FLAG , SORTING , READLINE 
CHARACTER'S  TIME 
CHARACTER* 80  LINES (3500) 

REAL  TIMES, TIMEARRY(3S00, 3), TEMPI, TEMP2.TEMP3 
OPEN (UNIT-10 , FILE" * ldermi . tr ’ .STATUS-* UNKNOWN') 

OPEN (UNIT- 1 1 , FILE- * ldermi . atf » , STATUS- ’ UNKNOWN  * ) 

OPEN (UNIT- 12 , FILE- ’ lderm . atf  * , STATUS- *  UNKNOWN » ) 

1-0 

NUMLINES-0 

LINE-0 

FUG-0 

10  READ  ( 10 ,  FMT- 100 ,  END-200)  TIME 

100  FORMAT(AS) 

LINE  -  LINE  *  1 
NUMLINES  -  HUMLINES  4  1 
IF(TIME  .EQ.  *  TIME*)  THEN 

IF  (FUG  .Eq.  0)  THEN 
FUG  -  1 
NUMLINES  -  0 
GO  TO  10 
ENDIF 
I  ■  I  ♦  1 

TIMEARRY(I.l)  -  TIMES 
TIMEARRY(I ,2)  -  NUMLINES  -  1 
TIMEARRY(I ,3)  -  LINE  -  NUMLINES 
BACKSPACE  10 

READ(10,FMT-101 ,ERR-195)TIMES 

101  F0RMAT(6X,F9.4) 

NUMLINES-0 

ENDIF 
GO  TO  10 

195  WRITE(ll.llO) 

1 10  FORMAT ( IX , *  BACKSPACE  ERROR ' ) 

200  I  ■  I  4  1 
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TIMEARRY(I.l)  ■  TIMES 
TIMEARRY (I, 2)  ■  NUHLINES  -  1 
TIMEARRY (I, 3)  -  LINE  -  NUMLINES 
C 

REWIND  10 
HEADLINE  -  1 

16  READ(10,FMT-116,END-300)LINES(READLINE) 

116  F0RMAT(A60) 

READLINE  -  HEADLINE  ♦  1 
GO  TO  16 
300  CONTINUE 
C 

DO  26  RUNTHRU  «  1,1 
DO  20  SORTING  •  1,1-1 

IF (TIMEARRY (SORTING, 1)  .GT.  TIMEARRY (SORTING+1 , 1) )  THEN 
TEMPI  ■  TIMEARRY (SORTING+1 , 1) 

TIMEARRY (SORTING+1, 1)  «  TIMEARRY(SORTING.l) 

TIMEARRY (SORTING, 1)  ■  TEMPI 
TEMP2  »  TIMEARRY (SORTING+1, 2) 

TIMEARRY(S0RTING+1 ,2)  -  TIMEARRY (SORTING, 2) 

TIMEARRY (SORTING ,2)  -  TEMP2 
TEMP3  ■  TIMEARRY (SORTING+1 , 3) 

TIMEARRY(S0RTING+1,3)  •  TIMEARRY (SORTING, 3) 

TIMEARRY (SORTING ,3)  »  TEHP3 
ENDIF 

20  CONTINUE 

26  CONTINUE 
C 

DO  36  J  ■  1,1 
LINE  •  TIMEARRY(J,3) 

NUMLINES  -  TIMEARRY(J,2) 

DO  30  K  -  LINE, LINE  +  NUMLINES 
WRITE(11 ,116)LINES(K) 

WRITE(12, 116) LINES (K) 

30  CONTINUE 

36  CONTINUE 

WRITE  (11, 126) 

WRITE(12,126) 

126  FORMAT (IX, ‘END’) 

C 

CLOSE (10) 

CLOSE(ll) 

CLOSE (12) 

END 
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